Skip to content

固体物理期末复习资料 ​

依据课件:第 1 章(1.1–1.12)、第 2 章(2.1–2.6)、第 3 章(3.1–3.5)、第 4 章(5.1–5.10 及恒定磁场、太阳能电池专题)、第 5 章(6.1–6.3),并参考「总结」文件夹中的各章总结。

编号对照:课件文件名 3.x 对应第 3 章,5.x 对应第 4 章,6.x 对应第 5 章。

标记约定:「补充」表示课件之外、为把推导补完整而增加的内容;「待核对」表示课件表述不清或前后不一致之处。

符号约定:ℏ=h/2π\hbar=h/2\pi;kBk_{\mathrm B} 为玻尔兹曼常量;qq 在晶格振动中为波矢,在电子论中为电荷,按上下文区分;VV 在热容部分为晶体体积,VcV_c 为原胞体积,按上下文区分。

公式说明约定:为保持简洁,推导链中的中间等式不逐条重复单位说明;凡标有「重要公式」的结论式,均在其后给出符号含义、单位与适用条件;其余中间式新引入的符号在推导文字中说明。

目录 ​

第 0 章 绪论 ​

研究对象:固体物理研究的不是单个原子,而是大量原子组成固体后表现出的集体性质;重点在于用原子结构、电子结构解释固体的宏观性质,并据此控制、改善这些性质。

固体的分类:

  • 晶体:至少在微米量级范围内原子排列具有周期性,即长程有序。单晶体整块由一个晶粒组成;多晶体由许多单晶(晶粒)通过晶界结合而成。
  • 非晶体:不具有长程周期性,但保留短程序(近邻原子数目、种类、键长、键角与晶体相近)。
  • 准晶体:有长程取向性,但没有长程平移对称性;1984 年 Shechtman 在快速冷却的 AlMn 合金中发现,其电子衍射图具有五重对称性。

晶体的宏观特性:长程有序性(周期性、平移对称性)、自限性、各向异性、解理性、均匀性、对称性、固定熔点。晶体的宏观特性由内部结构的周期性决定,是微观特性的反映。

检测手段:X 射线衍射、电子衍射。晶体有尖锐的衍射峰;非晶的衍射峰宽缓,俗称「馒头峰」。

易错点:

  • 晶体外形(晶面大小、形状)受生长条件影响,不是晶体品种的特征因素;晶面角守恒才是品种的特征因素。
  • 多晶体整体上不表现出单晶那样的各向异性。
  • 准晶不是「有长程序的晶体」:它有取向序,但无平移周期性。

0.1 易错点与常考题型 ​

易错点

  1. 把「长程有序」等同于「所有方向都周期性」;准晶有长程取向序但无平移周期序。
  2. 把多晶的整体性质当作单晶的各向异性。
  3. 把晶体外形当作晶体品种的判据;品种判据是晶面角守恒和衍射花样。

常考题型

  1. 判断晶体、非晶体、准晶体,并说明判据。
  2. 列举晶体的宏观特性,解释其微观来源。
  3. 说明单晶、多晶、液晶、非晶的区别。

第 1 章 晶体结构 ​

1.1 晶体的宏观特性与分类 ​

长程有序性(周期性):晶体可看作平移单元在空间重复堆砌而成。一维为等间距点列;二维为平行四边形或六边形网格;三维为平行六面体网格。

自限性:晶体自发形成封闭多面体的能力。表现为有规则的外形。同一品种的晶体,对应晶面间的夹角恒定,称为晶面角守恒。

各向异性:不同方向上原子排列不同,物理性质不同。表现为解理性、双折射等。解理面是面间结合较弱的晶面;解理面通常比较光滑、平整。

固定熔点:晶体熔化时温度出现平台;非晶没有固定熔点。

分类与判据:

类型原子排列衍射特征
晶体(单晶/多晶)长程有序尖锐衍射峰
非晶体短程有序、长程无序宽缓的「馒头峰」
准晶体有长程取向序,无长程平移序五重等非晶体学对称衍射

1.2 密堆积与致密度 ​

配位数(CN):一个原子周围最近邻原子的数目。排列越紧密,配位数越大,结合能越低。

最密堆积:全同原子球的最紧密堆积。同层内每个球与 6 个球相切;上下两层各有 3 个球相切,故 CN =12=12。

两种最密堆积:

  • 立方密积(面心立方,fcc):密排面按 ABCABC⋯ABCABC\cdots 顺序堆垛。典型晶体:Cu、Ag、Au、Al、Pb、Ni、γ\gamma-Fe。
  • 六角密积(hcp):密排面按 ABAB⋯ABAB\cdots 顺序堆垛。典型晶体:Be、Mg、Zn、Ti、Cd、Zr。

致密度(堆积因数):

ρ=vV\rho=\frac{v}{V}

其中 ρ\rho 无量纲;vv 为晶胞内原子球总体积,单位 m3\mathrm{m^3};VV 为晶胞体积,单位 m3\mathrm{m^3}。

晶胞中原子数:

N=ni+12nf+14ne+18ncN=n_i+\frac{1}{2}n_f+\frac{1}{4}n_e+\frac{1}{8}n_c

其中 ni,nf,ne,ncn_i,n_f,n_e,n_c 分别为内部、面上、棱上、顶角上的原子数;贡献因子依次为 1,12,14,181,\tfrac12,\tfrac14,\tfrac18。若讨论二维平面,则内部为 11、边界为 12\tfrac12、顶角为 14\tfrac14。

重要公式:面心立方致密度

N=8×18+6×12=4,4R=2 aN=8\times\frac18+6\times\frac12=4,\qquad 4R=\sqrt2\,a

ρ=N⋅43πR3a3=4⋅43π(24a)3a3=2 π6≈0.7405\rho=\frac{N\cdot\frac43\pi R^3}{a^3} =\frac{4\cdot\frac43\pi\left(\frac{\sqrt2}{4}a\right)^3}{a^3} =\frac{\sqrt2\,\pi}{6}\approx 0.7405

其中 aa 为晶格常数(m\mathrm{m}),RR 为原子球半径(m\mathrm{m})。适用条件:硬球模型、原子相切。

重要公式:体心立方致密度

N=8×18+1=2,4R=3 aN=8\times\frac18+1=2,\qquad 4R=\sqrt3\,a

ρ=2⋅43π(34a)3a3=3 π8≈0.680\rho=\frac{2\cdot\frac43\pi\left(\frac{\sqrt3}{4}a\right)^3}{a^3} =\frac{\sqrt3\,\pi}{8}\approx 0.680

适用条件:硬球模型、体对角线方向相切。

重要公式:六角密积致密度

N=2,a=2R,c=223 aN=2,\qquad a=2R,\qquad c=2\sqrt{\frac23}\,a

ρ=2⋅43πR332a2c=2 π6≈0.7405\rho=\frac{2\cdot\frac43\pi R^3}{\frac{\sqrt3}{2}a^2c} =\frac{\sqrt2\,\pi}{6}\approx 0.7405

其中 cc 为六角原胞高(m\mathrm{m}),aa 为底面边长(m\mathrm{m})。适用条件:理想 hcp,c/a=8/3≈1.633c/a=\sqrt{8/3}\approx1.633。

理想 hcp 的 c/ac/a 证明:在密堆积中,底层三个球与上层一个球构成正四面体,球心间距为 aa。上层球心在底面上的投影位于底面正三角形中心,到底面顶点的距离为 23⋅32a=33a\tfrac23\cdot\tfrac{\sqrt3}{2}a=\tfrac{\sqrt3}{3}a。故半高

(c2)2=a2−(33a)2=23a2⟹c=223 a=83 a\left(\frac c2\right)^2=a^2-\left(\frac{\sqrt3}{3}a\right)^2=\frac23a^2 \quad\Longrightarrow\quad c=2\sqrt{\frac23}\,a=\sqrt{\frac83}\,a

其中 cc 的单位为 m\mathrm m,aa 的单位为 m\mathrm m。适用条件:理想 hcp、硬球相切。

三种立方格子的比较:

项目简立方 sc体心立方 bcc面心立方 fcc
晶胞体积a3a^3a3a^3a3a^3
每晶胞格点数112244
原胞体积a3a^3a3/2a^3/2a3/4a^3/4
最近邻数66881212
最近邻距离aa32a\frac{\sqrt3}{2}a22a\frac{\sqrt2}{2}a
次近邻数12126666
次近邻距离2 a\sqrt2\,aaaaa

易错点:

  • fcc 晶胞含 44 个格点,原胞只含 11 个;bcc 晶胞含 22 个,原胞含 11 个。
  • 计算致密度时,RR 与 aa 的关系必须由「最近邻原子相切」确定:sc 为 a=2Ra=2R,bcc 为 3a=4R\sqrt3a=4R,fcc 为 2a=4R\sqrt2a=4R。
  • hcp 的 c/ac/a 不是任意值;理想密堆积要求 c/a=8/3c/a=\sqrt{8/3}。

1.3 空间点阵、基元、原胞与晶胞 ​

布拉伐空间点阵:晶体内部结构可看成相同的「点子」在空间作规则、周期性、无限分布;这些点子称为结点(格点)。格点的位矢为

Rl=l1a1+l2a2+l3a3\mathbf R_l=l_1\mathbf a_1+l_2\mathbf a_2+l_3\mathbf a_3

其中 l1,l2,l3l_1,l_2,l_3 为整数;a1,a2,a3\mathbf a_1,\mathbf a_2,\mathbf a_3 为三个基矢,单位 m\mathrm{m};Rl\mathbf R_l 为格点位矢(正格矢),单位 m\mathrm{m}。

晶体与基元的关系:

晶体 = 点阵(晶格)+ 基元

基元是晶体的最小重复结构单元。简单晶格的基元是一个位置等价的原子;复式晶格的基元含两个以上位置不等价的原子。

原胞与晶胞:

  • 原胞:晶格中体积最小的重复单元,每个原胞含且只含一个格点;基矢记为 a1,a2,a3\mathbf a_1,\mathbf a_2,\mathbf a_3,体积

    Ω=a1⋅(a2×a3)\Omega=\mathbf a_1\cdot(\mathbf a_2\times\mathbf a_3)

    单位 m3\mathrm{m^3}。
  • 晶胞(单胞):为反映晶格对称性而选取的较大周期单元,可能含多个格点;基矢记为 a,b,c\mathbf a,\mathbf b,\mathbf c。

四种基本原胞基矢:

简立方:

a1=ai,a2=aj,a3=ak,Ω=a3\mathbf a_1=a\mathbf i,\quad \mathbf a_2=a\mathbf j,\quad \mathbf a_3=a\mathbf k,\qquad \Omega=a^3

面心立方:

a1=a2(j+k),a2=a2(k+i),a3=a2(i+j),Ω=a34\mathbf a_1=\frac a2(\mathbf j+\mathbf k),\quad \mathbf a_2=\frac a2(\mathbf k+\mathbf i),\quad \mathbf a_3=\frac a2(\mathbf i+\mathbf j),\qquad \Omega=\frac{a^3}{4}

体心立方:

a1=a2(−i+j+k),a2=a2(i−j+k),a3=a2(i+j−k),Ω=a32\mathbf a_1=\frac a2(-\mathbf i+\mathbf j+\mathbf k),\quad \mathbf a_2=\frac a2(\mathbf i-\mathbf j+\mathbf k),\quad \mathbf a_3=\frac a2(\mathbf i+\mathbf j-\mathbf k),\qquad \Omega=\frac{a^3}{2}

六角密积(第一种取法):

a1=12a i+32a j,a2=−a i,a3=c k,Ω=32a2c\mathbf a_1=\frac12a\,\mathbf i+\frac{\sqrt3}{2}a\,\mathbf j,\quad \mathbf a_2=-a\,\mathbf i,\quad \mathbf a_3=c\,\mathbf k,\qquad \Omega=\frac{\sqrt3}{2}a^2c

其中 aa 为晶格常数或六角底面边长(m\mathrm{m}),cc 为六角原胞高(m\mathrm{m}),i,j,k\mathbf i,\mathbf j,\mathbf k 为直角坐标单位矢量。

晶胞与原胞的关系:

Vcell=n VprimV_{\mathrm{cell}}=n\,V_{\mathrm{prim}}

其中 VcellV_{\mathrm{cell}} 为晶胞体积,VprimV_{\mathrm{prim}} 为原胞体积(单位均为 m3\mathrm{m^3});nn 为一个晶胞所含的格点数。例如 fcc:Vcell=a3V_{\mathrm{cell}}=a^3,Vprim=a3/4V_{\mathrm{prim}}=a^3/4。

易错点:

  • 原胞不一定取立方体;fcc、bcc 的原胞是平行六面体,不是立方体。
  • 「晶胞含几个原子」与「原胞含几个格点」不同:金刚石晶胞含 88 个原子,但原胞含 22 个原子。
  • 同一晶格的基矢取法不唯一,但原胞体积唯一。

1.4 典型晶体结构 ​

结构型单胞中原子数原子位置(以 aa 为单位)最近邻距离配位数
fcc(Cu)44(000)(000),(01212)(0\tfrac12\tfrac12),(12012)(\tfrac12 0\tfrac12),(12120)(\tfrac12\tfrac12 0)22a\frac{\sqrt2}{2}a1212
bcc(W)22(000)(000),(121212)(\tfrac12\tfrac12\tfrac12)32a\frac{\sqrt3}{2}a88
CsClCs$^+$ 11,Cl$^-$ 11(000)(000),(121212)(\tfrac12\tfrac12\tfrac12)32a\frac{\sqrt3}{2}a88
NaClNa$^+$ 44,Cl$^-$ 44Na 占 fcc 位置;Cl 占 (1200)(\tfrac12 00) 等位置a2\frac a266
金刚石88fcc 位置加 (141414)(\tfrac14\tfrac14\tfrac14) 等四组位置34a\frac{\sqrt3}{4}a44
ZnS(闪锌矿)Zn 44,S 44两套 fcc 沿体对角线相对位移 14\tfrac1434a\frac{\sqrt3}{4}a44

结构描述要点:

  • NaCl:Na$^+$ 和 Cl$^-$ 各组成一套 fcc 子晶格,两套沿棱方向相对位移 a/2a/2 套构。
  • CsCl:Cs$^+$ 和 Cl$^-$ 各组成一套简立方子晶格,沿体对角线相对位移 12\tfrac12 套构;CsCl 本身不是 bcc 点阵,而是复式晶格。
  • ZnS(闪锌矿):两套 fcc 子晶格沿体对角线相对位移 14\tfrac14 套构。
  • 金刚石:具有 ZnS 结构,由同种碳原子组成;每个晶胞含 88 个碳原子,每个碳原子 sp3sp^3 杂化,与 44 个近邻形成共价键,键角 109∘28′109^\circ28'。
  • 钙钛矿 ABO$_3$:A 为二价或一价金属,B 为四价或五价金属;BO$_3$ 形成氧八面体基团,是钙钛矿结构的标志。典型:CaTiO$_3$、BaTiO$_3$、PbZrO$_3$、LiNbO$_3$、LiTaO$_3$。

易错点:

  • CsCl 的配位数是 88,但不能说「CsCl 是体心立方点阵」;它是两套简立方套构的复式晶格。
  • NaCl 的配位数是 66,最近邻距离是 a/2a/2,不是 2a/2\sqrt2a/2。
  • 金刚石与 ZnS 的原子位置相同,但金刚石是单质,ZnS 是两种原子。

1.5 晶列指数与晶面指数 ​

1.5.1 晶列指数 ​

晶列:所有格点分列在一系列相互平行的直线上,这些直线称为晶列;一族平行晶列定义一个晶向。

原胞晶列指数:设格点位矢

Rl=l1a1+l2a2+l3a3\mathbf R_l=l_1\mathbf a_1+l_2\mathbf a_2+l_3\mathbf a_3

若 l1,l2,l3l_1,l_2,l_3 为互质整数,则用 [l1l2l3][l_1l_2l_3] 表示该晶列方向。负指数在数字上方加横线,如 [1ˉ11][\bar1 11]。

晶胞晶列指数:以 a,b,c\mathbf a,\mathbf b,\mathbf c 为基矢,格点位矢

R=ma+nb+pc\mathbf R=m\mathbf a+n\mathbf b+p\mathbf c

若 m,n,pm,n,p 为互质整数,则用 [mnp][mnp] 表示晶列方向。

求法:选定原点 → 写出待求晶列上某格点的位矢 → 化为互质整数 → 填入方括号,负数戴帽。等效晶向用 ⟨mnp⟩\langle mnp\rangle 表示。

1.5.2 晶面指数(密勒指数) ​

晶面族:通过任一格点可作一族平行、等距、格点分布相同的晶面;通过一个格点可作无限多族晶面。

密勒指数求法(以晶胞基矢 a,b,c\mathbf a,\mathbf b,\mathbf c 为坐标轴):

  1. 求晶面在三个坐标轴上的截距(以 a,b,ca,b,c 为单位);
  2. 取截距的倒数;
  3. 化为互质整数 h,k,lh,k,l;
  4. 写成 (hkl)(hkl)。

用截距表示,即

1r:1s:1t=h:k:l\frac1r:\frac1s:\frac1t=h:k:l

其中 r,s,tr,s,t 为晶面在三个坐标轴上的截距(以 a,b,ca,b,c 为单位的无量纲数);h,k,lh,k,l 为互质整数。

第二种方法(法线方向余弦法):晶面法线与三个基矢的夹角余弦之比等于三个晶面指数之比:

cos⁡α1:cos⁡α2:cos⁡α3=h1:h2:h3\cos\alpha_1:\cos\alpha_2:\cos\alpha_3=h_1:h_2:h_3

其中 αi\alpha_i 为晶面法线与第 ii 个基矢的夹角(rad\mathrm{rad});hih_i 为晶面指数(无量纲)。适用条件:以原胞基矢定义的晶面指数。

若晶面与某坐标轴平行,则该方向截距为无穷大,对应指数为 00。

示例:

  • 截距为 4a,b,c4a,b,c 的晶面:倒数为 14,1,1\tfrac14,1,1,密勒指数为 (144)(144)。
  • 截距为 −3a,−b,2c-3a,-b,2c 的晶面:倒数为 −13,−1,12-\tfrac13,-1,\tfrac12,严格按定义应化为 (2ˉ6ˉ3)(\bar2\bar6 3);课件此处写作 (263)(263),待核对(课件可能省略了负号)。
  • 立方体中 AEGAEG 面:截距为 a,b,ca,b,c,指数 (111)(111);ABCDABCD 面:与 cc 轴平行,指数 (001)(001);DIHGDIHG 面:截距为 a2,b,∞\tfrac a2,b,\infty,指数 (120)(120)。

晶面指数与面间距:密勒指数越简单(如 (100)(100)、(110)(110)),晶面间距越大,原子面密度越大,晶面间结合力越弱,越容易解理。

重要公式:面间距

dh1h2h3=2π∣Gh1h2h3∣=2π∣h1b1+h2b2+h3b3∣d_{h_1h_2h_3}=\frac{2\pi}{\lvert \mathbf G_{h_1h_2h_3}\rvert} =\frac{2\pi}{\lvert h_1\mathbf b_1+h_2\mathbf b_2+h_3\mathbf b_3\rvert}

其中 dd 为面间距(m\mathrm{m});Gh1h2h3\mathbf G_{h_1h_2h_3} 为对应晶面的倒格矢(m−1\mathrm{m^{-1}});h1,h2,h3h_1,h_2,h_3 为晶面指数。适用条件:以原胞基矢为坐标轴定义的晶面指数;若用晶胞基矢定义密勒指数,则应使用对应的倒格子基矢。

易错点:

  • 晶列指数用方括号 [ ][\,],晶面指数用圆括号 ( )(\,),晶面族用大括号 { }\{\,\},等效晶向用尖括号 ⟨ ⟩\langle\,\rangle。
  • 求密勒指数时先取截距的倒数,再化整;不要直接把截距化整。
  • 负指数戴帽表示,不要写成负号在括号外。
  • fcc、bcc 的晶面指数可用晶胞基矢(密勒指数)或原胞基矢定义,两者数值不同,不能混用面间距公式。

1.6 倒格子 ​

定义:由正格子原胞基矢 a1,a2,a3\mathbf a_1,\mathbf a_2,\mathbf a_3 定义三个倒格子基矢

b1=2πa2×a3a1⋅(a2×a3),b2=2πa3×a1a1⋅(a2×a3),b3=2πa1×a2a1⋅(a2×a3)\mathbf b_1=2\pi\frac{\mathbf a_2\times\mathbf a_3}{\mathbf a_1\cdot(\mathbf a_2\times\mathbf a_3)},\quad \mathbf b_2=2\pi\frac{\mathbf a_3\times\mathbf a_1}{\mathbf a_1\cdot(\mathbf a_2\times\mathbf a_3)},\quad \mathbf b_3=2\pi\frac{\mathbf a_1\times\mathbf a_2}{\mathbf a_1\cdot(\mathbf a_2\times\mathbf a_3)}

其中 bi\mathbf b_i 的单位为 m−1\mathrm{m^{-1}};分母 Ω=a1⋅(a2×a3)\Omega=\mathbf a_1\cdot(\mathbf a_2\times\mathbf a_3) 为正格子原胞体积(m3\mathrm{m^3})。

重要公式:正倒格子基矢的正交关系

ai⋅bj=2πδij={2π,i=j0,i≠j\mathbf a_i\cdot\mathbf b_j=2\pi\delta_{ij} =\begin{cases}2\pi,&i=j\\[2pt]0,&i\ne j\end{cases}

其中 δij\delta_{ij} 为克罗内克符号,无量纲。该式是倒格子定义的核心,也是所有倒格子计算的基础。

倒格矢:

Kh=h1b1+h2b2+h3b3\mathbf K_h=h_1\mathbf b_1+h_2\mathbf b_2+h_3\mathbf b_3

其中 h1,h2,h3h_1,h_2,h_3 为整数;Kh\mathbf K_h 的单位为 m−1\mathrm{m^{-1}}。

重要公式:正格矢与倒格矢的正交关系

Rl⋅Kh=2πμ,μ∈Z\mathbf R_l\cdot\mathbf K_h=2\pi\mu,\qquad \mu\in\mathbb Z

其中 Rl\mathbf R_l 为正格矢(m\mathrm{m}),Kh\mathbf K_h 为倒格矢(m−1\mathrm{m^{-1}}),μ\mu 为整数。

倒格子原胞体积:

Ω∗=b1⋅(b2×b3)=(2π)3Ω\Omega^*=\mathbf b_1\cdot(\mathbf b_2\times\mathbf b_3)=\frac{(2\pi)^3}{\Omega}

其中 Ω∗\Omega^* 的单位为 m−3\mathrm{m^{-3}},Ω\Omega 的单位为 m3\mathrm{m^3}。

证明:由 bi\mathbf b_i 的定义,b1⋅(b2×b3)\mathbf b_1\cdot(\mathbf b_2\times\mathbf b_3) 的分子为 (2π)3(2\pi)^3 乘以三个叉乘的混合积,分母为 Ω3\Omega^3;利用 (a2×a3)⋅[(a3×a1)×(a1×a2)]=Ω3(\mathbf a_2\times\mathbf a_3)\cdot[(\mathbf a_3\times\mathbf a_1)\times(\mathbf a_1\times\mathbf a_2)]=\Omega^3,即得 Ω∗=(2π)3/Ω\Omega^*=(2\pi)^3/\Omega。

倒格矢的方向与面间距:

重要公式:倒格矢垂直于对应晶面

Gh1h2h3⊥(h1h2h3)\mathbf G_{h_1h_2h_3}\perp (h_1h_2h_3)

即倒格矢垂直于对应的 (h1h2h3)(h_1h_2h_3) 晶面族。其中 Gh1h2h3=h1b1+h2b2+h3b3\mathbf G_{h_1h_2h_3}=h_1\mathbf b_1+h_2\mathbf b_2+h_3\mathbf b_3,单位 m−1\mathrm{m^{-1}}。

证明:在晶面族 (h1h2h3)(h_1h_2h_3) 中取离原点最近的晶面,它与三个坐标轴的交点分别为 a1/h1\mathbf a_1/h_1、a2/h2\mathbf a_2/h_2、a3/h3\mathbf a_3/h_3。晶面内两个矢量可取

CA→=a1h1−a3h3,CB→=a2h2−a3h3\overrightarrow{CA}=\frac{\mathbf a_1}{h_1}-\frac{\mathbf a_3}{h_3},\qquad \overrightarrow{CB}=\frac{\mathbf a_2}{h_2}-\frac{\mathbf a_3}{h_3}

于是

Gh1h2h3⋅CA→=(h1b1+h2b2+h3b3)⋅(a1h1−a3h3)=2π−2π=0\mathbf G_{h_1h_2h_3}\cdot\overrightarrow{CA} =(h_1\mathbf b_1+h_2\mathbf b_2+h_3\mathbf b_3)\cdot\left(\frac{\mathbf a_1}{h_1}-\frac{\mathbf a_3}{h_3}\right) =2\pi-2\pi=0

同理 Gh1h2h3⋅CB→=0\mathbf G_{h_1h_2h_3}\cdot\overrightarrow{CB}=0。因此 Gh1h2h3\mathbf G_{h_1h_2h_3} 垂直于该晶面。

重要公式:面间距与倒格矢长度

dh1h2h3=2π∣Gh1h2h3∣d_{h_1h_2h_3}=\frac{2\pi}{\lvert\mathbf G_{h_1h_2h_3}\rvert}

其中 dd 的单位为 m\mathrm{m},G\mathbf G 的单位为 m−1\mathrm{m^{-1}}。适用条件:晶面指数按原胞基矢定义。

证明:原点到最近晶面的距离等于 a1/h1\mathbf a_1/h_1 在晶面法线方向 G/∣G∣\mathbf G/\lvert\mathbf G\rvert 上的投影:

d=a1h1⋅G∣G∣=1h1∣G∣a1⋅(h1b1+h2b2+h3b3)=2π∣G∣d=\frac{\mathbf a_1}{h_1}\cdot\frac{\mathbf G}{\lvert\mathbf G\rvert} =\frac{1}{h_1\lvert\mathbf G\rvert}\mathbf a_1\cdot(h_1\mathbf b_1+h_2\mathbf b_2+h_3\mathbf b_3) =\frac{2\pi}{\lvert\mathbf G\rvert}

常用结果:

  • 简立方:bi=2πaii\mathbf b_i=\frac{2\pi}{a}\mathbf i_i,倒格子仍为简立方。
  • 面心立方:倒格子为体心立方;反之,体心立方的倒格子为面心立方。
  • 倒格子原胞体积 Ω∗=(2π)3/Ω\Omega^*=(2\pi)^3/\Omega;正、倒格子互为倒逆。

例:证明面心立方的倒格子为体心立方

面心立方原胞基矢为

a1=a2(j+k),a2=a2(k+i),a3=a2(i+j),Ω=a34\mathbf a_1=\frac a2(\mathbf j+\mathbf k),\quad \mathbf a_2=\frac a2(\mathbf k+\mathbf i),\quad \mathbf a_3=\frac a2(\mathbf i+\mathbf j),\qquad \Omega=\frac{a^3}{4}

由倒格子基矢定义:

b1=2πa2×a3Ω=2πa(−i+j+k)\mathbf b_1=2\pi\frac{\mathbf a_2\times\mathbf a_3}{\Omega} =\frac{2\pi}{a}(-\mathbf i+\mathbf j+\mathbf k)

b2=2πa3×a1Ω=2πa(i−j+k)\mathbf b_2=2\pi\frac{\mathbf a_3\times\mathbf a_1}{\Omega} =\frac{2\pi}{a}(\mathbf i-\mathbf j+\mathbf k)

b3=2πa1×a2Ω=2πa(i+j−k)\mathbf b_3=2\pi\frac{\mathbf a_1\times\mathbf a_2}{\Omega} =\frac{2\pi}{a}(\mathbf i+\mathbf j-\mathbf k)

这三个矢量正是边长为 4π/a4\pi/a 的体心立方格子的三个原胞基矢,故面心立方的倒格子为体心立方。其中 i,j,k\mathbf i,\mathbf j,\mathbf k 为直角坐标单位矢量;aa 为晶格常数(m\mathrm m)。同理可证体心立方的倒格子为面心立方。

易错点:

  • 倒格子基矢定义中的 2π2\pi 不能漏;若采用晶体学惯例(不带 2π2\pi),面间距公式相应改为 d=1/∣G∣d=1/\lvert\mathbf G\rvert。
  • b1\mathbf b_1 不一定平行于 a1\mathbf a_1;它垂直于 a2,a3\mathbf a_2,\mathbf a_3 所在平面。
  • 面间距公式 d=2π/∣G∣d=2\pi/\lvert\mathbf G\rvert 只适用于与倒格矢对应的那一族晶面;不要对任意方向乱用。

1.7 晶体的对称性 ​

对称性:经过某种动作后,晶体能够自身重合的特性。对称操作:使晶体自身重合的动作。对称素:对称操作所依赖的几何要素(点、线、面)。

四种基本点对称操作:

  1. 旋转对称 CnC_n:绕固定轴转 2π/n2\pi/n 后重合,该轴为 nn 度旋转轴。绕 x1x_1 轴转 θ\theta 的矩阵为

    A=(1000cos⁡θ−sin⁡θ0sin⁡θcos⁡θ)A=\begin{pmatrix} 1&0&0\\ 0&\cos\theta&-\sin\theta\\ 0&\sin\theta&\cos\theta \end{pmatrix}

    其中 θ\theta 的单位为 rad\mathrm{rad};AA 无量纲。
  2. 中心反演 ii:(x1,x2,x3)→(−x1,−x2,−x3)(x_1,x_2,x_3)\to(-x_1,-x_2,-x_3),

    A=(−1000−1000−1),det⁡A=−1A=\begin{pmatrix}-1&0&0\\0&-1&0\\0&0&-1\end{pmatrix},\qquad \det A=-1

  3. 镜象 mm:以 x3=0x_3=0 面为对称面,

    A=(10001000−1),det⁡A=−1A=\begin{pmatrix}1&0&0\\0&1&0\\0&0&-1\end{pmatrix},\qquad \det A=-1

  4. 旋转—反演 nˉ\bar n:绕轴转 2π/n2\pi/n 后再作中心反演。

重要公式:晶体允许的旋转轴

2cos⁡θ=m,m=0,±1,±22\cos\theta=m,\qquad m=0,\pm1,\pm2

θ=2πn,n=1,2,3,4,6\theta=\frac{2\pi}{n},\qquad n=1,2,3,4,6

其中 θ\theta 为旋转角(rad\mathrm{rad}),mm 为整数,nn 为旋转轴度数。适用条件:三维晶体的平移周期性。

证明:在晶面内取相邻格点 A,O,BA,O,B,∣OA∣=∣OB∣\lvert OA\rvert=\lvert OB\rvert。绕 OO 转 θ\theta 后,A→A′A\to A',B→B′B\to B'。由于晶格周期性,A′B′A'B' 也必须是格点列上的整数倍:

∣A′B′∣=2cos⁡θ ∣AO∣=m∣AO∣\lvert A'B'\rvert=2\cos\theta\,\lvert AO\rvert=m\lvert AO\rvert

故 cos⁡θ=m/2\cos\theta=m/2。由 −1≤cos⁡θ≤1-1\le\cos\theta\le1 得 m=0,±1,±2m=0,\pm1,\pm2,对应 θ=π/3,π/2,2π/3,π,2π\theta=\pi/3,\pi/2,2\pi/3,\pi,2\pi,即 n=6,4,3,2,1n=6,4,3,2,1。五重轴对应 cos⁡θ=±12\cos\theta=\pm\tfrac12 以外的值,不能与平移周期性相容。

独立点对称操作:1,2,3,4,6,i,m,4ˉ1,2,3,4,6,i,m,\bar4 共 88 种。其余旋转—反演操作可由它们组合得到,例如 1ˉ=i\bar1=i,2ˉ=m\bar2=m,3ˉ=3+i\bar3=3+i,6ˉ=3+m\bar6=3+m。

点群与空间群:

  • 全部点对称操作构成点群;晶体只有 3232 种点群。
  • 点对称操作加平移操作构成空间群;晶体共有 230230 种空间群。
  • 平移对称操作还包括螺旋轴和滑移反映面。

立方体与正四面体的对称操作数:

  • 立方体:33 个立方轴 C4C_4,每个轴有 π/2,π,3π/2\pi/2,\pi,3\pi/2 三个操作,共 3×3=93\times3=9;66 个面对角线 C2C_2,每个 11 个操作,共 66;44 个体对角线 C3C_3,每个 2π/3,4π/32\pi/3,4\pi/3 两个操作,共 4×2=84\times2=8;再加不动操作。纯旋转共 2424 个;加上中心反演得到另外 2424 个,总计 4848 个。
  • 正四面体:纯旋转 1212 个;含中心反演的操作 1212 个;总计 2424 个。正四面体既无四度轴也无对称心。

易错点:

  • 晶体中不存在五度旋转轴,但准晶体可以有五重衍射对称。
  • 旋转—反演轴 nˉ\bar n 不都是独立对称素;独立的只有 88 种。
  • 立方体有 4848 个对称操作,正四面体有 2424 个,不要混淆。
  • 点群描述宏观对称性;空间群还包含平移操作。

1.8 七大晶系与十四种布拉伐格子 ​

取布拉伐原胞三个基矢 a,b,c\mathbf a,\mathbf b,\mathbf c,长度记为 a,b,ca,b,c,夹角记为 α\alpha(b\mathbf b 与 c\mathbf c)、β\beta(a\mathbf a 与 c\mathbf c)、γ\gamma(a\mathbf a 与 b\mathbf b)。

晶系基矢与夹角条件布拉伐格子
三斜a≠b≠ca\ne b\ne c,α≠β≠γ≠90∘\alpha\ne\beta\ne\gamma\ne90^\circ简单三斜(11 种)
单斜a≠b≠ca\ne b\ne c,α=γ=90∘≠β\alpha=\gamma=90^\circ\ne\beta简单单斜、底心单斜(22 种)
正交a≠b≠ca\ne b\ne c,α=β=γ=90∘\alpha=\beta=\gamma=90^\circ简单、底心、体心、面心正交(44 种)
三角a=b=ca=b=c,α=β=γ<120∘≠90∘\alpha=\beta=\gamma<120^\circ\ne90^\circ三角(11 种)
四角(正方)a=b≠ca=b\ne c,α=β=γ=90∘\alpha=\beta=\gamma=90^\circ简单四角、体心四角(22 种)
六角a=b≠ca=b\ne c,α=β=90∘\alpha=\beta=90^\circ,γ=120∘\gamma=120^\circ六角(11 种)
立方a=b=ca=b=c,α=β=γ=90∘\alpha=\beta=\gamma=90^\circ简立方、体心立方、面心立方(33 种)

合计 1+2+4+1+2+1+3=141+2+4+1+2+1+3=14 种布拉伐格子。

为什么没有底心四方、面心四方、体心单斜、面心单斜? 因为这些格子可以通过重新选取基矢化归为更简单的格子,不是新的布拉伐格子。课件在「补充内容」中给出了图示说明。

易错点:

  • 七大晶系是按点群对称性划分的;十四种布拉伐格子是在此基础上再考虑带心方式得到的。
  • 三角晶系与六角晶系不同:三角的 a=b=ca=b=c,六角的 a=b≠ca=b\ne c。
  • 立方晶系不等于「立方体外形」;关键是 a=b=ca=b=c 且三个夹角均为 90∘90^\circ。

1.9 易错点与常考题型 ​

易错点

  1. 把晶体外形当作品种特征;实际上晶面大小、形状随生长条件变化,晶面角守恒才是特征。
  2. fcc、bcc、hcp 的配位数与最近邻距离记混:fcc 为 1212、2a/2\sqrt2a/2;bcc 为 88、3a/2\sqrt3a/2;hcp 为 1212、aa。
  3. 致密度计算中漏掉晶胞原子数贡献因子(面 1/21/2、棱 1/41/4、顶角 1/81/8)。
  4. 把 CsCl 说成体心立方点阵;把 NaCl 说成简立方点阵。
  5. 晶列指数与晶面指数括号混用;密勒指数忘记取倒数。
  6. 倒格子定义漏 2π2\pi;把 bi\mathbf b_i 当作与 ai\mathbf a_i 平行。
  7. 面间距公式 d=2π/∣G∣d=2\pi/\lvert\mathbf G\rvert 与晶面指数定义不匹配(原胞指数 vs 密勒指数)。
  8. 认为晶体可以有五度旋转轴;把 nˉ\bar n 都当作独立对称素。
  9. 把「点群 3232 种」与「空间群 230230 种」混淆。

常考题型

  1. 给定晶格常数,求 fcc/bcc/hcp 的致密度、配位数、最近邻距离。
  2. 给定晶胞中的原子坐标,判断结构型(NaCl、CsCl、ZnS、金刚石)。
  3. 求给定晶列/晶面的指数;由指数画晶面。
  4. 证明正格子与倒格子的关系;求 fcc/bcc 的倒格子;证明面间距公式。
  5. 判断简单晶格/复式晶格;求原胞基矢与原胞体积。
  6. 列举立方体/正四面体的全部对称操作并计数。
  7. 判断七种晶系、十四种布拉伐格子;解释为什么没有某种带心格子。

第 2 章 晶体的结合 ​

2.1 原子量子数与电负性 ​

2.1.1 四个量子数 ​

量子数符号取值物理意义
主量子数nnn=1,2,3,⋯n=1,2,3,\cdots主要决定能量,表示电子离核的平均距离,对应 K,L,M,N,⋯K,L,M,N,\cdots 壳层
角量子数lll=0,1,2,⋯ ,n−1l=0,1,2,\cdots,n-1决定轨道角动量大小和电子云形状,对应 s,p,d,f,⋯s,p,d,f,\cdots
磁量子数mlm_lml=0,±1,±2,⋯ ,±lm_l=0,\pm1,\pm2,\cdots,\pm l决定轨道角动量在外磁场方向的分量
自旋量子数msm_sms=±12m_s=\pm\tfrac12决定自旋角动量在外磁场方向的分量

重要公式:氢原子能级

En=−13.6n2 eVE_n=-\frac{13.6}{n^2}\ \mathrm{eV}

其中 EnE_n 的单位为 eV\mathrm{eV},nn 为主量子数。适用条件:氢原子或类氢离子;多电子原子中能量还取决于 ll。

重要公式:角动量及其分量

∣L∣=ℏl(l+1),Lz=mlℏ\lvert\mathbf L\rvert=\hbar\sqrt{l(l+1)},\qquad L_z=m_l\hbar

其中 ∣L∣\lvert\mathbf L\rvert 的单位为 J⋅s\mathrm{J\cdot s},ℏ=1.0546×10−34 J⋅s\hbar=1.0546\times10^{-34}\ \mathrm{J\cdot s};ll 为角量子数,mlm_l 为磁量子数。适用条件:单电子原子;多电子原子需考虑耦合。

壳层与支壳层电子数:

Nn=∑l=0n−12(2l+1)=2n2N_n=\sum_{l=0}^{n-1}2(2l+1)=2n^2

其中 NnN_n 为第 nn 壳层最多容纳的电子数,无量纲;2(2l+1)2(2l+1) 为给定 ll 的支壳层状态数。适用条件:泡利不相容原理。

核外电子排布三原则:

  1. 能量最低原理:电子优先占据能量最低的能级。
  2. 泡利不相容原理:一个原子轨道最多容纳 22 个电子,且自旋相反。
  3. 洪特规则:总自旋 SS 越大能量越低;SS 相同时总轨道角动量 LL 越大能量越低;S,LS,L 都相同时,未满半壳层或刚好半壳层取 JJ 越小能量越低(正序),否则 JJ 越大能量越低(反序)。

待核对:课件第 12 页在写磁量子数时同时出现「mlm_l 取绝对值不大于 ll 的整数」与「−(n−1),⋯ ,n−1-(n-1),\cdots,n-1」两种说法,后者与前者矛盾。正确范围是 ml=0,±1,⋯ ,±lm_l=0,\pm1,\cdots,\pm l,共 2l+12l+1 个。

2.1.2 电离能、亲和能与电负性 ​

  • 电离能 II:原子失去一个电子所需能量。第一电离能从中性原子移去第一个电子;第二电离能从正一价离子再移去一个电子。II 越大,电子越难失去。
  • 亲和能 AA:原子捕获一个电子放出的能量。AA 越大,越易得电子。亲和过程不是电离过程的逆过程。
  • 电负性 χ\chi:表征原子束缚电子能力的量。χ\chi 大表示吸引电子能力强,χ\chi 小表示易失去电子。

重要公式:穆力肯电负性

χ=0.18 (I+A)\chi=0.18\,(I+A)

其中 II 为电离能,AA 为亲和能;χ\chi 无量纲(按 Pauling 标度)。适用条件:课件给出的经验式;I,AI,A 通常取 eV\mathrm{eV}。待核对:课件未注明 I,AI,A 的单位,若取 kJ/mol\mathrm{kJ/mol} 则系数不同。

周期趋势:同一周期从左到右电负性增大;同一族从下到上电负性增大。电负性可综合衡量元素的金属性与非金属性。

易错点:

  • 电离能越大越难失电子,不是越易失电子。
  • 亲和能是「放出」的能量,符号约定要与教材统一。
  • 电负性没有单位;不要把 Pauling 标度与 Mulliken 标度的数值混用。

2.2 结合能与结合力的一般理论 ​

结合能:自由原子(或分子)结合成晶体时放出的能量。

重要公式:结合能

W=E0−ENW=E_0-E_N

其中 WW 为结合能(J\mathrm{J} 或 eV\mathrm{eV});E0E_0 为自由原子状态的总能量;ENE_N 为晶体平衡态的总能量。适用条件:忽略温度、压强对能量的影响。

晶体总相互作用能:以第 ii 个原子为中心,其与其他原子的相互作用势能之和为 uiu_i;由于每对相互作用只能计算一次,

U=12N∑j≠iu(rij)=12NuiU=\frac12N\sum_{j\ne i}u(r_{ij})=\frac12N u_i

其中 UU 为晶体总势能(J\mathrm{J});NN 为原子总数;u(rij)u(r_{ij}) 为第 i,ji,j 个原子间的相互作用势(J\mathrm{J});rijr_{ij} 为原子间距(m\mathrm{m})。因子 12\tfrac12 是因为每对相互作用被两个原子各计算了一次。

两原子间的相互作用势与力:

重要公式:伦纳德—琼斯型相互作用势

u(r)=−Arm+Brnu(r)=-\frac{A}{r^m}+\frac{B}{r^n}

其中 uu 的单位为 J\mathrm{J};rr 为原子间距(m\mathrm{m});A,BA,B 为大于零的常数,AA 的单位为 J⋅mm\mathrm{J\cdot m^m},BB 的单位为 J⋅mn\mathrm{J\cdot m^n};m,nm,n 为大于零的幂指数。第一项为吸引项(负),第二项为排斥项(正)。适用条件:惰性气体、离子晶体等近似球对称的闭壳层相互作用。

重要公式:相互作用力

f(r)=−dudr=−mArm+1+nBrn+1f(r)=-\frac{\mathrm du}{\mathrm dr} =-\frac{mA}{r^{m+1}}+\frac{nB}{r^{n+1}}

其中 ff 的单位为 N\mathrm{N}(沿 rr 方向,正为排斥、负为吸引);rr 的单位为 m\mathrm{m}。适用条件:有心力、准静态。

平衡位置与稳定性:

dudr∣r=r0=0⟹mAr0m+1=nBr0n+1⟹r0n−m=nBmA\left.\frac{\mathrm du}{\mathrm dr}\right|_{r=r_0}=0 \quad\Longrightarrow\quad \frac{mA}{r_0^{m+1}}=\frac{nB}{r_0^{n+1}} \quad\Longrightarrow\quad r_0^{n-m}=\frac{nB}{mA}

其中 r0r_0 为平衡原子间距(m\mathrm{m})。在 r=r0r=r_0 处,u(r)u(r) 取最小值,f(r0)=0f(r_0)=0。

重要公式:平衡时的势能与 n>mn>m

u(r0)=−Ar0m(1−mn)u(r_0)=-\frac{A}{r_0^m}\left(1-\frac mn\right)

d2udr2∣r0=m(n−m)Ar0m+2>0⟹n>m\left.\frac{\mathrm d^2u}{\mathrm dr^2}\right|_{r_0} =\frac{m(n-m)A}{r_0^{m+2}}>0 \quad\Longrightarrow\quad n>m

其中 u(r0)u(r_0) 的单位为 J\mathrm{J};r0r_0 的单位为 m\mathrm{m}。n>mn>m 表明排斥作用随距离衰减得比吸引作用更快,即排斥作用是短程的。

证明:由平衡条件 B=mAnr0n−mB=\dfrac{mA}{n}r_0^{n-m},代入 u(r0)u(r_0):

u(r0)=−Ar0m+Br0n=−Ar0m+mAnr0n−mr0n=−Ar0m(1−mn)u(r_0)=-\frac{A}{r_0^m}+\frac{B}{r_0^n} =-\frac{A}{r_0^m}+\frac{mA}{n}\frac{r_0^{n-m}}{r_0^n} =-\frac{A}{r_0^m}\left(1-\frac mn\right)

二阶导数:

u′′(r)=−m(m+1)Arm+2+n(n+1)Brn+2u''(r)=-\frac{m(m+1)A}{r^{m+2}}+\frac{n(n+1)B}{r^{n+2}}

代入 BB:

u′′(r0)=Ar0m+2[−m(m+1)+m(n+1)]=m(n−m)Ar0m+2u''(r_0)=\frac{A}{r_0^{m+2}}\left[-m(m+1)+m(n+1)\right] =\frac{m(n-m)A}{r_0^{m+2}}

稳定平衡要求 u′′(r0)>0u''(r_0)>0,故 n>mn>m。

最大有效引力:在 r=rmr=r_m 处引力最大,条件为

dfdr∣rm=0⟺d2udr2∣rm=0\left.\frac{\mathrm df}{\mathrm dr}\right|_{r_m}=0 \quad\Longleftrightarrow\quad \left.\frac{\mathrm d^2u}{\mathrm dr^2}\right|_{r_m}=0

其中 rmr_m 为最大引力对应的间距(m\mathrm{m})。通常 rm<r0r_m<r_0。

体积弹性模量与压缩系数:

重要公式:体积弹性模量

K=−V(∂P∂V)T=Vd2UdV2∣V0K=-V\left(\frac{\partial P}{\partial V}\right)_T =V\left.\frac{\mathrm d^2U}{\mathrm dV^2}\right|_{V_0}

其中 KK 的单位为 Pa\mathrm{Pa}(N/m2\mathrm{N/m^2});VV 为晶体体积(m3\mathrm{m^3});V0V_0 为平衡体积;PP 为压强(Pa\mathrm{Pa});UU 为总势能(J\mathrm{J})。适用条件:绝热或等温过程近似、小形变。

重要公式:压缩系数

κ=−1V(∂V∂P)T=1K\kappa=-\frac1V\left(\frac{\partial V}{\partial P}\right)_T=\frac1K

其中 κ\kappa 的单位为 Pa−1\mathrm{Pa^{-1}};KK 的单位为 Pa\mathrm{Pa}。适用条件:各向同性压缩、小形变。

证明:由热力学第一定律 dU=TdS−PdV\mathrm dU=T\mathrm dS-P\mathrm dV,在 T=0T=0 或绝热近似下 dU=−PdV\mathrm dU=-P\mathrm dV,故 P=−dU/dVP=-\mathrm dU/\mathrm dV。于是

K=−VdPdV=Vd2UdV2K=-V\frac{\mathrm dP}{\mathrm dV} =V\frac{\mathrm d^2U}{\mathrm dV^2}

易错点:

  • 总势能 UU 前的因子 12\tfrac12 不能漏;uiu_i 是「一个原子与所有其他原子」的势能和。
  • u(r)u(r) 中吸引项为负、排斥项为正;f=−du/drf=-du/dr,负号不能丢。
  • 平衡条件是 u′(r0)=0u'(r_0)=0,不是 u(r0)=0u(r_0)=0。
  • n>mn>m 是稳定性要求,不是任意假设。
  • KK 的单位是 Pa\mathrm{Pa},κ\kappa 的单位是 Pa−1\mathrm{Pa^{-1}},二者互为倒数。

2.3 五种结合类型 ​

键型晶体类型结合能(kcal/mol\mathrm{kcal/mol})熔点硬度导电性方向性/饱和性典型晶体特性
离子键离子晶体150150–370370高高固态差,熔融或水解后导电无方向性(球对称)NaCl、CsCl易碎,红外有特征峰,可见光透明
共价键原子晶体125125–300300高高绝缘体或半导体有方向性、饱和性C(金刚石)、Si、Ge硬而脆
金属键金属晶体2525–200200有高有低有高有低好无方向性Cu、Ag、Au、Na、W延展性、金属光泽
范德瓦尔斯键分子晶体<5<5低低绝缘有极性惰性元素固体、CO$_2$结合能低、易挥发
氢键氢键晶体55–1515低低有有饱和性H$_2$O(冰)、HF结合能强于范德瓦尔斯力

2.3.1 离子晶体 ​

  • 结构:正负离子交替排列,典型为 NaCl、CsCl;II-VI 族化合物(CdS、ZnS)和 I-VII 族化合物(NaCl、CsCl)常视为离子晶体。
  • 结合力:正负离子间的静电库仑力;电子组态与惰性气体相同,具有球形对称性,价键无方向性。
  • 性质:结合能大(150150–370 kcal/mol370\ \mathrm{kcal/mol}),结构稳定;导电性差、熔点高、硬度高、热膨胀系数小;红外区有特征吸收峰,可见光透明;受力时易沿解理面碎裂(脆性)。

脆性解释:受力前近邻离子电性相反,引力为主;受力后近邻离子电性相同,斥力为主,故易碎。

2.3.2 原子晶体(共价晶体) ​

  • 典型:IV 族元素 C(金刚石)、Si、Ge、Sn(灰锡)。
  • 结合力:共价键,即共用电子对;具有饱和性和方向性。
  • 饱和性:一个原子能形成的共价键数目有最大值。价电子壳层未达半满(少于 44 个)时,成键数等于价电子数;价电子数 N≥4N\ge4 时,成键数符合 (8−N)(8-N) 规则。
  • 方向性:在电子云重叠密度大的方向成键;头对头形成 σ\sigma 键(牢固),肩并肩形成 π\pi 键(易断)。
  • 杂化:成键时能量相近的价电子轨道重新组合,形成杂化轨道。杂化前后轨道数目不变,但伸展方向和形状改变。金刚石为 sp3sp^3 杂化,44 个杂化轨道指向正四面体顶点,键角 109∘28′109^\circ28';石墨中碳为 sp2sp^2 杂化,形成平面三角形,未杂化的 pp 轨道垂直于平面,形成 π\pi 键。
  • 性质:熔点高、导电性差、硬而脆。金刚石是良好绝缘体;Si、Ge 在极低温下才是绝缘体,电阻率随温度升高迅速下降,是典型半导体。

2.3.3 金属晶体 ​

  • 典型:I、II 族元素及过渡元素。
  • 结合力:原子实与公有化电子云之间的静电库仑力;价键无方向性,故要求排列最紧密。
  • 结构:常见 fcc(Cu、Ag、Au、Al,CN =12=12)、hcp(Be、Mg、Zn、Cd,CN =12=12)、bcc(Li、Na、K、Rb、Cs、Mo、W,CN =8=8)。
  • 性质:导电性好、导热性好、有金属光泽、延展性好。延展性来源于金属键无方向性,原子层间可相对滑动而不易断裂。

2.3.4 分子晶体与范德瓦尔斯力 ​

  • 典型:惰性元素固体、CO$_2$ 等。
  • 结合力:非极性或极性分子间的瞬时/固有偶极相互作用,即范德瓦尔斯力。包括:
    1. 葛生(Keesom)力:固有偶极矩间的取向作用;
    2. 德拜(Debye)力:诱导偶极矩间的诱导力;
    3. 伦敦(London)力:瞬态偶极矩间的色散力。
  • 结构:以分子(或原子)为单元密堆积。
  • 性质:结合能很低(<5 kcal/mol<5\ \mathrm{kcal/mol}),熔点低、硬度低、导电性差、可见光透明。

2.3.5 氢键晶体 ​

  • 典型:冰(H$_2$O)、HF。
  • 结合形式:X−H⋯YX-\mathrm H\cdots Y,其中 X,YX,Y 为 F、O、N 等电负性大、原子半径小的非金属原子。
  • 特点:氢键具有饱和性;结合能一般低于化学键但强于范德瓦尔斯力(55–15 kcal/mol15\ \mathrm{kcal/mol})。
  • 反常性质:HF 的沸点比 HCl、HBr、HI 高,原因是存在氢键。

2.3.6 混合键 ​

  • 石墨:层内为共价键(sp2sp^2 杂化,形成六角网格),层间为范德瓦尔斯力,同时具有类似金属的导电性(π\pi 电子可在层内运动),是三种键共存的混合键晶体。
  • C$_{60}$:分子间为范德瓦尔斯力,分子内为共价键。

易错点:

  • 离子晶体在固态不导电,但熔融或溶于水后导电;不能说「离子晶体导电性差」而不加条件。
  • 共价晶体不全是绝缘体;Si、Ge 是半导体。
  • 金属键无方向性、无饱和性;共价键有方向性、有饱和性。
  • 范德瓦尔斯力包括取向力、诱导力、色散力三种,不只是一种。
  • 氢键不是完整的化学键,但比范德瓦尔斯力强,且有饱和性。

2.4 马德隆常数与离子晶体结合能 ​

离子间相互作用能:设第 ii 个离子与第 jj 个离子的间距为 rijr_{ij},则

u(rij)=±q24πε0rij+brijnu(r_{ij})=\pm\frac{q^2}{4\pi\varepsilon_0 r_{ij}}+\frac{b}{r_{ij}^n}

其中 qq 为离子电荷量(C\mathrm C);ε0=8.854×10−12 F/m\varepsilon_0=8.854\times10^{-12}\ \mathrm{F/m};rijr_{ij} 为离子间距(m\mathrm m);bb 为排斥项系数(J⋅mn\mathrm{J\cdot m^n});nn 为排斥指数。第一项中异号离子取正(吸引),同号离子取负(排斥);第二项为重叠排斥能。适用条件:刚性离子近似、球对称、非极化。

马德隆常数:设最近邻距离为 rr,任意离子间距可写成

rij=ajrr_{ij}=a_j r

其中 aja_j 为无量纲的几何因子。则总内能

重要公式:离子晶体总内能

U=N∑j≠i(±q24πε0rij+brijn)=−N[αq24πε0r−Brn]U=N\sum_{j\ne i}\left(\pm\frac{q^2}{4\pi\varepsilon_0 r_{ij}}+\frac{b}{r_{ij}^n}\right) =-N\left[\frac{\alpha q^2}{4\pi\varepsilon_0 r}-\frac{B}{r^n}\right]

α=∑j≠i(±1aj),B=∑j≠ibajn\alpha=\sum_{j\ne i}\left(\pm\frac{1}{a_j}\right),\qquad B=\sum_{j\ne i}\frac{b}{a_j^n}

其中 UU 的单位为 J\mathrm J;NN 为离子对数;α\alpha 为马德隆常数(无量纲);BB 的单位为 J⋅mn\mathrm{J\cdot m^n};rr 为最近邻距离(m\mathrm m)。符号约定:同号离子取负、异号离子取正(课件写作「同号为负,异号为正」)。适用条件:离子晶体、刚性离子模型。

马德隆常数的计算(埃夫琴晶胞法):把晶体分成许多大的中性离子组(埃夫琴单胞),每组正负离子数相同;以某一离子为参考离子,只计算本组内各离子对它的能量贡献,求和无量纲因子。

重要公式:NaCl 的马德隆常数

αNaCl=6−122+83−62+245−⋯≈1.7476\alpha_{\mathrm{NaCl}} =6-\frac{12}{\sqrt2}+\frac{8}{\sqrt3}-\frac{6}{2}+\frac{24}{\sqrt5}-\cdots \approx1.7476

其中各项依次为最近邻(66 个异号)、次近邻(1212 个同号)、第三近邻(88 个异号)等的贡献;α\alpha 无量纲。适用条件:NaCl 型结构、无限大晶体、级数收敛求和。

常用马德隆常数:

αNaCl≈1.7476,αCsCl≈1.7627,αZnS≈1.6381\alpha_{\mathrm{NaCl}}\approx1.7476,\qquad \alpha_{\mathrm{CsCl}}\approx1.7627,\qquad \alpha_{\mathrm{ZnS}}\approx1.6381

其中 α\alpha 均无量纲;分别对应 NaCl、CsCl、ZnS 结构。

补充:离子晶体的平衡间距与结合能(由课件通式推出)

由 dU/dr=0\mathrm dU/\mathrm dr=0:

r0n−1=nBαq2/(4πε0)r_0^{n-1}=\frac{nB}{\alpha q^2/(4\pi\varepsilon_0)}

U(r0)=−Nαq24πε0r0(1−1n)U(r_0)=-\frac{N\alpha q^2}{4\pi\varepsilon_0 r_0}\left(1-\frac1n\right)

其中 r0r_0 为平衡最近邻距离(m\mathrm m);U(r0)U(r_0) 为平衡结合能(J\mathrm J);其余符号同上。适用条件:离子晶体、n>1n>1、忽略零点能。

易错点:

  • 马德隆常数 α\alpha 是无量纲的几何因子,与 q,rq,r 无关;不同结构数值不同。
  • 总内能中的 NN 是离子对数;若用离子总数,则需注意因子。
  • 求和时必须按「异号为正、同号为负」处理,且各近邻的数目与距离因子不能漏。
  • 埃夫琴法要求选取中性离子组,否则会引入额外电荷,收敛很慢。
  • 马德隆级数条件收敛,直接按球壳求和结果依赖求和顺序;实际用中性组或 Evjen 方法。

2.5 易错点与常考题型 ​

易错点

  1. 把结合能 W=E0−ENW=E_0-E_N 的符号写反;结合能是放出的能量,应为正值。
  2. 总势能 U=12NuiU=\tfrac12N u_i 漏掉 12\tfrac12。
  3. 平衡条件写成 u(r0)=0u(r_0)=0;应为 u′(r0)=0u'(r_0)=0。
  4. 由 u′′(r0)>0u''(r_0)>0 推出 n>mn>m 时,符号或代数出错。
  5. 把离子晶体的「固态不导电」误写成「不导电」。
  6. 把 CsCl 当作 bcc 点阵;把 NaCl 当作简立方。
  7. 马德隆常数计算中近邻数、距离因子、正负号三者之一出错。
  8. 电负性、电离能、亲和能的概念混淆;亲和过程不是电离过程的逆过程。

常考题型

  1. 由 u(r)=−A/rm+B/rnu(r)=-A/r^m+B/r^n 求 r0r_0、u(r0)u(r_0)、最大引力位置 rmr_m,证明 n>mn>m。
  2. 由 U(V)U(V) 求体积弹性模量 KK 与压缩系数 κ\kappa。
  3. 计算 NaCl 型晶体的马德隆常数(给定近邻表)。
  4. 判断晶体类型并解释熔点、硬度、导电性、脆性、延展性。
  5. 用杂化轨道理论解释金刚石 sp3sp^3、石墨 sp2sp^2 的结构与性质。
  6. 比较五种结合类型的结合能、方向性、饱和性和典型晶体。

第 3 章 晶格振动与晶格热容 ​

3.1 一维单原子链 ​

3.1.1 模型与受力 ​

设一维单原子链有 NN 个原子,原子质量 mm,晶格常数 aa,第 nn 个原子偏离平衡位置的位移为 unu_n。只考虑最近邻相互作用。

弹性系数的来源:设两原子在平衡位置的相互作用势为 U(a)U(a),相对位移 δ=un+1−un\delta=u_{n+1}-u_n。在平衡位置附近作泰勒展开:

U(a+δ)=U(a)+(dUdδ)aδ+12(d2Udδ2)aδ2+⋯U(a+\delta)=U(a)+\left(\frac{\mathrm dU}{\mathrm d\delta}\right)_a\delta +\frac12\left(\frac{\mathrm d^2U}{\mathrm d\delta^2}\right)_a\delta^2+\cdots

平衡位置势能取极小,故

(dUdδ)a=0\left(\frac{\mathrm dU}{\mathrm d\delta}\right)_a=0

于是恢复力

重要公式:简谐恢复力与弹性系数

f=−(dUdδ)=−(d2Udδ2)aδ=−βδf=-\left(\frac{\mathrm dU}{\mathrm d\delta}\right) =-\left(\frac{\mathrm d^2U}{\mathrm d\delta^2}\right)_a\delta =-\beta\delta

β=(d2Udδ2)a\beta=\left(\frac{\mathrm d^2U}{\mathrm d\delta^2}\right)_a

其中 ff 为恢复力(N\mathrm N);δ\delta 为相对位移(m\mathrm m);β\beta 为弹性系数(N/m\mathrm{N/m},也写作 J/m2\mathrm{J/m^2})。适用条件:平衡位置附近的小振动、简谐近似。

第 nn 个原子的受力:右邻拉伸贡献 +β(un+1−un)+\beta(u_{n+1}-u_n),左邻拉伸贡献 −β(un−un−1)-\beta(u_n-u_{n-1}),故

重要公式:一维单原子链运动方程

md2undt2=β(un+1+un−1−2un)m\frac{\mathrm d^2u_n}{\mathrm dt^2} =\beta(u_{n+1}+u_{n-1}-2u_n)

其中 mm 为原子质量(kg\mathrm{kg});unu_n 为第 nn 个原子的位移(m\mathrm m);β\beta 为弹性系数(N/m\mathrm{N/m});tt 为时间(s\mathrm s)。适用条件:最近邻近似、简谐近似、一维单原子链。

3.1.2 试解与色散关系 ​

设试解为行波(格波):

un=Aei(ωt−qna)u_n=Ae^{i(\omega t-qna)}

其中 AA 为振幅(m\mathrm m);ω\omega 为角频率(rad/s\mathrm{rad/s});qq 为波矢(rad/m\mathrm{rad/m},一维情形常写为 m−1\mathrm{m^{-1}});nana 为第 nn 个原子的平衡位置(m\mathrm m)。

代入运动方程:

−mω2Aei(ωt−qna)=βAei(ωt−qna)(e−iqa+eiqa−2)-m\omega^2Ae^{i(\omega t-qna)} =\beta A e^{i(\omega t-qna)}\left(e^{-iqa}+e^{iqa}-2\right)

约去公因子,得

−mω2=2β(cos⁡qa−1)=−4βsin⁡2qa2-m\omega^2=2\beta(\cos qa-1)=-4\beta\sin^2\frac{qa}{2}

重要公式:一维单原子链色散关系

ω2=2βm(1−cos⁡qa)=4βmsin⁡2qa2\omega^2=\frac{2\beta}{m}(1-\cos qa) =\frac{4\beta}{m}\sin^2\frac{qa}{2}

ω(q)=2βm∣sin⁡qa2∣\omega(q)=2\sqrt{\frac{\beta}{m}}\left\lvert\sin\frac{qa}{2}\right\rvert

其中 ω\omega 为角频率(rad/s\mathrm{rad/s});β\beta 为弹性系数(N/m\mathrm{N/m});mm 为原子质量(kg\mathrm{kg});qq 为波矢(m−1\mathrm{m^{-1}});aa 为晶格常数(m\mathrm m)。适用条件:一维单原子链、最近邻、简谐近似。

重要公式:最大频率

ωmax⁡=2βm\omega_{\max}=2\sqrt{\frac{\beta}{m}}

其中 ωmax⁡\omega_{\max} 的单位为 rad/s\mathrm{rad/s}。适用条件:一维单原子链;在 qa=±πqa=\pm\pi 处取到。

色散关系的性质:

  1. 周期性:ω(q+2π/a)=ω(q)\omega(q+2\pi/a)=\omega(q)。因此 qq 只需取在一个周期内。
  2. 对称性:ω(−q)=ω(q)\omega(-q)=\omega(q)。
  3. 通常把 qq 限制在第一布里渊区:

    −πa<q≤πa-\frac{\pi}{a}<q\le\frac{\pi}{a}

    其中 qq 的单位为 m−1\mathrm{m^{-1}},aa 的单位为 m\mathrm m。

长波极限(q→0q\to0,即 qa≪1qa\ll1):

sin⁡qa2≈qa2⟹ω≈aβm q\sin\frac{qa}{2}\approx\frac{qa}{2} \quad\Longrightarrow\quad \omega\approx a\sqrt{\frac{\beta}{m}}\,q

重要公式:长波极限与声速

ω=vsq,vs=aβm\omega=v_s q,\qquad v_s=a\sqrt{\frac{\beta}{m}}

其中 vsv_s 为声速(m/s\mathrm{m/s});aa 为晶格常数(m\mathrm m);β\beta 为弹性系数(N/m\mathrm{N/m});mm 为原子质量(kg\mathrm{kg})。适用条件:qa≪1qa\ll1,即长波极限。此时格波可看作连续介质中的弹性波。

短波极限(q=±π/aq=\pm\pi/a,第一布里渊区边界):

ω=ωmax⁡=2βm,qa=π\omega=\omega_{\max}=2\sqrt{\frac{\beta}{m}},\qquad qa=\pi

相邻原子相位差为 π\pi,振动方向相反;波长 λ=2a\lambda=2a。

物理图像:长波极限下相邻原子相位差趋于零,同向振动;短波极限下相邻原子反相振动。格波不是某一个原子的振动,而是所有原子以相同频率、固定相位关系参与的集体振动。

3.1.3 玻恩—卡门周期性边界条件 ​

将 NN 个原子首尾相接成环,要求

un+N=unu_{n+N}=u_n

代入试解:

Aei[ωt−q(n+N)a]=Aei(ωt−qna)⟹e−iqaN=1Ae^{i[\omega t-q(n+N)a]}=Ae^{i(\omega t-qna)} \quad\Longrightarrow\quad e^{-iqaN}=1

重要公式:波矢的量子化

q=2πhNa,h∈Zq=\frac{2\pi h}{Na},\qquad h\in\mathbb Z

其中 qq 为波矢(m−1\mathrm{m^{-1}});NN 为原子数;aa 为晶格常数(m\mathrm m);hh 为整数。适用条件:一维 NN 原子环链、玻恩—卡门周期性边界条件。

由 −π/a<q≤π/a-\pi/a<q\le\pi/a 得

−N2<h≤N2-\frac N2<h\le\frac N2

故 qq 共有 NN 个取值,相邻 qq 的间隔为

Δq=2πNa=2πL\Delta q=\frac{2\pi}{Na}=\frac{2\pi}{L}

其中 L=NaL=Na 为链长(m\mathrm m)。每个 qq 对应一个格波模式,故一维单原子链的格波数、振动模式数均等于 NN,等于其自由度数。

易错点:

  • 色散关系是 ω2\omega^2 的表达式,开方后取绝对值;ω\omega 不能为负。
  • qq 的范围是 −π/a<q≤π/a-\pi/a<q\le\pi/a(或 −π/a≤q<π/a-\pi/a\le q<\pi/a),长度为一个倒格矢 2π/a2\pi/a,不是 4π/a4\pi/a。
  • 长波极限的声速是 aβ/ma\sqrt{\beta/m},不是 β/m\sqrt{\beta/m}。
  • 短波极限相邻原子反相,波长 2a2a,不是 aa。

3.2 一维双原子链 ​

3.2.1 模型与运动方程 ​

一维复式格子,原胞含两个原子:质量 MM 的原子编号为 2n−2,2n,2n+2,⋯2n-2,2n,2n+2,\cdots,质量 mm 的原子编号为 2n−1,2n+1,⋯2n-1,2n+1,\cdots。原胞长度(晶格常数)为 2a2a,最近邻弹性系数为 β\beta。

运动方程:

{Md2u2ndt2=β(u2n+1+u2n−1−2u2n)md2u2n+1dt2=β(u2n+2+u2n−2u2n+1)\begin{cases} M\dfrac{\mathrm d^2u_{2n}}{\mathrm dt^2} =\beta(u_{2n+1}+u_{2n-1}-2u_{2n})\\[8pt] m\dfrac{\mathrm d^2u_{2n+1}}{\mathrm dt^2} =\beta(u_{2n+2}+u_{2n}-2u_{2n+1}) \end{cases}

其中 M,mM,m 为原子质量(kg\mathrm{kg});uu 为位移(m\mathrm m);β\beta 为弹性系数(N/m\mathrm{N/m})。适用条件:最近邻、简谐近似、一维双原子链。

3.2.2 试解与色散关系 ​

设

u2n=Aei(ωt−2nqa),u2n+1=Bei[ωt−(2n+1)qa]u_{2n}=Ae^{i(\omega t-2nqa)},\qquad u_{2n+1}=Be^{i[\omega t-(2n+1)qa]}

其中 A,BA,B 分别为轻重原子的振幅(m\mathrm m);qq 为波矢(m−1\mathrm{m^{-1}});ω\omega 为角频率(rad/s\mathrm{rad/s})。

代入运动方程并约去公因子,得

{(2β−Mω2)A−2βcos⁡(qa)B=0−2βcos⁡(qa)A+(2β−mω2)B=0\begin{cases} (2\beta-M\omega^2)A-2\beta\cos(qa)B=0\\[4pt] -2\beta\cos(qa)A+(2\beta-m\omega^2)B=0 \end{cases}

A,BA,B 有非零解的条件是系数行列式为零:

∣2β−Mω2−2βcos⁡(qa)−2βcos⁡(qa)2β−mω2∣=0\begin{vmatrix} 2\beta-M\omega^2 & -2\beta\cos(qa)\\ -2\beta\cos(qa) & 2\beta-m\omega^2 \end{vmatrix}=0

展开:

Mmω4−2β(M+m)ω2+4β2sin⁡2(qa)=0Mm\omega^4-2\beta(M+m)\omega^2+4\beta^2\sin^2(qa)=0

重要公式:一维双原子链色散关系

ω±2=βM+mMm{1±[1−4Mm(M+m)2sin⁡2(qa)]1/2}\omega_{\pm}^2 =\beta\frac{M+m}{Mm} \left\{1\pm\left[1-\frac{4Mm}{(M+m)^2}\sin^2(qa)\right]^{1/2}\right\}

其中 ω+\omega_+ 对应光学支,ω−\omega_- 对应声学支;ω\omega 的单位为 rad/s\mathrm{rad/s};M,mM,m 为原子质量(kg\mathrm{kg});β\beta 为弹性系数(N/m\mathrm{N/m});qq 为波矢(m−1\mathrm{m^{-1}});aa 为半原胞长度(m\mathrm m)。适用条件:一维双原子链、最近邻、简谐近似。

折合质量:

1μ=1M+1m⟹μ=MmM+m\frac1\mu=\frac1M+\frac1m \quad\Longrightarrow\quad \mu=\frac{Mm}{M+m}

其中 μ\mu 为折合质量(kg\mathrm{kg})。

支的命名:

  • ω−\omega_-:q→0q\to0 时 ω−→0\omega_-\to0,称为声学支(声频支)。
  • ω+\omega_+:q=0q=0 时 ω+≠0\omega_+\ne0,称为光学支(光频支)。

边界值:

  • q=0q=0:

    ω−=0,ω+=2βμ=2β(1M+1m)\omega_-=0,\qquad \omega_+=\sqrt{\frac{2\beta}{\mu}} =\sqrt{2\beta\left(\frac1M+\frac1m\right)}

  • q=±π2aq=\pm\dfrac{\pi}{2a}(第一布里渊区边界):

    ω+=2βmin⁡(M,m),ω−=2βmax⁡(M,m)\omega_+=\sqrt{\frac{2\beta}{\min(M,m)}},\qquad \omega_-=\sqrt{\frac{2\beta}{\max(M,m)}}

    两支之间出现频率间隙。

波矢范围:由玻恩—卡门条件 un+N=unu_{n+N}=u_n(NN 为原胞数),

q=2πsNa,−N2<s≤N2q=\frac{2\pi s}{Na},\qquad -\frac{N}{2}<s\le\frac{N}{2}

故

−π2a<q≤π2a-\frac{\pi}{2a}<q\le\frac{\pi}{2a}

其中 qq 的单位为 m−1\mathrm{m^{-1}}。qq 有 NN 个取值,每个 qq 对应两个频率(两支),故总模式数为 2N2N,等于 2N2N 个原子的自由度数。

3.2.3 长波极限与振幅关系 ​

长声学波(q→0q\to0):

ω−≈a2βM+m q\omega_-\approx a\sqrt{\frac{2\beta}{M+m}}\,q

振幅比:

(BA)ω−=2β−Mω−22βcos⁡(qa)≈1\left(\frac BA\right)_{\omega_-} =\frac{2\beta-M\omega_-^2}{2\beta\cos(qa)} \approx1

即长声学波中,原胞内轻重原子运动完全一致,振幅和相位无差别;长声学波代表原胞质心的运动。

长光学波(q→0q\to0):

ω+≈2βμ=2β(1M+1m)\omega_+\approx\sqrt{\frac{2\beta}{\mu}} =\sqrt{2\beta\left(\frac1M+\frac1m\right)}

振幅比:

(BA)ω+=2β−Mω+22βcos⁡(qa)≈−Mm\left(\frac BA\right)_{\omega_+} =\frac{2\beta-M\omega_+^2}{2\beta\cos(qa)} \approx-\frac{M}{m}

即

MA+mB=0MA+mB=0

长光学波中,原胞内轻重原子振动方向相反,质心近似不动。

易错点:

  • 双原子链的晶格常数是 2a2a,故布里渊区边界是 ±π/(2a)\pm\pi/(2a),不是 ±π/a\pm\pi/a。
  • 声学支在 q=0q=0 时 ω=0\omega=0;光学支在 q=0q=0 时 ω≠0\omega\ne0。
  • 光学支频率不一定总高于声学支的任意点,但在同一 qq 下 ω+>ω−\omega_+>\omega_-;两支之间有禁带。
  • 振幅比公式中的 A,BA,B 是复振幅,负号表示反相。

3.3 三维晶格振动、声子与格波支数 ​

3.3.1 三维晶格振动 ​

设晶体有 N=N1N2N3N=N_1N_2N_3 个原胞,每个原胞含 nn 个原子。第 ll 个原胞中第 ss 个原子沿 α\alpha 方向的位移为

ul,s,α=As,αei(ωt−q⋅Rl)u_{l,s,\alpha}=A_{s,\alpha}e^{i(\omega t-\mathbf q\cdot\mathbf R_l)}

其中 As,αA_{s,\alpha} 为振幅(m\mathrm m);q\mathbf q 为波矢(m−1\mathrm{m^{-1}});Rl\mathbf R_l 为第 ll 个原胞的位矢(m\mathrm m);ω\omega 为角频率(rad/s\mathrm{rad/s})。

由玻恩—卡门边界条件:

q=h1N1b1+h2N2b2+h3N3b3\mathbf q=\frac{h_1}{N_1}\mathbf b_1+\frac{h_2}{N_2}\mathbf b_2+\frac{h_3}{N_3}\mathbf b_3

其中 h1,h2,h3h_1,h_2,h_3 为整数;b1,b2,b3\mathbf b_1,\mathbf b_2,\mathbf b_3 为倒格子基矢(m−1\mathrm{m^{-1}});NiN_i 为三个方向的元胞数。q\mathbf q 共有 NN 个取值,在简约布里渊区内准连续。

波矢密度:每个波矢点占据的体积为

ΔVq=(2π)3NΩ=(2π)3Vc\Delta V_q=\frac{(2\pi)^3}{N\Omega}=\frac{(2\pi)^3}{V_c}

其中 ΔVq\Delta V_q 的单位为 m−3\mathrm{m^{-3}};Ω\Omega 为原胞体积(m3\mathrm{m^3});Vc=NΩV_c=N\Omega 为晶体体积(m3\mathrm{m^3})。故波矢密度

ρq=Vc(2π)3\rho_q=\frac{V_c}{(2\pi)^3}

单位为 m−3\mathrm{m^{-3}}。

格波支数与模式数:

重要公式:三维晶格的格波支数

总支数:3n3n;声学支:33;光学支:3(n−1)3(n-1)。

总模式数:3nN3nN。

其中 nn 为原胞内原子数(无量纲);NN 为原胞数;3nN3nN 为晶体总自由度数。适用条件:三维晶体、每个原胞 nn 个原子、简谐近似。

例:金刚石:每个原胞含 22 个原子,n=2n=2,故有 66 支格波:33 支声学波(11 纵 22 横)、33 支光学波(11 纵 22 横);振动模式数为 6N6N。金刚石结构存在简并现象。

3.3.2 声子 ​

定义:采用谐振子模型描述晶格振动时,把格波的能量量子称为声子。

重要公式:声子的能量与准动量

E=ℏωq,p=ℏqE=\hbar\omega_{\mathbf q},\qquad \mathbf p=\hbar\mathbf q

其中 EE 为声子能量(J\mathrm J);ωq\omega_{\mathbf q} 为格波角频率(rad/s\mathrm{rad/s});p\mathbf p 为准动量(kg⋅m/s\mathrm{kg\cdot m/s});q\mathbf q 为波矢(m−1\mathrm{m^{-1}});ℏ=1.0546×10−34 J⋅s\hbar=1.0546\times10^{-34}\ \mathrm{J\cdot s}。适用条件:简谐近似、单声子过程。

谐振子能量:

En=(n+12)ℏωq,n=0,1,2,⋯E_n=\left(n+\frac12\right)\hbar\omega_{\mathbf q},\qquad n=0,1,2,\cdots

其中 EnE_n 为第 nn 个激发态能量(J\mathrm J);nn 为声子数(无量纲);零点能 12ℏωq\tfrac12\hbar\omega_{\mathbf q}。激发过程中产生 nn 个频率为 ωq\omega_{\mathbf q} 的声子,能量改变 ΔE=nℏωq\Delta E=n\hbar\omega_{\mathbf q}。

平均声子数(玻色—爱因斯坦分布):

重要公式:平均声子数

n(ω,T)=1eℏω/kBT−1n(\omega,T)=\frac{1}{e^{\hbar\omega/k_{\mathrm B}T}-1}

其中 nn 为平均声子数(无量纲);ω\omega 为角频率(rad/s\mathrm{rad/s});TT 为绝对温度(K\mathrm K);kB=1.38×10−23 J/Kk_{\mathrm B}=1.38\times10^{-23}\ \mathrm{J/K}。适用条件:热平衡、玻色子、化学势为零。

声子的性质:

  1. 声子是准粒子,不能脱离晶格单独存在。
  2. 声子能量量子为 ℏωq\hbar\omega_{\mathbf q}。
  3. 声子具有准动量 ℏq\hbar\mathbf q,但不是真实动量;所有声子的动量和为零,与外界作用时表现为单个声子的能量、动量交换。
  4. 声子可产生与湮灭。
  5. 声子服从玻色—爱因斯坦统计。
  6. 声子体现晶体中原子集体振动的能量量子化,是固体中的一种元激发。

非简谐效应:简谐近似下声子间无相互作用,是理想玻色气体;非简谐项(泰勒展开的三阶及以上项)引入声子间碰撞,使声子气体达到热平衡。非简谐效应还导致热膨胀、有限热导率等。

易错点:

  • 声子不是真实粒子,不能在真空中传播。
  • 声子准动量守恒只在相差一个倒格矢的意义上成立,即 ∑qi=Kh\sum\mathbf q_i=\mathbf K_h 也可。
  • 声子数不守恒;温度升高,平均声子数增加。
  • 光学支的声子不一定「能量更高」;在 q=0q=0 处光学支频率高于声学支,但在整个布里渊区要具体比较。

3.4 晶格振动谱的实验测定 ​

晶格振动谱(色散关系):ω(q)\omega(\mathbf q) 的关系。主要实验方法:中子的非弹性散射、X 射线散射、光的散射。

3.4.1 中子的非弹性散射 ​

入射中子质量 MnM_n,动量 P\mathbf P,能量 E=P2/2MnE=P^2/2M_n;散射后动量 P′\mathbf P',能量 E′=P′2/2MnE'=P'^2/2M_n。散射过程满足能量守恒和准动量守恒:

重要公式:中子散射的能量与准动量守恒

P22Mn−P′22Mn=±ℏω(q)\frac{P^2}{2M_n}-\frac{P'^2}{2M_n}=\pm\hbar\omega(\mathbf q)

P−P′=±ℏq+ℏKh\mathbf P-\mathbf P'=\pm\hbar\mathbf q+\hbar\mathbf K_h

其中 P,P′P,P' 为中子动量(kg⋅m/s\mathrm{kg\cdot m/s});MnM_n 为中子质量(kg\mathrm{kg});ω(q)\omega(\mathbf q) 为声子角频率(rad/s\mathrm{rad/s});q\mathbf q 为声子波矢(m−1\mathrm{m^{-1}});Kh\mathbf K_h 为倒格矢(m−1\mathrm{m^{-1}});「++」表示吸收一个声子,「−-」表示发射一个声子。适用条件:单声子非弹性散射、相干散射。

实验中固定入射中子动量 P\mathbf P,测量不同散射方向上的 P′\mathbf P',即可确定 ω(q)\omega(\mathbf q)。

3.4.2 光的散射 ​

入射光子波矢 k\mathbf k、能量 ℏΩ\hbar\Omega;散射光子波矢 k′\mathbf k'、能量 ℏΩ′\hbar\Omega'。能量与准动量守恒:

重要公式:光子散射的能量与准动量守恒

Ω−Ω′=±ω(q)\Omega-\Omega'=\pm\omega(\mathbf q)

k−k′=±q+Kh\mathbf k-\mathbf k'=\pm\mathbf q+\mathbf K_h

其中 Ω,Ω′\Omega,\Omega' 为入射、散射光子角频率(rad/s\mathrm{rad/s});k,k′\mathbf k,\mathbf k' 为光子波矢(m−1\mathrm{m^{-1}});ω(q)\omega(\mathbf q) 为声子角频率;q\mathbf q 为声子波矢。适用条件:单声子非弹性散射、透明介质。

可见光波矢约为 105 cm−110^5\ \mathrm{cm^{-1}} 量级,故相互作用的声子波矢也在 105 cm−110^5\ \mathrm{cm^{-1}} 量级,只是布里渊区中心附近很小一部分区域的声子,即长波声子。

光散射分类:

  • 光波与长声学波相互作用:布里渊散射;
  • 光波与长光学波相互作用:喇曼散射。
  • 喇曼散射又分:斯托克斯散射(ωs′<ωi\omega'_{\mathrm s}<\omega_{\mathrm i},发射声子)和反斯托克斯散射(ωs′>ωi\omega'_{\mathrm s}>\omega_{\mathrm i},吸收声子)。其中 ωi\omega_{\mathrm i} 为入射光角频率,ωs′\omega'_{\mathrm s} 为散射光角频率(rad/s\mathrm{rad/s})。

3.4.3 X 射线散射 ​

X 光光子能量约 104 eV10^4\ \mathrm{eV},声子能量约 102 eV10^2\ \mathrm{eV},能量变化很小,不易测量。但 X 射线波矢与倒格矢同数量级,测定范围可遍及整个布里渊区,可测大波矢范围的振动谱。

易错点:

  • 中子散射中「++」对应吸收声子,「−-」对应发射声子;不要与能量增减混淆。
  • 光子散射的准动量守恒可相差一个倒格矢 Kh\mathbf K_h(Umklapp 过程)。
  • 可见光只能探测布里渊区中心附近的长波声子;要测整个布里渊区需用中子或 X 射线。

3.5 晶格热容 ​

3.5.1 定义与实验规律 ​

晶体的比热:

CV=(∂E∂T)VC_V=\left(\frac{\partial E}{\partial T}\right)_V

其中 CVC_V 为定容热容(J/K\mathrm{J/K});EE 为晶体平均内能(J\mathrm J);TT 为温度(K\mathrm K);VV 为体积(m3\mathrm{m^3})。

一般情况下

CV=CVa+CVeC_V=C_V^{\mathrm a}+C_V^{\mathrm e}

其中 CVaC_V^{\mathrm a} 为晶格振动(声子)贡献,CVeC_V^{\mathrm e} 为电子贡献;通常 CVe≪CVaC_V^{\mathrm e}\ll C_V^{\mathrm a},本节只讨论晶格振动。

实验规律:

  1. 高温时,CV=3nNkBC_V=3nNk_{\mathrm B}。其中 nn 为原胞中原子数,NN 为原胞数,nNnN 为晶体中原子总数,kB=1.38×10−23 J/Kk_{\mathrm B}=1.38\times10^{-23}\ \mathrm{J/K}。
  2. 低温时,CVC_V 按 T3T^3 趋于零。

3.5.2 经典理论:杜隆—珀替定律 ​

由能量均分定理,每个自由度的平均能量为 kBTk_{\mathrm B}T。三维晶体有 3nN3nN 个自由度,故

E=3nNkBTE=3nNk_{\mathrm B}T

重要公式:杜隆—珀替定律

CV=(∂E∂T)V=3nNkBC_V=\left(\frac{\partial E}{\partial T}\right)_V=3nNk_{\mathrm B}

其中 CVC_V 的单位为 J/K\mathrm{J/K};nn 为原胞中原子数;NN 为原胞数;kBk_{\mathrm B} 为玻尔兹曼常量。适用条件:高温(T≫θDT\gg\theta_{\mathrm D})、经典极限。低温时经典理论失效。

3.5.3 量子理论 ​

晶体在简谐近似下可看成 3nN3nN 个相互独立的简谐振子。第 ii 个谐振子的能量为

Ei=(ni+12)ℏωiE_i=\left(n_i+\frac12\right)\hbar\omega_i

忽略零点能:

Ei=niℏωiE_i=n_i\hbar\omega_i

热平衡下平均声子数

nˉi=1eℏωi/kBT−1\bar n_i=\frac{1}{e^{\hbar\omega_i/k_{\mathrm B}T}-1}

故

重要公式:晶格振动总能量

E=∑i=13nNℏωieℏωi/kBT−1E=\sum_{i=1}^{3nN}\frac{\hbar\omega_i}{e^{\hbar\omega_i/k_{\mathrm B}T}-1}

其中 EE 为总振动能(J\mathrm J);ωi\omega_i 为第 ii 个模式的角频率(rad/s\mathrm{rad/s});nNnN 为原子总数;3nN3nN 为模式数。适用条件:简谐近似、忽略零点能、热平衡。

3.5.4 爱因斯坦模型 ​

模型假设:

  1. 晶体中原子的振动相互独立;
  2. 所有原子都具有同一频率 ω\omega。

晶体有 3nN3nN 个频率为 ω\omega 的振动。

重要公式:爱因斯坦模型的总能量与热容

E=3nN(ℏωeℏω/kBT−1+12ℏω)E=3nN\left(\frac{\hbar\omega}{e^{\hbar\omega/k_{\mathrm B}T}-1}+\frac12\hbar\omega\right)

CV=3nNkBeℏω/kBT(eℏω/kBT−1)2(ℏωkBT)2C_V=3nNk_{\mathrm B} \frac{e^{\hbar\omega/k_{\mathrm B}T}} {\left(e^{\hbar\omega/k_{\mathrm B}T}-1\right)^2} \left(\frac{\hbar\omega}{k_{\mathrm B}T}\right)^2

其中 EE 的单位为 J\mathrm J;CVC_V 的单位为 J/K\mathrm{J/K};nNnN 为原子总数;ω\omega 为爱因斯坦频率(rad/s\mathrm{rad/s})。适用条件:所有振子同一频率、简谐近似、忽略零点能。

爱因斯坦温度:

kBθE=ℏω⟹θE=ℏωkBk_{\mathrm B}\theta_{\mathrm E}=\hbar\omega \quad\Longrightarrow\quad \theta_{\mathrm E}=\frac{\hbar\omega}{k_{\mathrm B}}

其中 θE\theta_{\mathrm E} 为爱因斯坦温度(K\mathrm K)。大多数固体的 θE\theta_{\mathrm E} 在 100100–300 K300\ \mathrm K 范围。

爱因斯坦比热函数:

f(θET)=(θET)2eθE/T(eθE/T−1)2f\left(\frac{\theta_{\mathrm E}}{T}\right) =\left(\frac{\theta_{\mathrm E}}{T}\right)^2 \frac{e^{\theta_{\mathrm E}/T}}{\left(e^{\theta_{\mathrm E}/T}-1\right)^2}

CV=3nNkBf(θET)C_V=3nNk_{\mathrm B}f\left(\frac{\theta_{\mathrm E}}{T}\right)

高低温极限(令 x=θE/Tx=\theta_{\mathrm E}/T,f=x2ex/(ex−1)2f=x^2e^x/(e^x-1)^2):

  • 高温 T≫θET\gg\theta_{\mathrm E},即 x≪1x\ll1:ex≈1+x+x2/2e^x\approx1+x+x^2/2,故 f≈1f\approx1,

    CV≈3nNkBC_V\approx3nNk_{\mathrm B}

    与杜隆—珀替定律一致。
  • 低温 T≪θET\ll\theta_{\mathrm E},即 x≫1x\gg1:ex≫1e^x\gg1,故 f≈x2e−xf\approx x^2e^{-x},

    CV≈3nNkB(θET)2e−θE/TC_V\approx3nNk_{\mathrm B} \left(\frac{\theta_{\mathrm E}}{T}\right)^2 e^{-\theta_{\mathrm E}/T}

    CVC_V 比 T3T^3 更快地趋于零。

低温不符合实验的原因:爱因斯坦频率 ωE∼1013 Hz\omega_{\mathrm E}\sim10^{13}\ \mathrm{Hz},处于远红外光频区,相当于长光学波极限。甚低温度下,格波频率很低,属于长声学波;晶体的比热主要由长声学波决定。爱因斯坦模型把所有模式取为同一高频,忽略了声学支的低频模式,故低温下与实验不符。

3.5.5 模式密度(频率分布函数) ​

定义:单位频率间隔内的振动模式数。

ρ(ω)=lim⁡Δω→0ΔnΔω\rho(\omega)=\lim_{\Delta\omega\to0}\frac{\Delta n}{\Delta\omega}

其中 ρ(ω)\rho(\omega) 的单位为 s/rad\mathrm{s/rad}(或 s\mathrm{s},视 ω\omega 用角频率还是频率);Δn\Delta n 为 [ω,ω+Δω][\omega,\omega+\Delta\omega] 内的模式数。

归一化条件:

∫0ωmρ(ω) dω=3nN\int_0^{\omega_m}\rho(\omega)\,\mathrm d\omega=3nN

其中 ωm\omega_m 为最高频率(截止频率,rad/s\mathrm{rad/s});nNnN 为原子总数。

重要公式:模式密度的一般表达式

ρ(ω)=∑α=13nVc(2π)3∫Sωds∣∇qωα(q)∣\rho(\omega)=\sum_{\alpha=1}^{3n} \frac{V_c}{(2\pi)^3} \int_{S_\omega}\frac{\mathrm ds}{\left\lvert\nabla_{\mathbf q}\omega_\alpha(\mathbf q)\right\rvert}

其中 ρ(ω)\rho(\omega) 的单位为 s/rad\mathrm{s/rad};VcV_c 为晶体体积(m3\mathrm{m^3});α\alpha 为格波支指标;SωS_\omega 为等频面;ds\mathrm ds 为等频面上的面积元(m−2\mathrm{m^{-2}});∇qωα\nabla_{\mathbf q}\omega_\alpha 为频率梯度(m/s\mathrm{m/s})。适用条件:三维晶体、等频面非奇异(∇qω≠0\nabla_{\mathbf q}\omega\ne0)。

推导思路:两等频面 ω\omega 与 ω+dω\omega+\mathrm d\omega 之间的波矢数为

Δn=Vc(2π)3∫Sωds dq\Delta n=\frac{V_c}{(2\pi)^3}\int_{S_\omega}\mathrm ds\,\mathrm dq

由 dω=∣∇qω∣ dq\mathrm d\omega=\lvert\nabla_{\mathbf q}\omega\rvert\,\mathrm dq 得 dq=dω/∣∇qω∣\mathrm dq=\mathrm d\omega/\lvert\nabla_{\mathbf q}\omega\rvert,故

ρ(ω)=Δndω=Vc(2π)3∫Sωds∣∇qω∣\rho(\omega)=\frac{\Delta n}{\mathrm d\omega} =\frac{V_c}{(2\pi)^3}\int_{S_\omega} \frac{\mathrm ds}{\lvert\nabla_{\mathbf q}\omega\rvert}

一维单原子链的模式密度:

ω=ωm∣sin⁡qa2∣,ωm=2βm\omega=\omega_m\left\lvert\sin\frac{qa}{2}\right\rvert,\qquad \omega_m=2\sqrt{\frac{\beta}{m}}

∣dωdq∣=a2ωm2−ω2\left\lvert\frac{\mathrm d\omega}{\mathrm dq}\right\rvert =\frac a2\sqrt{\omega_m^2-\omega^2}

重要公式:一维单原子链模式密度

ρ(ω)=2Nπ(ωm2−ω2)−1/2,0≤ω<ωm\rho(\omega)=\frac{2N}{\pi}\left(\omega_m^2-\omega^2\right)^{-1/2}, \qquad 0\le\omega<\omega_m

其中 ρ(ω)\rho(\omega) 的单位为 s/rad\mathrm{s/rad};NN 为原子数;ωm\omega_m 为截止频率(rad/s\mathrm{rad/s});ω\omega 为角频率(rad/s\mathrm{rad/s})。适用条件:一维单原子链、最近邻、简谐近似。

推导:

ρ(ω)=∣2×L2πdqdω∣=Lπ∣dqdω∣=Naπ2aωm2−ω2=2Nπωm2−ω2\rho(\omega) =\left\lvert2\times\frac{L}{2\pi}\frac{\mathrm dq}{\mathrm d\omega}\right\rvert =\frac{L}{\pi}\left\lvert\frac{\mathrm dq}{\mathrm d\omega}\right\rvert =\frac{Na}{\pi}\frac{2}{a\sqrt{\omega_m^2-\omega^2}} =\frac{2N}{\pi\sqrt{\omega_m^2-\omega^2}}

其中 L=NaL=Na 为链长(m\mathrm m)。因子 22 来自 ±q\pm q 两支。

低维/三维模式密度表:

维度波矢密度等频面波矢数目模式密度
一维L/(2π)L/(2\pi)两点2 dq2\,\mathrm dqLπ∣1∇qω∣\dfrac{L}{\pi}\left\lvert\dfrac{1}{\nabla_q\omega}\right\rvert
二维S/(2π)2S/(2\pi)^2圆的周长2πq dq2\pi q\,\mathrm dqSq2π∣1∇qω∣\dfrac{Sq}{2\pi}\left\lvert\dfrac{1}{\nabla_q\omega}\right\rvert
三维Vc/(2π)3V_c/(2\pi)^3球面积4πq2 dq4\pi q^2\,\mathrm dq4πq2Vc(2π)3∣1∇qω∣\dfrac{4\pi q^2V_c}{(2\pi)^3}\left\lvert\dfrac{1}{\nabla_q\omega}\right\rvert

表中 LL 为长度(m\mathrm m),SS 为面积(m2\mathrm{m^2}),VcV_c 为体积(m3\mathrm{m^3}),qq 为波矢(m−1\mathrm{m^{-1}})。注意:表中结果是关于 qq 的函数,模式密度最终要以 ω\omega 为单位,必须通过色散关系转换。

3.5.6 德拜模型 ​

模型假设:

  1. 晶体视为连续介质,格波视为弹性波;
  2. 有两种振动类型:一个纵波、两个横波;
  3. 晶格振动频率在 0∼ωD0\sim\omega_{\mathrm D} 之间,ωD\omega_{\mathrm D} 为德拜频率。

色散关系:

ω=vq\omega=vq

其中 vv 为波速(m/s\mathrm{m/s}),∇qω=v\nabla_{\mathbf q}\omega=v。

由一般模式密度公式:

ρα(ω)=Vc(2π)3∫Sωdsv=Vc(2π)34πq2v=Vc2π2ω2v3\rho_\alpha(\omega)=\frac{V_c}{(2\pi)^3}\int_{S_\omega}\frac{\mathrm ds}{v} =\frac{V_c}{(2\pi)^3}\frac{4\pi q^2}{v} =\frac{V_c}{2\pi^2}\frac{\omega^2}{v^3}

重要公式:德拜模型的模式密度

ρ(ω)=3Vc2π2v3ω2=9nNωD3ω2,0≤ω≤ωD\rho(\omega)=\frac{3V_c}{2\pi^2v^3}\omega^2 =\frac{9nN}{\omega_{\mathrm D}^3}\omega^2,\qquad 0\le\omega\le\omega_{\mathrm D}

其中 ρ(ω)\rho(\omega) 的单位为 s/rad\mathrm{s/rad};VcV_c 为晶体体积(m3\mathrm{m^3});vv 为平均波速(m/s\mathrm{m/s});nNnN 为原子总数;ωD\omega_{\mathrm D} 为德拜频率(rad/s\mathrm{rad/s})。适用条件:连续介质近似、ω=vq\omega=vq、三支声学波(11 纵 22 横)。

德拜频率的确定:由归一化条件

∫0ωD9nNωD3ω2 dω=3nN\int_0^{\omega_{\mathrm D}}\frac{9nN}{\omega_{\mathrm D}^3}\omega^2\,\mathrm d\omega=3nN

自动满足。若用 ρ=3Vcω2/(2π2v3)\rho=3V_c\omega^2/(2\pi^2v^3),则

Vc2π2v3ωD33=3nN⟹ωD3=6π2nNVcv3\frac{V_c}{2\pi^2v^3}\frac{\omega_{\mathrm D}^3}{3}=3nN \quad\Longrightarrow\quad \omega_{\mathrm D}^3=6\pi^2\frac{nN}{V_c}v^3

其中 nN/VcnN/V_c 为原子数密度(m−3\mathrm{m^{-3}})。补充:若记原子数密度为 na=nN/Vcn_{\mathrm a}=nN/V_c,则 ωD=v(6π2na)1/3\omega_{\mathrm D}=v(6\pi^2n_{\mathrm a})^{1/3}。

德拜温度:

θD=ℏωDkB\theta_{\mathrm D}=\frac{\hbar\omega_{\mathrm D}}{k_{\mathrm B}}

其中 θD\theta_{\mathrm D} 的单位为 K\mathrm K。

重要公式:德拜模型的热容

CV=3nNkBf(θDT)C_V=3nNk_{\mathrm B}f\left(\frac{\theta_{\mathrm D}}{T}\right)

f(θDT)=3(TθD)3∫0θD/Tex(ex−1)2x4 dxf\left(\frac{\theta_{\mathrm D}}{T}\right) =3\left(\frac{T}{\theta_{\mathrm D}}\right)^3 \int_0^{\theta_{\mathrm D}/T} \frac{e^x}{(e^x-1)^2}x^4\,\mathrm dx

其中 CVC_V 的单位为 J/K\mathrm{J/K};nNnN 为原子总数;θD\theta_{\mathrm D} 为德拜温度(K\mathrm K);x=ℏω/(kBT)x=\hbar\omega/(k_{\mathrm B}T) 为无量纲积分变量。适用条件:连续介质近似、三支声学波、TT 不接近熔化温度。

高低温极限:

  • 高温 T≫θDT\gg\theta_{\mathrm D}:x≪1x\ll1,ex≈1+xe^x\approx1+x,被积函数 ex(ex−1)2x4≈x2\dfrac{e^x}{(e^x-1)^2}x^4\approx x^2,故

    f≈3(TθD)3∫0θD/Tx2 dx=3(TθD)3⋅13(θDT)3=1,CV≈3nNkBf\approx3\left(\frac{T}{\theta_{\mathrm D}}\right)^3\int_0^{\theta_{\mathrm D}/T}x^2\,\mathrm dx =3\left(\frac{T}{\theta_{\mathrm D}}\right)^3\cdot\frac13\left(\frac{\theta_{\mathrm D}}{T}\right)^3=1 ,\qquad C_V\approx3nNk_{\mathrm B}

    与实验高温值一致。
  • 低温 T≪θDT\ll\theta_{\mathrm D}:积分上限趋于无穷,∫0∞ex(ex−1)2x4 dx=4π415\displaystyle\int_0^\infty\frac{e^x}{(e^x-1)^2}x^4\,\mathrm dx=\dfrac{4\pi^4}{15},故

    CV≈3nNkB⋅3(TθD)3⋅4π415=12π45nNkB(TθD)3C_V\approx3nNk_{\mathrm B}\cdot 3\left(\frac{T}{\theta_{\mathrm D}}\right)^3\cdot\frac{4\pi^4}{15} =\frac{12\pi^4}{5}nNk_{\mathrm B} \left(\frac{T}{\theta_{\mathrm D}}\right)^3

    即德拜 T3T^3 定律。

不同维度的德拜结果:

维度模式密度高温热容低温热容
一维Lπv\dfrac{L}{\pi v}nNkBnNk_{\mathrm B}ATAT
二维Sω2πv2\dfrac{S\omega}{2\pi v^2}2nNkB2nNk_{\mathrm B}AT2AT^2
三维Vcω22π2v3\dfrac{V_c\omega^2}{2\pi^2v^3}3nNkB3nNk_{\mathrm B}AT3AT^3

其中 AA 为与材料有关的常数;nNnN 为原子总数。低维系统在低温下的热容不再按 T3T^3 变化。

易错点:

  • 爱因斯坦模型的高温极限是 3nNkB3nNk_{\mathrm B},低温极限是指数衰减,不是 T3T^3。
  • 德拜模型的高温极限是 3nNkB3nNk_{\mathrm B},低温极限是 T3T^3。
  • 模式密度归一化给出的是 3nN3nN(总模式数),不是 NN。
  • 德拜模型的 ω=vq\omega=vq 只适用于长波(声学支);光学支不包含在简单德拜模型中。
  • 德拜温度 θD\theta_{\mathrm D} 与爱因斯坦温度 θE\theta_{\mathrm E} 的定义都是 ℏω/kB\hbar\omega/k_{\mathrm B},但物理含义不同。

3.6 易错点与常考题型 ​

易错点

  1. 一维单原子链色散关系写成 ω=2βm(1−cos⁡qa)\omega=\frac{2\beta}{m}(1-\cos qa),漏开方或漏因子。
  2. 把 qq 的范围写成 −π/a≤q≤π/a-\pi/a\le q\le\pi/a 的长度当成 2π/a2\pi/a 以外的范围;实际长度就是一个倒格矢。
  3. 双原子链的晶格常数是 2a2a,布里渊区边界是 ±π/(2a)\pm\pi/(2a)。
  4. 声学支与光学支的命名混淆;声学支 q=0q=0 时 ω=0\omega=0,光学支 q=0q=0 时 ω≠0\omega\ne0。
  5. 长声学波振幅比约为 11(同相),长光学波振幅比约为 −M/m-M/m(反相)。
  6. 三维晶格支数写成 3n3n 支声学波;应为 33 支声学、3(n−1)3(n-1) 支光学。
  7. 声子能量写成 ℏω/2\hbar\omega/2;单声子能量是 ℏω\hbar\omega,零点能才是 ℏω/2\hbar\omega/2。
  8. 爱因斯坦模型低温热容误写成 T3T^3;它是指数衰减。
  9. 德拜模型低温热容误写成指数衰减;它是 T3T^3。
  10. 模式密度归一化写成 NN 而不是 3nN3nN。

常考题型

  1. 由运动方程推导一维单原子链色散关系,求 ωmax⁡\omega_{\max}、声速、长波/短波极限。
  2. 用玻恩—卡门条件求 qq 的取值个数与间距。
  3. 推导一维双原子链色散关系,讨论声学支/光学支、频率间隙、振幅关系。
  4. 判断给定晶体的格波支数、声学支数、光学支数、总模式数。
  5. 计算声子能量、准动量、平均声子数。
  6. 由模式密度定义推导一维链或德拜模型的 ρ(ω)\rho(\omega)。
  7. 用爱因斯坦模型/德拜模型计算 CVC_V,讨论高低温极限。
  8. 解释爱因斯坦模型低温失效的原因。

第 4 章 能带论 ​

4.1 能带论的基础与三种近似 ​

电子共有化:价电子不再束缚于某一原子,可以在不同原子的相似轨道间转移,从而在整个固体中运动。

能带的起源:孤立原子有分立能级(1s,2s,2p,⋯1s,2s,2p,\cdots)。NN 个原子组成晶体时,由于泡利不相容原理,原来相同的能级分裂成 NN 个能量接近的新能级。NN 很大时,相邻能级能量差仅约 10−22 eV10^{-22}\ \mathrm{eV},几乎连续,形成能量范围约 eV\mathrm{eV} 的能带。

能带的一般规律:

  1. 外层电子共有化程度显著,能带较宽;内层电子轨道重叠少,能带较窄。
  2. 点阵间距越小,能带越宽。
  3. 两能带可能重叠。
  4. 能带在布里渊区边界处产生能隙;禁带宽度取决于周期势场变化的剧烈程度。
  5. 晶体中电子波函数是布洛赫函数,兼有共有化运动和绕原子核运动的特点。

三种近似:

  1. 绝热近似:原子核质量远大于电子,电子运动速度远快于核;把核视为不动,电子在固定核势场中运动。
  2. 单电子近似:把多电子体系化为单电子问题;每个电子在离子实和其他电子形成的平均势场中运动。
  3. 周期场近似:晶体中大量原子规则排列,形成周期性势场 V(r)=V(r+R)V(\mathbf r)=V(\mathbf r+\mathbf R)。

两大模型:

  • 近自由电子近似:价电子能量高、晶体势场较弱,电子行为接近自由电子;把周期性势场的微弱起伏作为微扰。
  • 紧束缚近似:内层电子能量低、晶体势场较强,电子基本围绕原子核运动;把其他原子的作用作为微扰,用孤立原子轨道的线性组合(LCAO)构成晶体波函数。

4.2 布洛赫定理 ​

4.2.1 定理内容 ​

晶体中单电子哈密顿量

H^=−ℏ22m∇2+V(r)\hat H=-\frac{\hbar^2}{2m}\nabla^2+V(\mathbf r)

其中 H^\hat H 的单位为 J\mathrm J;mm 为电子质量(kg\mathrm{kg});ℏ\hbar 为约化普朗克常量(J⋅s\mathrm{J\cdot s});V(r)V(\mathbf r) 为周期性势场(J\mathrm J),满足

V(r+Rn)=V(r)V(\mathbf r+\mathbf R_n)=V(\mathbf r)

Rn\mathbf R_n 为任意正格矢(m\mathrm m)。

重要公式:布洛赫定理

ψk(r)=eik⋅ruk(r),uk(r+Rn)=uk(r)\psi_{\mathbf k}(\mathbf r)=e^{i\mathbf k\cdot\mathbf r}u_{\mathbf k}(\mathbf r), \qquad u_{\mathbf k}(\mathbf r+\mathbf R_n)=u_{\mathbf k}(\mathbf r)

等价形式:

ψk(r+Rn)=eik⋅Rnψk(r)\psi_{\mathbf k}(\mathbf r+\mathbf R_n)=e^{i\mathbf k\cdot\mathbf R_n}\psi_{\mathbf k}(\mathbf r)

其中 ψk\psi_{\mathbf k} 为电子波函数;k\mathbf k 为简约波矢(m−1\mathrm{m^{-1}});uku_{\mathbf k} 为具有晶格周期性的调制函数;Rn\mathbf R_n 为正格矢(m\mathrm m)。适用条件:单电子近似、周期性势场。

物理意义:

  • eik⋅re^{i\mathbf k\cdot\mathbf r} 为平面波,反映电子共有化特征;
  • uk(r)u_{\mathbf k}(\mathbf r) 为周期性函数,反映电子受周期势场束缚的特征;
  • ℏk\hbar\mathbf k 称为准动量,k\mathbf k 是标志电子状态的量子数;不同 k\mathbf k 表示不同状态、不同能量,其物理意义是不同原胞之间波函数的位相差。

简约波矢的取值范围:

k=l1N1b1+l2N2b2+l3N3b3\mathbf k=\frac{l_1}{N_1}\mathbf b_1+\frac{l_2}{N_2}\mathbf b_2+\frac{l_3}{N_3}\mathbf b_3

其中 lil_i 为整数;NiN_i 为三个方向的元胞数;bi\mathbf b_i 为倒格子基矢(m−1\mathrm{m^{-1}})。k\mathbf k 通常取在简约布里渊区(第一布里渊区)内。

能带的周期性:

E(k)=E(k+Kh)E(\mathbf k)=E(\mathbf k+\mathbf K_h)

其中 Kh\mathbf K_h 为任意倒格矢(m−1\mathrm{m^{-1}})。k\mathbf k 与 k+Kh\mathbf k+\mathbf K_h 描述同一电子态。

4.2.2 布洛赫定理的证明 ​

第一步:平移算符及其性质。定义平移算符

Tαf(r)=f(r+aα),α=1,2,3T_{\alpha}f(\mathbf r)=f(\mathbf r+\mathbf a_\alpha),\qquad \alpha=1,2,3

其中 aα\mathbf a_\alpha 为三个正格基矢(m\mathrm m)。

性质:

  1. TαTβ=TβTαT_\alpha T_\beta=T_\beta T_\alpha;
  2. TαTβ=Tα+βT_\alpha T_\beta=T_{\alpha+\beta};
  3. [H,Tα]=0[H,T_\alpha]=0。

证明性质 3:

HTαψ(r)=Hψ(r+aα)=[−ℏ22m∇2+V(r+aα)]ψ(r+aα)HT_\alpha\psi(\mathbf r) =H\psi(\mathbf r+\mathbf a_\alpha) =\left[-\frac{\hbar^2}{2m}\nabla^2+V(\mathbf r+\mathbf a_\alpha)\right]\psi(\mathbf r+\mathbf a_\alpha)

由 V(r+aα)=V(r)V(\mathbf r+\mathbf a_\alpha)=V(\mathbf r),且 ∇2\nabla^2 对 r\mathbf r 与 r+aα\mathbf r+\mathbf a_\alpha 相同,

HTαψ(r)=[−ℏ22m∇2+V(r)]ψ(r+aα)=TαHψ(r)HT_\alpha\psi(\mathbf r) =\left[-\frac{\hbar^2}{2m}\nabla^2+V(\mathbf r)\right]\psi(\mathbf r+\mathbf a_\alpha) =T_\alpha H\psi(\mathbf r)

故 HTα=TαHHT_\alpha=T_\alpha H。

第二步:确定本征值形式。由于 [H,Tα]=0[H,T_\alpha]=0,HH 与 TαT_\alpha 有共同本征函数:

Hψ=Eψ,Tαψ=λαψH\psi=E\psi,\qquad T_\alpha\psi=\lambda_\alpha\psi

由 TαTβ=Tα+βT_\alpha T_\beta=T_{\alpha+\beta} 得

λα+β=λαλβ\lambda_{\alpha+\beta}=\lambda_\alpha\lambda_\beta

只有指数函数满足这一关系,故可设

λα=eik⋅aα\lambda_\alpha=e^{i\mathbf k\cdot\mathbf a_\alpha}

其中 k\mathbf k 为待定矢量。推广到任意正格矢 Rm=m1a1+m2a2+m3a3\mathbf R_m=m_1\mathbf a_1+m_2\mathbf a_2+m_3\mathbf a_3:

TRmψ(r)=ψ(r+Rm)=eik⋅Rmψ(r)T_{\mathbf R_m}\psi(\mathbf r)=\psi(\mathbf r+\mathbf R_m) =e^{i\mathbf k\cdot\mathbf R_m}\psi(\mathbf r)

第三步:由边界条件确定 k\mathbf k。设晶体在三个方向分别有 N1,N2,N3N_1,N_2,N_3 个原胞,采用玻恩—卡门周期性边界条件:

ψ(r+N1a1)=ψ(r),ψ(r+N2a2)=ψ(r),ψ(r+N3a3)=ψ(r)\psi(\mathbf r+N_1\mathbf a_1)=\psi(\mathbf r),\quad \psi(\mathbf r+N_2\mathbf a_2)=\psi(\mathbf r),\quad \psi(\mathbf r+N_3\mathbf a_3)=\psi(\mathbf r)

代入平移本征值:

eiN1k⋅a1=1,eiN2k⋅a2=1,eiN3k⋅a3=1e^{iN_1\mathbf k\cdot\mathbf a_1}=1,\quad e^{iN_2\mathbf k\cdot\mathbf a_2}=1,\quad e^{iN_3\mathbf k\cdot\mathbf a_3}=1

故

Nik⋅ai=2πli,li∈ZN_i\mathbf k\cdot\mathbf a_i=2\pi l_i,\qquad l_i\in\mathbb Z

利用 ai⋅bj=2πδij\mathbf a_i\cdot\mathbf b_j=2\pi\delta_{ij},得

k=l1N1b1+l2N2b2+l3N3b3\mathbf k=\frac{l_1}{N_1}\mathbf b_1+\frac{l_2}{N_2}\mathbf b_2+\frac{l_3}{N_3}\mathbf b_3

第四步:化为布洛赫形式。令

uk(r)=e−ik⋅rψk(r)u_{\mathbf k}(\mathbf r)=e^{-i\mathbf k\cdot\mathbf r}\psi_{\mathbf k}(\mathbf r)

则

uk(r+Rn)=e−ik⋅(r+Rn)ψk(r+Rn)=e−ik⋅re−ik⋅Rneik⋅Rnψk(r)=uk(r)u_{\mathbf k}(\mathbf r+\mathbf R_n) =e^{-i\mathbf k\cdot(\mathbf r+\mathbf R_n)}\psi_{\mathbf k}(\mathbf r+\mathbf R_n) =e^{-i\mathbf k\cdot\mathbf r}e^{-i\mathbf k\cdot\mathbf R_n} e^{i\mathbf k\cdot\mathbf R_n}\psi_{\mathbf k}(\mathbf r) =u_{\mathbf k}(\mathbf r)

故 uku_{\mathbf k} 具有晶格周期性,ψk=eik⋅ruk\psi_{\mathbf k}=e^{i\mathbf k\cdot\mathbf r}u_{\mathbf k}。

易错点:

  • 布洛赫定理说的是「波函数乘以平面波因子后具有晶格周期性」,不是「波函数本身具有晶格周期性」。
  • k\mathbf k 是简约波矢,不是真实动量;ℏk\hbar\mathbf k 是准动量。
  • k\mathbf k 与 k+Kh\mathbf k+\mathbf K_h 描述同一状态,能量相同;不要把 k\mathbf k 的唯一取值与状态唯一对应混淆。
  • 证明中 [H,Tα]=0[H,T_\alpha]=0 的关键是势场的周期性;若势场非周期,布洛赫定理不成立。

4.3 近自由电子近似 ​

4.3.1 零级近似 ​

取周期性势场的平均值为零(或并入能量零点),零级哈密顿量

H^0=−ℏ22md2dx2+V0\hat H_0=-\frac{\hbar^2}{2m}\frac{\mathrm d^2}{\mathrm dx^2}+V_0

零级能量与本征函数:

Ek0=ℏ2k22m+V0,ψk0(x)=1LeikxE_k^0=\frac{\hbar^2k^2}{2m}+V_0,\qquad \psi_k^0(x)=\frac{1}{\sqrt L}e^{ikx}

其中 Ek0E_k^0 的单位为 J\mathrm J;kk 为波矢(m−1\mathrm{m^{-1}});mm 为电子质量(kg\mathrm{kg});L=NaL=Na 为晶体长度(m\mathrm m);NN 为原胞数;aa 为晶格常数(m\mathrm m);V0V_0 为平均势能(J\mathrm J)。

由周期性边界条件:

k=2πlNa,l∈Zk=\frac{2\pi l}{Na},\qquad l\in\mathbb Z

4.3.2 周期势场的傅里叶展开 ​

V(x)=V0+∑n≠0Vnei2πnaxV(x)=V_0+\sum_{n\ne0}V_n e^{i\frac{2\pi n}{a}x}

Vn=1a∫−a/2a/2V(x)e−i2πnax dxV_n=\frac1a\int_{-a/2}^{a/2}V(x)e^{-i\frac{2\pi n}{a}x}\,\mathrm dx

其中 VnV_n 为第 nn 个傅里叶系数(J\mathrm J);aa 为晶格常数(m\mathrm m);V(x)V(x) 为周期势场(J\mathrm J)。

4.3.3 非简并微扰 ​

一级能量修正:

Ek(1)=⟨k∣V−V0∣k⟩=1L∫0LV(x) dx=V0E_k^{(1)}=\langle k\lvert V-V_0\rvert k\rangle =\frac1L\int_0^L V(x)\,\mathrm dx=V_0

若将 V0V_0 取为零,则 Ek(1)=0E_k^{(1)}=0。

二级能量修正:

重要公式:近自由电子的二级微扰能量

Ek=Ek0+∑n≠0∣Vn∣2Ek0−Ek−2πna0E_k=E_k^0+\sum_{n\ne0} \frac{\lvert V_n\rvert^2} {E_k^0-E_{k-\frac{2\pi n}{a}}^0}

Ek=ℏ2k22m+V0+∑n≠0∣Vn∣2ℏ22m[k2−(k−2πna)2]E_k=\frac{\hbar^2k^2}{2m}+V_0 +\sum_{n\ne0} \frac{\lvert V_n\rvert^2} {\frac{\hbar^2}{2m}\left[k^2-\left(k-\frac{2\pi n}{a}\right)^2\right]}

其中 EkE_k 为电子能量(J\mathrm J);kk 为波矢(m−1\mathrm{m^{-1}});VnV_n 为周期势场的傅里叶系数(J\mathrm J);mm 为电子质量(kg\mathrm{kg});aa 为晶格常数(m\mathrm m)。适用条件:非简并微扰,即 kk 远离布里渊区边界,Ek0≠Ek′0E_k^0\ne E_{k'}^0。

微扰后的波函数:

重要公式:近自由电子波函数

ψk(x)=1Leikx+1L∑n≠0VnEk0−Ek−2πna0ei(k−2πna)x\psi_k(x)=\frac{1}{\sqrt L}e^{ikx} +\frac{1}{\sqrt L}\sum_{n\ne0} \frac{V_n}{E_k^0-E_{k-\frac{2\pi n}{a}}^0} e^{i\left(k-\frac{2\pi n}{a}\right)x}

其中 ψk\psi_k 为电子波函数;L=NaL=Na 为晶体长度(m\mathrm m);kk 为波矢(m−1\mathrm{m^{-1}});VnV_n 为势场傅里叶系数(J\mathrm J);Ek0E_k^0 为零级能量(J\mathrm J);aa 为晶格常数(m\mathrm m)。第一项为波矢为 kk 的前进平面波,第二项为周期势场作用产生的散射波。适用条件:非简并微扰、kk 远离布里渊区边界。

适用性:当

k−2πna=−k⟹k=nπak-\frac{2\pi n}{a}=-k \quad\Longrightarrow\quad k=\frac{n\pi}{a}

即 kk 位于布里渊区边界时,Ek0=Ek′0E_k^0=E_{k'}^0,分母为零,非简并微扰失效,必须用简并微扰。

4.3.4 简并微扰与能带分裂 ​

在布里渊区边界 k=nπ/ak=n\pi/a 附近,取

k=−nπa+Δ,k′=nπa+Δk=-\frac{n\pi}{a}+\Delta,\qquad k'=\frac{n\pi}{a}+\Delta

其中 Δ\Delta 为小量(m−1\mathrm{m^{-1}})。零级波函数取两个简并态的线性组合:

ψ=Aψk0+Bψk′0\psi=A\psi_k^0+B\psi_{k'}^0

代入薛定谔方程,左乘 ψk0∗\psi_k^{0*}、ψk′0∗\psi_{k'}^{0*} 并积分,得

{(Ek0−E)A+VnB=0Vn∗A+(Ek′0−E)B=0\begin{cases} (E_k^0-E)A+V_nB=0\\[4pt] V_n^*A+(E_{k'}^0-E)B=0 \end{cases}

非零解条件:

∣Ek0−EVnVn∗Ek′0−E∣=0\begin{vmatrix} E_k^0-E & V_n\\ V_n^* & E_{k'}^0-E \end{vmatrix}=0

重要公式:简并微扰的能量分裂

E±=12(Ek0+Ek′0)±14(Ek0−Ek′0)2+∣Vn∣2E_{\pm}=\frac12\left(E_k^0+E_{k'}^0\right) \pm\sqrt{\frac14\left(E_k^0-E_{k'}^0\right)^2+\lvert V_n\rvert^2}

其中 E±E_\pm 为分裂后的两个能级(J\mathrm J);Ek0,Ek′0E_k^0,E_{k'}^0 为两个简并零级能量(J\mathrm J);VnV_n 为势场傅里叶系数(J\mathrm J)。适用条件:kk 在布里渊区边界附近、两态近简并。

在边界上 Δ=0\Delta=0,Ek0=Ek′0=TnE_k^0=E_{k'}^0=T_n,其中

Tn=ℏ22m(nπa)2T_n=\frac{\hbar^2}{2m}\left(\frac{n\pi}{a}\right)^2

为自由电子在 k=nπ/ak=n\pi/a 处的动能(J\mathrm J)。于是

E+=Tn+∣Vn∣,E−=Tn−∣Vn∣E_+=T_n+\lvert V_n\rvert,\qquad E_-=T_n-\lvert V_n\rvert

重要公式:禁带宽度

Eg=E+−E−=2∣Vn∣E_g=E_+-E_-=2\lvert V_n\rvert

其中 EgE_g 为禁带宽度(J\mathrm J 或 eV\mathrm{eV});VnV_n 为周期势场第 nn 个傅里叶系数。适用条件:近自由电子近似、布里渊区边界处。

物理图像:在布里渊区边界,前进波与布拉格反射波形成驻波;电子云分别集中在离子实之间(能量较低)和离子实上(能量较高),两种驻波的能量差就是禁带宽度。周期势场越强,∣Vn∣\lvert V_n\rvert 越大,禁带越宽。

易错点:

  • 非简并微扰在布里渊区边界失效,必须换用简并微扰。
  • 禁带宽度是 2∣Vn∣2\lvert V_n\rvert,不是 ∣Vn∣\lvert V_n\rvert。
  • 简并微扰中两个零级态是 kk 与 k′=k−2πn/ak'=k-2\pi n/a(或 k+2πn/ak+2\pi n/a),不是任意两个态。
  • 能带分裂后,原来能量较高的态升高、较低的态降低,称为能级「排斥」。

4.4 能带的三种表示方式 ​

  1. 扩展能区图式:把每个能带画在各自的 kk 区间内。第一能带 k∈(−π/a,π/a)k\in(-\pi/a,\pi/a),第二能带 k∈(π/a,2π/a)k\in(\pi/a,2\pi/a) 及 (−2π/a,−π/a)(-2\pi/a,-\pi/a),依此类推。
  2. 周期能区图式:利用 E(k)=E(k+Kh)E(\mathbf k)=E(\mathbf k+\mathbf K_h),把任意一条能量曲线通过倒格矢从一个布里渊区移到其他布里渊区,在每个布里渊区画出所有能带,构成 kk 空间中能量分布的完整图像。
  3. 简约能区图式:把每个能带都折叠到第一布里渊区内;每个能带在简约布里渊区都有各自的图像。简约布里渊区中的一个状态由「属于哪一个能带」和「简约波矢 k\mathbf k」共同标志。

折叠规则:

k′=k±2πak'=k\pm\frac{2\pi}{a}

其中 k′k' 为折叠后的简约波矢(m−1\mathrm{m^{-1}});kk 为原波矢(m−1\mathrm{m^{-1}});aa 为晶格常数(m\mathrm m)。通过平移整数个倒格矢,把 kk 移入第一布里渊区。

易错点:

  • 三种图式描述的是同一套 E(k)E(\mathbf k),只是 kk 的取法不同,不要把它们当成不同的物理。
  • 简约能区图中,同一 kk 点可能有多个能带值,必须标明能带指标。
  • 扩展图式中不同能带的 kk 区间不重叠;周期图式中每个布里渊区都画所有能带。

4.5 带隙条件与布拉格反射 ​

一维带隙条件:

k2=(k+2πna)2⟹k=nπak^2=\left(k+\frac{2\pi n}{a}\right)^2 \quad\Longrightarrow\quad k=\frac{n\pi}{a}

其中 kk 为波矢(m−1\mathrm{m^{-1}});aa 为晶格常数(m\mathrm m);nn 为整数。

三维带隙条件:

重要公式:三维带隙条件

∣k∣2=∣k+Gn∣2⟺2k⋅Gn+Gn2=0⟺k⋅Gn=−12Gn2\lvert\mathbf k\rvert^2=\lvert\mathbf k+\mathbf G_n\rvert^2 \quad\Longleftrightarrow\quad 2\mathbf k\cdot\mathbf G_n+\mathbf G_n^2=0 \quad\Longleftrightarrow\quad \mathbf k\cdot\mathbf G_n=-\frac12\mathbf G_n^2

其中 k\mathbf k 为波矢(m−1\mathrm{m^{-1}});Gn\mathbf G_n 为倒格矢(m−1\mathrm{m^{-1}})。适用条件:近自由电子近似、布里渊区边界附近。

布拉格反射条件:

2dsin⁡θ=nλ2d\sin\theta=n\lambda

其中 dd 为晶面间距(m\mathrm m);θ\theta 为入射角(与晶面的夹角,rad\mathrm{rad});λ\lambda 为波长(m\mathrm m);nn 为反射级数(整数)。

由 d=2π/∣Gn∣d=2\pi/\lvert\mathbf G_n\rvert 与 ∣Gn∣=2ksin⁡θ\lvert\mathbf G_n\rvert=2k\sin\theta,可得

2dsin⁡θ=nλ2d\sin\theta=n\lambda

即布拉格反射条件与带隙条件是等价的:布里渊区边界处满足布拉格反射条件,前进波与反射波相干叠加形成驻波,产生能隙。

易错点:

  • 三维带隙条件是 k⋅Gn=−12Gn2\mathbf k\cdot\mathbf G_n=-\tfrac12G_n^2,不是 k⋅Gn=0\mathbf k\cdot\mathbf G_n=0。
  • 布拉格公式中的 θ\theta 是入射光与晶面的夹角,不是与法线的夹角。
  • 带隙出现在布里渊区边界,不是布里渊区中心;中心处能量与自由电子接近。

4.6 布里渊区 ​

定义:在倒格子空间中,以某一格点为原点,作所有倒格矢 G\mathbf G 的垂直平分面。这些平面把倒易空间分割成许多包围原点的多面体,其中离原点最近的多面体称为第一布里渊区;离原点次近的多面体与第一布里渊区表面围成的区域称为第二布里渊区,依此类推。

二维正方格子:

倒格子基矢

b1=2πai,b2=2πaj\mathbf b_1=\frac{2\pi}{a}\mathbf i,\qquad \mathbf b_2=\frac{2\pi}{a}\mathbf j

倒格矢

G=2πa(h1i+h2j)\mathbf G= \frac{2\pi}{a}(h_1\mathbf i+h_2\mathbf j)

第一布里渊区为边长 2π/a2\pi/a 的正方形。特殊点:

Γ: k=(0,0),X: k=(πa,0),M: k=(πa,πa)\Gamma:\ \mathbf k=(0,0),\qquad X:\ \mathbf k=\left(\frac{\pi}{a},0\right),\qquad M:\ \mathbf k=\left(\frac{\pi}{a},\frac{\pi}{a}\right)

简单立方格子:第一布里渊区为边长 2π/a2\pi/a 的立方体。特殊点:

Γ: (0,0,0),X: (πa,0,0),M: (πa,πa,0),R: (πa,πa,πa)\Gamma:\ (0,0,0),\quad X:\ \left(\frac{\pi}{a},0,0\right),\quad M:\ \left(\frac{\pi}{a},\frac{\pi}{a},0\right),\quad R:\ \left(\frac{\pi}{a},\frac{\pi}{a},\frac{\pi}{a}\right)

面心立方格子:

  • 倒格子为体心立方;最近邻倒格点有 88 个,对应 88 个中垂面,围成八面体;但 66 个顶角被次近邻倒格点的中垂面截去,故第一布里渊区为十四面体(截角八面体)。
  • 第一布里渊区体积:

    VBZ=4(2πa)3V_{\mathrm{BZ}}=4\left(\frac{2\pi}{a}\right)^3

    其中 VBZV_{\mathrm{BZ}} 的单位为 m−3\mathrm{m^{-3}};aa 为晶格常数(m\mathrm m)。

体心立方格子:

  • 倒格子为面心立方;最近邻倒格点有 1212 个,对应 1212 个中垂面,围成十二面体(菱形十二面体);次近邻中垂面不切割它。
  • 第一布里渊区体积:

    VBZ=2(2πa)3V_{\mathrm{BZ}}=2\left(\frac{2\pi}{a}\right)^3

    其中 VBZV_{\mathrm{BZ}} 的单位为 m−3\mathrm{m^{-3}};aa 为晶格常数(m\mathrm m)。

易错点:

  • 布里渊区是倒格子空间中的区域,不是正空间中的晶胞。
  • fcc 的第一布里渊区是十四面体,不是八面体;bcc 的是十二面体,不是立方体。
  • 布里渊区体积等于倒格子原胞体积 (2π)3/Ω(2\pi)^3/\Omega;fcc/bcc 的具体值要用对应的 Ω\Omega 计算。
  • 特殊点符号 Γ,X,M,R\Gamma,X,M,R 的位置要记清;不同结构的特殊点不同。

4.7 紧束缚近似 ​

4.7.1 基本思想 ​

  • 晶格中原子间距较大,原子附近的电子受自身束缚较紧,不易共有化;
  • 在原子附近,电子行为接近孤立原子中的电子,晶体波函数接近孤立原子波函数;
  • 以孤立原子波函数为零级近似,把其他原子场的作用看成微扰;
  • 实质:把原子间相互作用看成微扰的简并微扰方法,微扰后的状态是 NN 个简并态的线性组合,即原子轨道线性组合法(LCAO)。

4.7.2 模型与微扰计算 ​

孤立原子本征方程:

[−ℏ22m∇2+V(r−Rm)]ϕi(r−Rm)=εiϕi(r−Rm)\left[-\frac{\hbar^2}{2m}\nabla^2+V(\mathbf r-\mathbf R_m)\right]\phi_i(\mathbf r-\mathbf R_m) =\varepsilon_i\phi_i(\mathbf r-\mathbf R_m)

其中 ϕi\phi_i 为孤立原子第 ii 个轨道波函数;εi\varepsilon_i 为孤立原子能级(J\mathrm J);Rm\mathbf R_m 为格点位矢(m\mathrm m)。

晶体中电子波动方程:

[−ℏ22m∇2+U(r)]ψ(r)=Eψ(r)\left[-\frac{\hbar^2}{2m}\nabla^2+U(\mathbf r)\right]\psi(\mathbf r)=E\psi(\mathbf r)

其中

U(r)=∑mV(r−Rm)U(\mathbf r)=\sum_m V(\mathbf r-\mathbf R_m)

为晶体周期性势场(J\mathrm J)。

令

ΔV(r)=U(r)−V(r−Rm)\Delta V(\mathbf r)=U(\mathbf r)-V(\mathbf r-\mathbf R_m)

为其他所有格点对电子的作用之和,视为微扰。

晶体波函数取原子轨道线性组合:

ψi(k,r)=∑mamϕi(r−Rm)\psi_i(\mathbf k,\mathbf r)=\sum_m a_m\phi_i(\mathbf r-\mathbf R_m)

其中 ama_m 为待定系数;k\mathbf k 为波矢(m−1\mathrm{m^{-1}})。由布洛赫定理,系数应取

am=Ceik⋅Rma_m=Ce^{i\mathbf k\cdot\mathbf R_m}

其中 CC 为归一化常数。

代入波动方程,左乘 ϕi∗(r−Rn)\phi_i^*(\mathbf r-\mathbf R_n) 并积分,得

(E−εi)an=−∑mJ(Rn−Rm)am(E-\varepsilon_i)a_n=-\sum_m J(\mathbf R_n-\mathbf R_m)a_m

其中

J(Rs)=−∫ϕi∗(ξ)[U(ξ)−V(ξ−Rs)]ϕi(ξ−Rs) dξJ(\mathbf R_s)=-\int\phi_i^*(\boldsymbol\xi) \left[U(\boldsymbol\xi)-V(\boldsymbol\xi-\mathbf R_s)\right] \phi_i(\boldsymbol\xi-\mathbf R_s)\,\mathrm d\boldsymbol\xi

为跃迁积分(J\mathrm J);Rs=Rn−Rm\mathbf R_s=\mathbf R_n-\mathbf R_m。

将 am=Ceik⋅Rma_m=Ce^{i\mathbf k\cdot\mathbf R_m} 代入,得

重要公式:紧束缚近似能带

Ei(k)=εi−J0−∑RsJ(Rs)eik⋅RsE_i(\mathbf k)=\varepsilon_i-J_0-\sum_{\mathbf R_s}J(\mathbf R_s)e^{i\mathbf k\cdot\mathbf R_s}

其中 Ei(k)E_i(\mathbf k) 为第 ii 个能带的能量(J\mathrm J);εi\varepsilon_i 为孤立原子能级(J\mathrm J);J0J_0 为同一格点的积分(J\mathrm J);J(Rs)J(\mathbf R_s) 为跃迁积分(J\mathrm J);Rs\mathbf R_s 为近邻格矢(m\mathrm m);k\mathbf k 为波矢(m−1\mathrm{m^{-1}})。适用条件:紧束缚近似、只考虑最近邻、原子轨道正交归一。

只考虑最近邻:

Ei(k)=εi−J0−∑Rs∈NNJ(Rs)eik⋅RsE_i(\mathbf k)=\varepsilon_i-J_0-\sum_{\mathbf R_s\in\mathrm{NN}}J(\mathbf R_s)e^{i\mathbf k\cdot\mathbf R_s}

其中 NN\mathrm{NN} 表示最近邻;J0J_0 和最近邻跃迁积分 J1J_1 通常为正值。

4.7.3 三个立方晶格的 ss 能带 ​

简立方:最近邻 66 个,Rs=(±a,0,0),(0,±a,0),(0,0,±a)\mathbf R_s=(\pm a,0,0),(0,\pm a,0),(0,0,\pm a)。

Es(k)=εs−J0−2J1(cos⁡kxa+cos⁡kya+cos⁡kza)E_s(\mathbf k)=\varepsilon_s-J_0-2J_1\left(\cos k_xa+\cos k_ya+\cos k_za\right)

其中 J1J_1 为最近邻跃迁积分(J\mathrm J);aa 为晶格常数(m\mathrm m)。

  • 带底:k=0\mathbf k=0(Γ\Gamma 点),Emin⁡=εs−J0−6J1E_{\min}=\varepsilon_s-J_0-6J_1;
  • 带顶:k=(π/a,π/a,π/a)\mathbf k=(\pi/a,\pi/a,\pi/a)(RR 点),Emax⁡=εs−J0+6J1E_{\max}=\varepsilon_s-J_0+6J_1;
  • 带宽:ΔE=12J1\Delta E=12J_1。

面心立方:最近邻 1212 个,Rs=a2(±1,±1,0)\mathbf R_s=\frac a2(\pm1,\pm1,0) 等。

Es(k)=εs−J0−4J1[cos⁡kxa2cos⁡kya2+cos⁡kya2cos⁡kza2+cos⁡kza2cos⁡kxa2]E_s(\mathbf k)=\varepsilon_s-J_0-4J_1\left[ \cos\frac{k_xa}{2}\cos\frac{k_ya}{2} +\cos\frac{k_ya}{2}\cos\frac{k_za}{2} +\cos\frac{k_za}{2}\cos\frac{k_xa}{2} \right]

  • 带底:k=0\mathbf k=0(Γ\Gamma 点),Emin⁡=εs−J0−12J1E_{\min}=\varepsilon_s-J_0-12J_1;
  • 带顶:k=(2π/a,0,0)\mathbf k=(2\pi/a,0,0) 等,Emax⁡=εs−J0+4J1E_{\max}=\varepsilon_s-J_0+4J_1;
  • 带宽:ΔE=16J1\Delta E=16J_1。

体心立方:最近邻 88 个,Rs=a2(±1,±1,±1)\mathbf R_s=\frac a2(\pm1,\pm1,\pm1)。

Es(k)=εs−J0−8J1cos⁡kxa2cos⁡kya2cos⁡kza2E_s(\mathbf k)=\varepsilon_s-J_0-8J_1\cos\frac{k_xa}{2}\cos\frac{k_ya}{2}\cos\frac{k_za}{2}

  • 带底:k=0\mathbf k=0(Γ\Gamma 点),Emin⁡=εs−J0−8J1E_{\min}=\varepsilon_s-J_0-8J_1;
  • 带顶:k=(2π/a,0,0)\mathbf k=(2\pi/a,0,0) 等,Emax⁡=εs−J0+8J1E_{\max}=\varepsilon_s-J_0+8J_1;
  • 带宽:ΔE=16J1\Delta E=16J_1。

能带宽度与重叠:带宽正比于 J1J_1,而 J1J_1 取决于近邻原子波函数的重叠程度;重叠越多,能带越宽。原子间距越小,重叠越大,能带越宽。

易错点:

  • 紧束缚近似中 k\mathbf k 的相位因子是 eik⋅Rse^{i\mathbf k\cdot\mathbf R_s},不是 e−ik⋅Rse^{-i\mathbf k\cdot\mathbf R_s};两者只差复共轭,能量相同。
  • 最近邻数目与格矢不能漏:sc 为 66,fcc 为 1212,bcc 为 88。
  • 带宽是 12J112J_1(sc)或 16J116J_1(fcc、bcc),不是 2J12J_1。
  • J1>0J_1>0 时,ss 带底在 Γ\Gamma 点;带顶位置随结构不同。
  • 紧束缚近似适用于内层电子或窄带;近自由电子近似适用于外层价电子或宽带。

4.8 电子的速度、加速度与有效质量 ​

4.8.1 波包与电子速度 ​

把电子看成准经典粒子,用波包描述。波包群速度即电子平均速度:

重要公式:电子速度

v(k)=1ℏ∇kE(k)\mathbf v(\mathbf k)=\frac1\hbar\nabla_{\mathbf k}E(\mathbf k)

其中 v\mathbf v 为电子速度(m/s\mathrm{m/s});E(k)E(\mathbf k) 为能带能量(J\mathrm J);k\mathbf k 为波矢(m−1\mathrm{m^{-1}});ℏ\hbar 为约化普朗克常量。适用条件:准经典近似、波包尺度远大于原胞、Δk≪2π/a\Delta k\ll2\pi/a。

波包中心位置:

r=r0+v(k0)t\mathbf r=\mathbf r_0+\mathbf v(\mathbf k_0)t

其中 r\mathbf r 为波包中心位置(m\mathrm m);v(k0)\mathbf v(\mathbf k_0) 为群速度(m/s\mathrm{m/s});tt 为时间(s\mathrm s)。

4.8.2 加速度与有效质量 ​

在外力 F\mathbf F 作用下,准经典运动方程:

ℏdkdt=F\hbar\frac{\mathrm d\mathbf k}{\mathrm dt}=\mathbf F

其中 F\mathbf F 为外力(N\mathrm N);k\mathbf k 为波矢(m−1\mathrm{m^{-1}});tt 为时间(s\mathrm s)。

由 a=dv/dt\mathbf a=\mathrm d\mathbf v/\mathrm dt 与链式法则:

a=1ℏ∇k[∇kE(k)]⋅dkdt=1ℏ2∇k[∇kE(k)]⋅F\mathbf a=\frac1\hbar\nabla_{\mathbf k}\left[\nabla_{\mathbf k}E(\mathbf k)\right]\cdot\frac{\mathrm d\mathbf k}{\mathrm dt} =\frac1{\hbar^2}\nabla_{\mathbf k}\left[\nabla_{\mathbf k}E(\mathbf k)\right]\cdot\mathbf F

重要公式:电子加速度

a=1ℏ2∇k[∇kE(k)]⋅F\mathbf a=\frac1{\hbar^2}\nabla_{\mathbf k}\left[\nabla_{\mathbf k}E(\mathbf k)\right]\cdot\mathbf F

其中 a\mathbf a 为加速度(m/s2\mathrm{m/s^2});F\mathbf F 为外力(N\mathrm N);E(k)E(\mathbf k) 为能带能量(J\mathrm J);∇k∇kE\nabla_{\mathbf k}\nabla_{\mathbf k}E 为 3×33\times3 张量(J⋅m2\mathrm{J\cdot m^2})。适用条件:准经典近似、单带、外力远小于晶格周期势场。

写成矩阵:

(axayaz)=1ℏ2(∂2E∂kx2∂2E∂kx∂ky∂2E∂kx∂kz∂2E∂ky∂kx∂2E∂ky2∂2E∂ky∂kz∂2E∂kz∂kx∂2E∂kz∂ky∂2E∂kz2)(FxFyFz)\begin{pmatrix}a_x\\a_y\\a_z\end{pmatrix} =\frac1{\hbar^2} \begin{pmatrix} \dfrac{\partial^2E}{\partial k_x^2}&\dfrac{\partial^2E}{\partial k_x\partial k_y}&\dfrac{\partial^2E}{\partial k_x\partial k_z}\\[6pt] \dfrac{\partial^2E}{\partial k_y\partial k_x}&\dfrac{\partial^2E}{\partial k_y^2}&\dfrac{\partial^2E}{\partial k_y\partial k_z}\\[6pt] \dfrac{\partial^2E}{\partial k_z\partial k_x}&\dfrac{\partial^2E}{\partial k_z\partial k_y}&\dfrac{\partial^2E}{\partial k_z^2} \end{pmatrix} \begin{pmatrix}F_x\\F_y\\F_z\end{pmatrix}

重要公式:倒有效质量张量

1mαβ∗=1ℏ2∂2E∂kα∂kβ\frac1{m^*_{\alpha\beta}} =\frac1{\hbar^2}\frac{\partial^2E}{\partial k_\alpha\partial k_\beta}

其中 mαβ∗m^*_{\alpha\beta} 为有效质量张量分量(kg\mathrm{kg});α,β=x,y,z\alpha,\beta=x,y,z;EE 的单位为 J\mathrm J;kk 的单位为 m−1\mathrm{m^{-1}}。适用条件:能带极值附近、抛物线近似。

主轴方向:选 kx,ky,kzk_x,k_y,k_z 沿张量主轴,则非对角分量为零,

1mα∗=1ℏ2∂2E∂kα2\frac1{m^*_\alpha}=\frac1{\hbar^2}\frac{\partial^2E}{\partial k_\alpha^2}

一维情形:

1m∗=1ℏ2d2Edk2\frac1{m^*}=\frac1{\hbar^2}\frac{\mathrm d^2E}{\mathrm dk^2}

一维紧束缚 ss 带:

E(k)=εs−J0−2J1cos⁡kaE(k)=\varepsilon_s-J_0-2J_1\cos ka

v(k)=1ℏdEdk=2aJ1ℏsin⁡kav(k)=\frac1\hbar\frac{\mathrm dE}{\mathrm dk} =\frac{2aJ_1}{\hbar}\sin ka

m∗=ℏ2d2E/dk2=ℏ22J1a2cos⁡kam^*=\frac{\hbar^2}{\mathrm d^2E/\mathrm dk^2} =\frac{\hbar^2}{2J_1a^2\cos ka}

其中 m∗m^* 为有效质量(kg\mathrm{kg});J1J_1 为最近邻跃迁积分(J\mathrm J);aa 为晶格常数(m\mathrm m);kk 为波矢(m−1\mathrm{m^{-1}})。

有效质量的特点:

  • 带底 k=0k=0:cos⁡ka=1\cos ka=1,m∗=ℏ2/(2J1a2)>0m^*=\hbar^2/(2J_1a^2)>0;
  • 带顶 k=±π/ak=\pm\pi/a:cos⁡ka=−1\cos ka=-1,m∗=−ℏ2/(2J1a2)<0m^*=-\hbar^2/(2J_1a^2)<0;
  • 拐点 cos⁡ka=0\cos ka=0:m∗→∞m^*\to\infty,此时电子平均加速度为零。

有效质量为负的原因:晶体中电子除受外力外,还受晶格作用力 Fl\mathbf F_l:

a=1m(F+Fl)\mathbf a=\frac1m(\mathbf F+\mathbf F_l)

定义有效质量 m∗m^* 使

a=1m∗F\mathbf a=\frac1{m^*}\mathbf F

则 m∗m^* 已概括了晶格的作用。当电子从外场获得的动量大于传递给晶格的动量时 m∗>0m^*>0;小于时 m∗<0m^*<0;全部交给晶格时 m∗→∞m^*\to\infty。

易错点:

  • 有效质量不是电子的真实质量,可为正、负或无穷大。
  • 有效质量是张量,只有在主轴方向或各向同性时才可写成标量。
  • 带底 m∗>0m^*>0,带顶 m∗<0m^*<0;不要把符号记反。
  • ℏdk/dt=F\hbar\mathrm d\mathbf k/\mathrm dt=\mathbf F 中的 F\mathbf F 是外力,不包括晶格力;晶格力已包含在 m∗m^* 中。
  • 准经典近似要求波包尺度远大于原胞,即 Δk≪2π/a\Delta k\ll2\pi/a;不能用来描述原子尺度的运动。

4.9 导体、半导体和绝缘体 ​

能带填充:

  • 满带:被电子占满的允许带;
  • 空带:每个能级都没有电子;
  • 价带:由价电子能级分裂形成的能带,通常是能量最高的能带,可能满、可能不满;
  • 导带:导体中部分价电子存在的不满带;
  • 禁带:相邻能带之间电子不允许具有的能量间隙。

基本观点:

  1. 满带电子不导电。
  2. 不满带中的电子对导电有贡献。

无外场时:E(k)=E(−k)E(\mathbf k)=E(-\mathbf k),v(k)=−v(−k)v(\mathbf k)=-v(-\mathbf k);电子填充 k\mathbf k 与 −k-\mathbf k 的几率相等,总电流为零。

有外场时:

ℏdkdt=F=−eε\hbar\frac{\mathrm d\mathbf k}{\mathrm dt}=\mathbf F=-e\boldsymbol\varepsilon

其中 F\mathbf F 为电场力(N\mathrm N);ee 为电子电荷量(C\mathrm C);ε\boldsymbol\varepsilon 为电场强度(V/m\mathrm{V/m});k\mathbf k 为波矢(m−1\mathrm{m^{-1}})。

  • 满带:k\mathbf k 轴上各点以相同速度移动,不改变均匀填满状态;从 A′A' 移出的电子同时从 AA 移入,仍保持满带,不产生电流。
  • 不满带:外场使电子分布向一方偏移,破坏原来的对称分布,电子电流只部分抵消,产生净电流。

三类固体的能带解释:

类型能带特征导电性
导体有导带(部分填充带),或能带重叠导电
半导体价带满,导带空,禁带较窄(约 2 eV2\ \mathrm{eV})热激发或掺杂后导电
绝缘体价带满,导带空,禁带较宽(33–10 eV10\ \mathrm{eV})几乎不导电

易错点:

  • 「满带不导电」不等于「满带中电子不动」;电子仍在运动,只是总电流为零。
  • 导体不一定「导带半满」;能带重叠也可使导体有部分填充的能带。
  • 半导体与绝缘体的区别在禁带宽度,不在是否有禁带。
  • 外场下 k\mathbf k 在 kk 空间匀速移动;满带电子会从布里渊区一侧移出、另一侧移入,保持满带。

4.10 能态密度 ​

定义:在 E∼E+ΔEE\sim E+\Delta E 能量范围内的能态数目为 ΔZ\Delta Z,则能态密度

N(E)=lim⁡ΔE→0ΔZΔE=dZdEN(E)=\lim_{\Delta E\to0}\frac{\Delta Z}{\Delta E}=\frac{\mathrm dZ}{\mathrm dE}

其中 N(E)N(E) 的单位为 J−1\mathrm{J^{-1}}(或 eV−1\mathrm{eV^{-1}},按能量单位而定);ZZ 为能态数(无量纲);EE 为能量(J\mathrm J)。

重要公式:能态密度的一般表达式

N(E)=2V(2π)3∫SEds∣∇kE(k)∣N(E)=\frac{2V}{(2\pi)^3}\int_{S_E}\frac{\mathrm ds}{\left\lvert\nabla_{\mathbf k}E(\mathbf k)\right\rvert}

其中 N(E)N(E) 的单位为 J−1\mathrm{J^{-1}};VV 为晶体体积(m3\mathrm{m^3});SES_E 为等能面;ds\mathrm ds 为等能面上的面积元(m−2\mathrm{m^{-2}});∇kE\nabla_{\mathbf k}E 为能量梯度(J⋅m\mathrm{J\cdot m});因子 22 来自电子自旋。适用条件:三维晶体、等能面非奇异。

推导思路:kk 空间中状态密度为 V/(2π)3V/(2\pi)^3(计入自旋为 2V/(2π)32V/(2\pi)^3);两等能面间体积元为 ds dk\mathrm ds\,\mathrm dk,且 dE=∣∇kE∣dk\mathrm dE=\lvert\nabla_{\mathbf k}E\rvert\mathrm dk,故

N(E)=2V(2π)3∫SEds∣∇kE∣N(E)=\frac{2V}{(2\pi)^3}\int_{S_E}\frac{\mathrm ds}{\lvert\nabla_{\mathbf k}E\rvert}

自由电子的能态密度:

E(k)=ℏ2k22m,∣∇kE∣=ℏ2kmE(\mathbf k)=\frac{\hbar^2k^2}{2m},\qquad \lvert\nabla_{\mathbf k}E\rvert=\frac{\hbar^2k}{m}

等能面为球面,半径 k=2mE/ℏk=\sqrt{2mE}/\hbar。

重要公式:自由电子三维能态密度

N(E)=V2π2(2mℏ2)3/2EN(E)=\frac{V}{2\pi^2}\left(\frac{2m}{\hbar^2}\right)^{3/2}\sqrt E

其中 N(E)N(E) 的单位为 J−1\mathrm{J^{-1}};VV 为晶体体积(m3\mathrm{m^3});mm 为电子质量(kg\mathrm{kg});EE 为能量(J\mathrm J)。适用条件:自由电子、三维、抛物线色散、E>0E>0。

低维自由电子能态密度:

维度N(E)N(E)特点
一维L(2m)1/24πℏE−1/2\dfrac{L(2m)^{1/2}}{4\pi\hbar}E^{-1/2}随 EE 减小而增大
二维L2mπℏ2\dfrac{L^2m}{\pi\hbar^2}与 EE 无关,为常数
三维L3(2m)3/22π2ℏ3E1/2\dfrac{L^3(2m)^{3/2}}{2\pi^2\hbar^3}E^{1/2}随 E1/2E^{1/2} 增大

表中 LL 为长度(m\mathrm m);mm 为电子质量(kg\mathrm{kg});EE 为能量(J\mathrm J)。待核对:课件给出的二维、三维公式已含自旋因子 22;一维公式按标准推导若含自旋应为 L(2m)1/22πℏE−1/2\dfrac{L(2m)^{1/2}}{2\pi\hbar}E^{-1/2},课件写作 L(2m)1/24πℏE−1/2\dfrac{L(2m)^{1/2}}{4\pi\hbar}E^{-1/2},相当于少了一个自旋因子,使用时请以课件为准并注意这一差异。

近自由电子的能态密度:在第一布里渊区内,等能面基本保持球面,N(E)∝E1/2N(E)\propto E^{1/2};接近布里渊区边界时,等能面外凸;当 EE 超过第二布里渊区最低能量 EBE_B 时,N(E)N(E) 由零迅速增大。能态密度在布里渊区边界附近出现异常,是能带结构的重要特征(补充:这种异常称为 van Hove 奇点)。

与模式密度的比较:

  • 模式密度 ρ(ω)\rho(\omega):单位频率间隔内的振动模式数;
  • 能态密度 N(E)N(E):单位能量间隔内的电子能态数;
  • 两者形式相似,都是「状态密度 × 等能面面积 / 梯度」。

易错点:

  • 能态密度公式中的因子 22 是自旋简并,不要漏。
  • 三维自由电子 N(E)∝E1/2N(E)\propto E^{1/2},二维为常数,一维为 E−1/2E^{-1/2};维数不同,规律不同。
  • N(E)N(E) 的单位取决于能量单位:用 J\mathrm J 时为 J−1\mathrm{J^{-1}},用 eV\mathrm{eV} 时为 eV−1\mathrm{eV^{-1}}。
  • 近自由电子的 N(E)N(E) 不是简单 E1/2E^{1/2};在布里渊区边界附近会偏离。

4.11 恒定磁场中电子的运动 ​

4.11.1 准经典运动 ​

恒定磁场中,准经典运动方程为

v(k)=1ℏ∇kE(k),ℏdkdt=q v(k)×B\mathbf v(\mathbf k)=\frac1\hbar\nabla_{\mathbf k}E(\mathbf k),\qquad \hbar\frac{\mathrm d\mathbf k}{\mathrm dt}=q\,\mathbf v(\mathbf k)\times\mathbf B

其中 v\mathbf v 为电子速度(m/s\mathrm{m/s});k\mathbf k 为波矢(m−1\mathrm{m^{-1}});qq 为电荷量(C\mathrm C,电子 q=−eq=-e);B\mathbf B 为磁感应强度(T\mathrm T);E(k)E(\mathbf k) 为能带能量(J\mathrm J)。

重要结论:

  1. 沿磁场方向的 kk 分量不随时间改变;
  2. 电子的能量 E(k)E(\mathbf k) 不随时间改变;
  3. 电子在 kk 空间的等能面上运动。

4.11.2 自由电子的量子理论 ​

无磁场时:

H^=p^22m=−ℏ22m∇2\hat H=\frac{\hat{\mathbf p}^2}{2m}=-\frac{\hbar^2}{2m}\nabla^2

有磁场时,把 p^\hat{\mathbf p} 换成 p^+qA\hat{\mathbf p}+q\mathbf A:

H^=12m(p^+qA)2\hat H=\frac1{2m}\left(\hat{\mathbf p}+q\mathbf A\right)^2

其中 A\mathbf A 为矢量势(T⋅m\mathrm{T\cdot m}),满足 B=∇×A\mathbf B=\nabla\times\mathbf A。取 B=(0,0,B)\mathbf B=(0,0,B),Landau 规范 A=(−By,0,0)\mathbf A=(-By,0,0)。

波函数取

ψ=ei(kxx+kzz)ϕ(y)\psi=e^{i(k_xx+k_zz)}\phi(y)

代入薛定谔方程,得 ϕ(y)\phi(y) 满足谐振子方程,其中心位置

y0=ℏkxqBy_0=\frac{\hbar k_x}{qB}

振动圆频率

ω0=qBm\omega_0=\frac{qB}{m}

重要公式:朗道能级

En=(n+12)ℏω0+ℏ2kz22m,n=0,1,2,⋯E_n=\left(n+\frac12\right)\hbar\omega_0+\frac{\hbar^2k_z^2}{2m}, \qquad n=0,1,2,\cdots

其中 EnE_n 为第 nn 个朗道能级(J\mathrm J);ω0=qB/m\omega_0=qB/m 为回旋频率(rad/s\mathrm{rad/s});qq 为电荷量(C\mathrm C);BB 为磁感应强度(T\mathrm T);mm 为电子质量(kg\mathrm{kg});kzk_z 为沿磁场方向的波矢(m−1\mathrm{m^{-1}})。适用条件:自由电子、均匀恒定磁场、非相对论。

物理意义:xx-yy 平面内的圆周运动对应简谐运动,能量量子化,形成朗道能级;沿磁场方向的运动仍为自由电子。

4.11.3 晶体中的有效质量近似 ​

晶体中电子在磁场中的哈密顿量可近似写成

H^=12m∗(p^+qA)2\hat H=\frac1{2m^*}\left(\hat{\mathbf p}+q\mathbf A\right)^2

其中 m∗m^* 为有效质量(kg\mathrm{kg})。电子(或空穴)做螺旋运动,回旋频率

ω0=qBm∗\omega_0=\frac{qB}{m^*}

回旋共振:在垂直磁场方向加频率为 ω\omega 的交变电场,当 ω=ω0\omega=\omega_0 时,交变电场能量被共振吸收。从量子观点看,相当于电子在朗道能级之间跃迁。

实验条件:

ω0τ≫1,ℏω0≫kBT\omega_0\tau\gg1,\qquad \hbar\omega_0\gg k_{\mathrm B}T

其中 τ\tau 为平均自由时间(s\mathrm s);TT 为温度(K\mathrm K)。前者要求高纯晶体、液氦温度,以获得足够长的自由时间;后者要求量子化效应显著。

旋转椭球等能面的回旋共振有效质量:

重要公式:回旋共振有效质量

ωc=eBmc∗\omega_c=\frac{eB}{m_c^*}

1mc∗2=cos⁡2θmt∗2+sin⁡2θml∗2\frac1{m_c^{*2}} =\frac{\cos^2\theta}{m_t^{*2}} +\frac{\sin^2\theta}{m_l^{*2}}

其中 ωc\omega_c 为回旋频率(rad/s\mathrm{rad/s});ee 为电子电荷量(C\mathrm C);BB 为磁感应强度(T\mathrm T);mc∗m_c^* 为回旋共振有效质量(kg\mathrm{kg});mt∗,ml∗m_t^*,m_l^* 分别为横向、纵向有效质量(kg\mathrm{kg});θ\theta 为磁场与椭球纵轴的夹角(rad\mathrm{rad})。适用条件:等能面为旋转椭球、单带近似。

实例:

  • Si:导带底在 ⟨100⟩\langle100\rangle 方向,共有 66 个等价点;mt∗=0.19mm_t^*=0.19m,ml∗=0.92mm_l^*=0.92m。
  • Ge:导带底在 ⟨111⟩\langle111\rangle 方向(L 点),共有 88 个等价点;mt∗=0.082mm_t^*=0.082m,ml∗=1.59mm_l^*=1.59m。

易错点:

  • 朗道能级中的 ω0=qB/m\omega_0=qB/m 是自由电子回旋频率;晶体中要换成 m∗m^*。
  • 回旋共振有效质量 mc∗m_c^* 不是简单的 mt∗m_t^* 或 ml∗m_l^*,而是与磁场方向有关。
  • 回旋共振条件 ω0τ≫1\omega_0\tau\gg1 和 ℏω0≫kBT\hbar\omega_0\gg k_{\mathrm B}T 不能漏。
  • 磁场沿 zz 方向时,kzk_z 守恒、能量守恒,电子在等能面上运动,不是在整个 kk 空间自由运动。

4.12 太阳能电池材料专题 ​

半导体与 PN 结:

  • 材料按导电性分为导体、半导体、绝缘体;太阳能电池材料为半导体。
  • PN 结形成:
    1. P 区空穴浓度高、N 区电子浓度高,浓度差引起扩散运动;
    2. 扩散形成空间电荷区(耗尽层);
    3. 空间电荷区产生内建电场,方向由 N 区指向 P 区,阻止多子扩散,称为阻挡层;
    4. 内建电场使少子发生漂移运动,方向与扩散相反;
    5. 扩散电流与漂移电流达到动态平衡时,PN 结总电流为零,空间电荷区宽度稳定。
  • 光生电动势:光照产生电子—空穴对,内建电场将其分离,光生电动势方向为 P 正 N 负;外电路中的光电流方向与内建电场方向相反。

易错点:

  • 内建电场方向是 N→P;光生电动势方向是 P→N。
  • 扩散运动与漂移运动方向相反;平衡时二者电流大小相等、方向相反。
  • 半导体导电能力随温度升高而增强,与金属相反。

4.13 易错点与常考题型 ​

易错点

  1. 布洛赫定理写成 ψ(r+R)=ψ(r)\psi(\mathbf r+\mathbf R)=\psi(\mathbf r);正确形式是乘 eik⋅Re^{i\mathbf k\cdot\mathbf R}。
  2. 把 ℏk\hbar\mathbf k 当作真实动量;它是准动量。
  3. 近自由电子近似在布里渊区边界仍用非简并微扰,导致分母为零。
  4. 禁带宽度写成 ∣Vn∣\lvert V_n\rvert;应为 2∣Vn∣2\lvert V_n\rvert。
  5. 三维带隙条件写成 k⋅Gn=0\mathbf k\cdot\mathbf G_n=0;应为 k⋅Gn=−12Gn2\mathbf k\cdot\mathbf G_n=-\tfrac12G_n^2。
  6. fcc 第一布里渊区说成八面体;应为十四面体。
  7. 紧束缚近似最近邻数或格矢漏项;sc、fcc、bcc 的带宽分别 12J1,16J1,16J112J_1,16J_1,16J_1。
  8. 有效质量在带顶为负,带底为正;符号记反。
  9. 能态密度漏自旋因子 22;三维自由电子 N(E)∝E1/2N(E)\propto E^{1/2}。
  10. 导体、半导体、绝缘体的区别只看「有无禁带」;实际还要看禁带宽度和填充情况。

常考题型

  1. 证明布洛赫定理;由边界条件求 k\mathbf k 的取值。
  2. 用近自由电子近似求一级、二级能量修正;在布里渊区边界用简并微扰求能隙。
  3. 判断能带在布里渊区边界的劈裂方向;解释禁带宽度与周期势场的关系。
  4. 画出扩展/周期/简约三种能带图,并进行折叠。
  5. 证明三维带隙条件与布拉格反射条件等价。
  6. 判断 fcc/bcc 的第一布里渊区形状与体积;标出 Γ,X,M,R\Gamma,X,M,R 点。
  7. 用紧束缚近似求 sc/fcc/bcc 的 ss 能带、带底/带顶/带宽。
  8. 由 E(k)E(k) 求 v(k)v(k)、m∗(k)m^*(k),讨论正负与无穷大。
  9. 解释满带不导电、不满带导电;判断导体/半导体/绝缘体。
  10. 推导自由电子能态密度;比较一维、二维、三维的 N(E)N(E)。
  11. 求朗道能级、回旋频率、回旋共振有效质量。

第 5 章 金属电子论 ​

5.1 经典模型:特鲁德理论 ​

模型要点:

  • 价电子组成电子气体,电子间无相互作用,电子与离子实间除碰撞外无相互作用;
  • 电子服从麦克斯韦—玻尔兹曼统计;
  • 碰撞使电子达到热平衡;碰撞后运动方向随机,初速度对平均速度无贡献;
  • 平均自由时间(弛豫时间)τ\tau 为连续两次碰撞之间的平均时间。

金属性质的经典解释:

  • 导电性:外电压下自由电子定向移动;
  • 导热性:电子碰撞传递能量;
  • 延展性:自由电子只起胶合作用,原子间易滑动而不断裂;
  • 金属光泽:金属吸收广谱光并重新反射,故不透明。

电子平均漂移速度:

vd=−eτmE\mathbf v_{\mathrm d}=-\frac{e\tau}{m}\mathbf E

其中 vd\mathbf v_{\mathrm d} 为漂移速度(m/s\mathrm{m/s});ee 为电子电荷量(C\mathrm C);τ\tau 为平均自由时间(s\mathrm s);mm 为电子质量(kg\mathrm{kg});E\mathbf E 为电场强度(V/m\mathrm{V/m})。

重要公式:特鲁德电导率与电阻率

σ=ne2τm,ρ=1σ=mne2τ\sigma=\frac{ne^2\tau}{m},\qquad \rho=\frac1\sigma=\frac{m}{ne^2\tau}

其中 σ\sigma 为电导率(S/m\mathrm{S/m});ρ\rho 为电阻率(Ω⋅m\mathrm{\Omega\cdot m});nn 为电子数密度(m−3\mathrm{m^{-3}});ee 为电子电荷量(C\mathrm C);τ\tau 为平均自由时间(s\mathrm s);mm 为电子质量(kg\mathrm{kg})。适用条件:经典自由电子气、各向同性、稳态。

平均自由程:

λ=vthτ\lambda=v_{\mathrm{th}}\tau

其中 λ\lambda 为平均自由程(m\mathrm m);vthv_{\mathrm{th}} 为平均热运动速度(m/s\mathrm{m/s});τ\tau 为平均自由时间(s\mathrm s)。典型金属 λ∼10−10\lambda\sim10^{-10}–10−9 m10^{-9}\ \mathrm m,τ∼10−15\tau\sim10^{-15}–10−14 s10^{-14}\ \mathrm s。

经典电子比热:

CVcl=32nkBC_V^{\mathrm{cl}}=\frac32nk_{\mathrm B}

其中 CVC_V 的单位为 J/(K⋅m3)\mathrm{J/(K\cdot m^3)}(单位体积);nn 为电子数密度(m−3\mathrm{m^{-3}});kBk_{\mathrm B} 为玻尔兹曼常量。按经典理论,电子比热应与晶格比热同量级。

成功与失败:

  • 成功:解释了金属直流电导、弛豫时间、平均自由程的量级。
  • 失败:电子比热理论值约 1.5nkB1.5nk_{\mathrm B},与实验相比过大(实验与经典之比约 0.010.01);无法解释电阻率随温度的变化;补充:经典模型也无法解释为什么只有费米面附近的电子参与输运。

易错点:

  • 经典模型把电子当经典粒子,违反泡利不相容原理,因此电子比热严重偏离实验。
  • τ\tau 是平均自由时间,不是平均自由程;λ=vτ\lambda=v\tau。
  • 电导率公式中的 nn 是电子数密度,不是原子数密度。

5.2 费米分布函数 ​

费米—狄拉克分布:

重要公式:费米分布函数

f(E)=1e(E−EF)/kBT+1f(E)=\frac{1}{e^{(E-E_{\mathrm F})/k_{\mathrm B}T}+1}

其中 f(E)f(E) 为能量为 EE 的本征态被电子占据的几率(无量纲);EE 为电子能量(J\mathrm J 或 eV\mathrm{eV});EFE_{\mathrm F} 为费米能级(化学势,J\mathrm J 或 eV\mathrm{eV});kBk_{\mathrm B} 为玻尔兹曼常量;TT 为绝对温度(K\mathrm K)。适用条件:热平衡、费米子、电子气。

物理意义:热平衡下,能量为 EE 的本征态被电子占据的几率,即平均占有数。

T=0 KT=0\ \mathrm K 的极限:

f(E)={1,E<EF12,E=EF0,E>EFf(E)= \begin{cases} 1,&E<E_{\mathrm F}\\[2pt] \dfrac12,&E=E_{\mathrm F}\\[2pt] 0,&E>E_{\mathrm F} \end{cases}

即在 E=EFE=E_{\mathrm F} 处发生陡变:在绝对零度,费米面以内的状态全被占据,费米面以外全空。

T≠0 KT\ne0\ \mathrm K:

  • 费米面附近 ±kBT\pm k_{\mathrm B}T 范围内 f(E)f(E) 快速变化;
  • 温度越高,f(E)f(E) 越偏离 0 K0\ \mathrm K 的阶跃线;
  • E=EFE=E_{\mathrm F} 时 f=1/2f=1/2;
  • E=EF−kBTE=E_{\mathrm F}-k_{\mathrm B}T 时 f≈0.73f\approx0.73,E=EF+kBTE=E_{\mathrm F}+k_{\mathrm B}T 时 f≈0.27f\approx0.27;
  • E=EF−2kBTE=E_{\mathrm F}-2k_{\mathrm B}T 时 f≈0.88f\approx0.88,E=EF+2kBTE=E_{\mathrm F}+2k_{\mathrm B}T 时 f≈0.12f\approx0.12。

电子总数:

N=∑if(Ei)=∫0∞f(E)N(E) dEN=\sum_i f(E_i)=\int_0^\infty f(E)N(E)\,\mathrm dE

其中 NN 为电子总数(无量纲);EiE_i 为本征态能量(J\mathrm J);N(E)N(E) 为能态密度(J−1\mathrm{J^{-1}});f(E)f(E) 为费米分布。适用条件:热平衡、能量连续近似。

数量级:室温 kBT≈2.6×10−2 eVk_{\mathrm B}T\approx2.6\times10^{-2}\ \mathrm{eV};金属费米能 EFE_{\mathrm F} 为几个 eV\mathrm{eV},故 kBT≪EFk_{\mathrm B}T\ll E_{\mathrm F}。

易错点:

  • 费米分布是「占据几率」,不是「电子数」;电子数还要乘能态密度再积分。
  • EFE_{\mathrm F} 是化学势,随温度略有变化;0 K 时的值记为 EF0E_{\mathrm F}^0。
  • 费米面不是「电子云的边界」,而是 kk 空间中占据态与空态的边界。
  • 费米分布与玻色分布、麦克斯韦—玻尔兹曼分布不同,不能混用。

5.3 费米能、费米面与费米温度 ​

5.3.1 0 K 的费米能 ​

自由电子能态密度:

N(E)=CE1/2,C=4πV(2mh2)3/2N(E)=CE^{1/2},\qquad C=4\pi V\left(\frac{2m}{h^2}\right)^{3/2}

其中 N(E)N(E) 的单位为 J−1\mathrm{J^{-1}};VV 为晶体体积(m3\mathrm{m^3});mm 为电子质量(kg\mathrm{kg});hh 为普朗克常量(J⋅s\mathrm{J\cdot s});CC 的单位为 J−3/2\mathrm{J^{-3/2}}。

0 K 时 f(E)=1f(E)=1(E<EF0E<E_{\mathrm F}^0),故

N=∫0EF0N(E) dE=23C(EF0)3/2N=\int_0^{E_{\mathrm F}^0}N(E)\,\mathrm dE =\frac23C\left(E_{\mathrm F}^0\right)^{3/2}

重要公式:0 K 费米能

EF0=ℏ22m(3π2n)2/3E_{\mathrm F}^0=\frac{\hbar^2}{2m}\left(3\pi^2n\right)^{2/3}

其中 EF0E_{\mathrm F}^0 为 0 K 费米能(J\mathrm J 或 eV\mathrm{eV});ℏ\hbar 为约化普朗克常量;mm 为电子质量(kg\mathrm{kg});n=N/Vn=N/V 为电子数密度(m−3\mathrm{m^{-3}})。适用条件:自由电子气、三维、0 K、抛物线色散。

证明:由 N=(2/3)C(EF0)3/2N=(2/3)C(E_{\mathrm F}^0)^{3/2} 与 C=4πV(2m/h2)3/2C=4\pi V(2m/h^2)^{3/2},

(EF0)3/2=3N2C=3n8π(h22m)3/2\left(E_{\mathrm F}^0\right)^{3/2} =\frac{3N}{2C} =\frac{3n}{8\pi}\left(\frac{h^2}{2m}\right)^{3/2}

两边取 2/32/3 次方,并利用 h=2πℏh=2\pi\hbar,得

EF0=ℏ22m(3π2n)2/3E_{\mathrm F}^0=\frac{\hbar^2}{2m}\left(3\pi^2n\right)^{2/3}

费米半径:

kF=(3π2n)1/3k_{\mathrm F}=\left(3\pi^2n\right)^{1/3}

其中 kFk_{\mathrm F} 为费米半径(m−1\mathrm{m^{-1}});nn 为电子数密度(m−3\mathrm{m^{-3}})。

费米能的其他写法:

EF0=ℏ2kF22mE_{\mathrm F}^0=\frac{\hbar^2k_{\mathrm F}^2}{2m}

费米动量与费米速度:

pF=ℏkF,vF=ℏkFmp_{\mathrm F}=\hbar k_{\mathrm F},\qquad v_{\mathrm F}=\frac{\hbar k_{\mathrm F}}{m}

其中 pFp_{\mathrm F} 为费米动量(kg⋅m/s\mathrm{kg\cdot m/s});vFv_{\mathrm F} 为费米速度(m/s\mathrm{m/s})。

费米温度:

TF=EF0kBT_{\mathrm F}=\frac{E_{\mathrm F}^0}{k_{\mathrm B}}

其中 TFT_{\mathrm F} 为费米温度(K\mathrm K);EF0E_{\mathrm F}^0 为 0 K 费米能(J\mathrm J);kBk_{\mathrm B} 为玻尔兹曼常量。典型金属 TF∼104T_{\mathrm F}\sim10^4–105 K10^5\ \mathrm K。

0 K 电子平均能量:

Eˉ=∫0EF0E N(E) dE∫0EF0N(E) dE=25C(EF0)5/223C(EF0)3/2=35EF0\bar E=\frac{\int_0^{E_{\mathrm F}^0}E\,N(E)\,\mathrm dE}{\int_0^{E_{\mathrm F}^0}N(E)\,\mathrm dE} =\frac{\frac25C(E_{\mathrm F}^0)^{5/2}}{\frac23C(E_{\mathrm F}^0)^{3/2}} =\frac35E_{\mathrm F}^0

重要公式:0 K 电子平均能量

Eˉ=35EF0\bar E=\frac35E_{\mathrm F}^0

其中 Eˉ\bar E 为电子平均能量(J\mathrm J);EF0E_{\mathrm F}^0 为 0 K 费米能(J\mathrm J)。适用条件:自由电子气、三维、0 K。

结论:绝对零度下电子仍有很大的平均能量,而经典理论给出零。原因是泡利不相容原理:电子不能全部填在最低能级,必须逐级填充到费米能。

5.3.2 有限温度的费米能 ​

重要公式:有限温度费米能

EF=EF0[1−π212(kBTEF0)2]E_{\mathrm F}=E_{\mathrm F}^0\left[1-\frac{\pi^2}{12}\left(\frac{k_{\mathrm B}T}{E_{\mathrm F}^0}\right)^2\right]

其中 EFE_{\mathrm F} 为温度 TT 下的费米能(J\mathrm J);EF0E_{\mathrm F}^0 为 0 K 费米能(J\mathrm J);TT 为温度(K\mathrm K);kBk_{\mathrm B} 为玻尔兹曼常量。适用条件:kBT≪EF0k_{\mathrm B}T\ll E_{\mathrm F}^0、自由电子气、三维。

结论:温度升高,费米能级略有下降;在通常温度下 kBT/EF0∼10−2k_{\mathrm B}T/E_{\mathrm F}^0\sim10^{-2},故 EF≈EF0E_{\mathrm F}\approx E_{\mathrm F}^0。

有限温度的平均能量:

Eˉ=35EF0[1+5π212(kBTEF0)2]\bar E=\frac35E_{\mathrm F}^0\left[1+\frac{5\pi^2}{12}\left(\frac{k_{\mathrm B}T}{E_{\mathrm F}^0}\right)^2\right]

其中 Eˉ\bar E 为平均能量(J\mathrm J)。适用条件同上。

易错点:

  • EF0E_{\mathrm F}^0 与 EFE_{\mathrm F} 不同:前者是 0 K0\ \mathrm K 值,后者是有限温度值;通常 EFE_{\mathrm F} 略低于 EF0E_{\mathrm F}^0。
  • 费米能不是「最高能量」的严格说法;0 K 时它是最高占据能级,有限温度时是化学势。
  • 费米温度 TFT_{\mathrm F} 不是实际温度,而是 EF/kBE_{\mathrm F}/k_{\mathrm B} 的能量标度。
  • 计算 EF0E_{\mathrm F}^0 时 nn 是电子数密度,不是原子数密度;对单价金属二者相等,对多价金属不等。

5.4 电子热容 ​

电子平均能量:

Eˉ=1N∫0∞Ef(E)N(E) dE\bar E=\frac1N\int_0^\infty E f(E)N(E)\,\mathrm dE

其中 Eˉ\bar E 为电子平均能量(J\mathrm J);NN 为电子总数;f(E)f(E) 为费米分布;N(E)N(E) 为能态密度(J−1\mathrm{J^{-1}})。

电子热容:

重要公式:电子热容

CVe=π22NkBkBTEF0=γTC_V^{\mathrm e}=\frac{\pi^2}{2}Nk_{\mathrm B}\frac{k_{\mathrm B}T}{E_{\mathrm F}^0} =\gamma T

γ=π2NkB22EF0\gamma=\frac{\pi^2Nk_{\mathrm B}^2}{2E_{\mathrm F}^0}

其中 CVeC_V^{\mathrm e} 为电子热容(J/K\mathrm{J/K});NN 为电子总数;kBk_{\mathrm B} 为玻尔兹曼常量;TT 为温度(K\mathrm K);EF0E_{\mathrm F}^0 为 0 K 费米能(J\mathrm J);γ\gamma 为电子热容系数(J/K2\mathrm{J/K^2})。适用条件:kBT≪EF0k_{\mathrm B}T\ll E_{\mathrm F}^0、自由电子气、三维。

推导要点:只有费米面附近约 kBTk_{\mathrm B}T 范围内的电子能被热激发;能态密度在该范围内近似为常数 N(EF)N(E_{\mathrm F}),对 Ef(E)N(E)E f(E)N(E) 作 Sommerfeld 展开,保留 T2T^2 项,即得 Eˉ=35EF0[1+5π212(kBT/EF0)2]\bar E=\frac35E_{\mathrm F}^0[1+\frac{5\pi^2}{12}(k_{\mathrm B}T/E_{\mathrm F}^0)^2],再对 TT 求导得 CVe=γTC_V^{\mathrm e}=\gamma T。

金属总热容:

CV=CVe+CVph=γT+bT3C_V=C_V^{\mathrm e}+C_V^{\mathrm ph}=\gamma T+bT^3

其中 CVe=γTC_V^{\mathrm e}=\gamma T 为电子贡献;CVph=bT3C_V^{\mathrm ph}=bT^3 为晶格(声子)贡献;γ,b\gamma,b 为材料常数。适用条件:低温(T≪θDT\ll\theta_{\mathrm D})。

低温下的主导项:

  • 极低温:γT\gamma T 可能大于 bT3bT^3,电子热容不可忽略;
  • 较高温度:bT3bT^3 增长更快,声子热容占主导;
  • 高温:CV→3nNkBC_V\to3nNk_{\mathrm B},电子贡献相对很小。

为什么电子热容很小:金属中大多数电子的能量远低于费米能,受泡利原理限制不能参与热激发;只有费米面附近约 kBTk_{\mathrm B}T 范围内的电子能被激发,这部分电子数占总数的比例约为 kBT/EF0∼10−2k_{\mathrm B}T/E_{\mathrm F}^0\sim10^{-2},故电子热容比经典值小约两个数量级。

易错点:

  • 电子热容与 TT 成正比,不是常数;经典理论给出的 32nkB\frac32nk_{\mathrm B} 是错误的。
  • 只有费米面附近的电子参与热激发,不是所有电子。
  • 低温金属热容写成 γT+bT3\gamma T+bT^3;不能只写 T3T^3。
  • γ\gamma 与 EF0E_{\mathrm F}^0 成反比:费米能越大,电子热容越小。

5.5 易错点与常考题型 ​

易错点

  1. 经典模型电子比热写成 32nkB\frac32nk_{\mathrm B} 并认为与实验相符;实际远大于实验。
  2. 费米分布中的 EFE_{\mathrm F} 当作 0 K 值;有限温度下 EFE_{\mathrm F} 略低。
  3. 电子总数写成 N=∫N(E)dEN=\int N(E)\mathrm dE,漏掉费米分布 f(E)f(E)。
  4. 费米能公式漏掉 (3π2n)2/3(3\pi^2n)^{2/3} 中的 33 或 2/32/3 次方。
  5. 费米速度写成 vF=kF/mv_{\mathrm F}=k_{\mathrm F}/m,漏 ℏ\hbar。
  6. 电子热容写成与 TT 无关;应为 γT\gamma T。
  7. 低温金属热容只写 T3T^3,漏掉电子项 γT\gamma T。
  8. 认为所有电子都参与热激发;实际只有费米面附近 kBTk_{\mathrm B}T 范围内的电子。

常考题型

  1. 由电子数密度 nn 求 EF0E_{\mathrm F}^0、kFk_{\mathrm F}、vFv_{\mathrm F}、TFT_{\mathrm F}。
  2. 证明 0 K 下 Eˉ=35EF0\bar E=\frac35E_{\mathrm F}^0。
  3. 由 N(E)N(E) 与 f(E)f(E) 求电子总数、平均能量。
  4. 推导或说明电子热容 CVe=γTC_V^{\mathrm e}=\gamma T,解释为什么远小于经典值。
  5. 解释金属低温热容 γT+bT3\gamma T+bT^3 的来源与温度依赖。
  6. 比较经典模型与量子模型对电导、比热的解释。

附录 A 物理量与公式大总结 ​

表中「章节」指本复习资料的章号;公式中的符号含义以正文对应章节为准。带「补充」的条目为课件外补齐的内容。

物理量名称符号单位对应公式适用条件所在章节
正格矢Rl\mathbf R_lm\mathrm mRl=l1a1+l2a2+l3a3\mathbf R_l=l_1\mathbf a_1+l_2\mathbf a_2+l_3\mathbf a_3布拉伐格子、li∈Zl_i\in\mathbb Z1.3
原胞体积Ω\Omegam3\mathrm{m^3}Ω=a1⋅(a2×a3)\Omega=\mathbf a_1\cdot(\mathbf a_2\times\mathbf a_3)三维晶格1.3
晶胞体积VVm3\mathrm{m^3}V=nΩV=n\Omega晶胞含 nn 个格点1.3
致密度ρ\rho无量纲ρ=N⋅43πR3V\rho=\dfrac{N\cdot\frac43\pi R^3}{V}硬球模型、原子相切1.2
配位数CN\mathrm{CN}无量纲最近邻原子数给定结构1.2
晶列指数[l1l2l3][l_1l_2l_3]无量纲R=l1a1+l2a2+l3a3\mathbf R=l_1\mathbf a_1+l_2\mathbf a_2+l_3\mathbf a_3原胞基矢、lil_i 互质1.5
密勒指数(hkl)(hkl)无量纲截距倒数化互质整数晶胞基矢1.5
面间距ddm\mathrm md=2π∣Gh∣d=\dfrac{2\pi}{\lvert\mathbf G_h\rvert}晶面指数按原胞基矢定义1.6
倒格子基矢bi\mathbf b_im−1\mathrm{m^{-1}}b1=2πa2×a3Ω\mathbf b_1=2\pi\dfrac{\mathbf a_2\times\mathbf a_3}{\Omega}三维晶格1.6
倒格矢Kh\mathbf K_hm−1\mathrm{m^{-1}}Kh=h1b1+h2b2+h3b3\mathbf K_h=h_1\mathbf b_1+h_2\mathbf b_2+h_3\mathbf b_3hi∈Zh_i\in\mathbb Z1.6
倒格子原胞体积Ω∗\Omega^*m−3\mathrm{m^{-3}}Ω∗=(2π)3Ω\Omega^*=\dfrac{(2\pi)^3}{\Omega}三维晶格1.6
允许旋转角θ\thetarad\mathrm{rad}2cos⁡θ=m2\cos\theta=m,m=0,±1,±2m=0,\pm1,\pm2晶体平移周期性1.7
结合能WWJ\mathrm JW=E0−ENW=E_0-E_N自由原子结合成晶体2.2
互作用势u(r)u(r)J\mathrm Ju=−Arm+Brnu=-\dfrac{A}{r^m}+\dfrac{B}{r^n}球对称闭壳层、n>mn>m2.2
互作用力f(r)f(r)N\mathrm Nf=−dudrf=-\dfrac{\mathrm du}{\mathrm dr}有心力、准静态2.2
平衡间距r0r_0m\mathrm mdudr∣r0=0\left.\dfrac{\mathrm du}{\mathrm dr}\right\rvert_{r_0}=0稳定平衡2.2
平衡势能u(r0)u(r_0)J\mathrm Ju(r0)=−Ar0m(1−mn)u(r_0)=-\dfrac{A}{r_0^m}\left(1-\dfrac mn\right)n>mn>m2.2
体积弹性模量KKPa\mathrm{Pa}K=Vd2UdV2∣V0K=V\left.\dfrac{\mathrm d^2U}{\mathrm dV^2}\right\rvert_{V_0}小形变、各向同性2.2
压缩系数κ\kappaPa−1\mathrm{Pa^{-1}}κ=−1V(∂V∂P)T=1K\kappa=-\dfrac1V\left(\dfrac{\partial V}{\partial P}\right)_T=\dfrac1K小形变2.2
马德隆常数α\alpha无量纲α=∑j≠i(±1aj)\alpha=\sum_{j\ne i}\left(\pm\dfrac1{a_j}\right)离子晶体、刚性离子2.4
电负性χ\chi无量纲χ=0.18(I+A)\chi=0.18(I+A)Mulliken 经验式,I,AI,A 取 eV\mathrm{eV}(待核对)2.1
弹性系数β\betaN/m\mathrm{N/m}f=−βδf=-\beta\delta,β=(d2Udδ2)a\beta=\left(\dfrac{\mathrm d^2U}{\mathrm d\delta^2}\right)_a简谐近似、平衡位置附近3.1
单原子链频率ω(q)\omega(q)rad/s\mathrm{rad/s}ω=2βm∣sin⁡qa2∣\omega=2\sqrt{\dfrac{\beta}{m}}\left\lvert\sin\dfrac{qa}{2}\right\rvert一维单原子链、最近邻3.1
最大频率ωmax⁡\omega_{\max}rad/s\mathrm{rad/s}ωmax⁡=2βm\omega_{\max}=2\sqrt{\dfrac{\beta}{m}}一维单原子链3.1
声速vsv_sm/s\mathrm{m/s}vs=aβmv_s=a\sqrt{\dfrac{\beta}{m}}长波极限 qa≪1qa\ll13.1
一维波矢qqm−1\mathrm{m^{-1}}q=2πhNaq=\dfrac{2\pi h}{Na}玻恩—卡门边界条件3.1
双原子链频率ω±\omega_\pmrad/s\mathrm{rad/s}ω±2=βM+mMm{1±[1−4Mm(M+m)2sin⁡2(qa)]1/2}\omega_\pm^2=\beta\dfrac{M+m}{Mm}\left\{1\pm\left[1-\dfrac{4Mm}{(M+m)^2}\sin^2(qa)\right]^{1/2}\right\}一维双原子链、最近邻3.2
折合质量μ\mukg\mathrm{kg}1μ=1M+1m\dfrac1\mu=\dfrac1M+\dfrac1m双原子链3.2
三维波矢q\mathbf qm−1\mathrm{m^{-1}}q=∑ihiNibi\mathbf q=\sum_i\dfrac{h_i}{N_i}\mathbf b_i三维晶格、玻恩—卡门条件3.3
波矢密度ρq\rho_qm−3\mathrm{m^{-3}}ρq=Vc(2π)3\rho_q=\dfrac{V_c}{(2\pi)^3}三维晶体3.3
声子能量EphE_{\mathrm{ph}}J\mathrm JEph=ℏωqE_{\mathrm{ph}}=\hbar\omega_{\mathbf q}简谐近似、单声子3.3
声子准动量p\mathbf pkg⋅m/s\mathrm{kg\cdot m/s}p=ℏq\mathbf p=\hbar\mathbf q简谐近似3.3
平均声子数nn无量纲n=1eℏω/kBT−1n=\dfrac{1}{e^{\hbar\omega/k_{\mathrm B}T}-1}热平衡、玻色统计3.3
模式密度ρ(ω)\rho(\omega)s/rad\mathrm{s/rad}ρ=∑αVc(2π)3∫Sωds∣∇qωα∣\rho=\sum_\alpha\dfrac{V_c}{(2\pi)^3}\displaystyle\int_{S_\omega}\dfrac{\mathrm ds}{\lvert\nabla_{\mathbf q}\omega_\alpha\rvert}三维、等频面非奇异3.5
德拜频率ωD\omega_{\mathrm D}rad/s\mathrm{rad/s}ωD=v(6π2na)1/3\omega_{\mathrm D}=v(6\pi^2n_{\mathrm a})^{1/3}(补充)连续介质、三支声学波3.5
德拜温度θD\theta_{\mathrm D}K\mathrm KθD=ℏωDkB\theta_{\mathrm D}=\dfrac{\hbar\omega_{\mathrm D}}{k_{\mathrm B}}德拜模型3.5
爱因斯坦温度θE\theta_{\mathrm E}K\mathrm KθE=ℏωkB\theta_{\mathrm E}=\dfrac{\hbar\omega}{k_{\mathrm B}}爱因斯坦模型3.5
晶格热容CVC_VJ/K\mathrm{J/K}CV=3nNkBf(θT)C_V=3nNk_{\mathrm B}f\left(\dfrac{\theta}{T}\right)爱因斯坦/德拜模型3.5
杜隆—珀替热容CVC_VJ/K\mathrm{J/K}CV=3nNkBC_V=3nNk_{\mathrm B}高温经典极限3.5
布洛赫波函数ψk\psi_{\mathbf k}—ψk=eik⋅ruk(r)\psi_{\mathbf k}=e^{i\mathbf k\cdot\mathbf r}u_{\mathbf k}(\mathbf r)周期势场、单电子4.2
简约波矢k\mathbf km−1\mathrm{m^{-1}}k=∑iliNibi\mathbf k=\sum_i\dfrac{l_i}{N_i}\mathbf b_i玻恩—卡门边界条件4.2
近自由电子能量EkE_kJ\mathrm JEk=Ek0+∑n≠0∣Vn∣2Ek0−Ek−2πn/a0E_k=E_k^0+\sum_{n\ne0}\dfrac{\lvert V_n\rvert^2}{E_k^0-E_{k-2\pi n/a}^0}非简并微扰、远离布里渊区边界4.3
禁带宽度EgE_gJ\mathrm JEg=2∣Vn∣E_g=2\lvert V_n\rvert近自由电子、布里渊区边界4.3
三维带隙条件k\mathbf km−1\mathrm{m^{-1}}2k⋅Gn+Gn2=02\mathbf k\cdot\mathbf G_n+\mathbf G_n^2=0近自由电子、布里渊区边界4.5
布拉格条件——2dsin⁡θ=nλ2d\sin\theta=n\lambda晶面反射4.5
紧束缚能带Ei(k)E_i(\mathbf k)J\mathrm JEi=εi−J0−∑RsJ(Rs)eik⋅RsE_i=\varepsilon_i-J_0-\sum_{\mathbf R_s}J(\mathbf R_s)e^{i\mathbf k\cdot\mathbf R_s}紧束缚、最近邻4.7
能带宽度ΔE\Delta EJ\mathrm Jsc:12J112J_1;fcc、bcc:16J116J_1紧束缚 ss 带4.7
电子速度v\mathbf vm/s\mathrm{m/s}v=1ℏ∇kE(k)\mathbf v=\dfrac1\hbar\nabla_{\mathbf k}E(\mathbf k)准经典近似4.8
有效质量m∗m^*kg\mathrm{kg}1mαβ∗=1ℏ2∂2E∂kα∂kβ\dfrac1{m^*_{\alpha\beta}}=\dfrac1{\hbar^2}\dfrac{\partial^2E}{\partial k_\alpha\partial k_\beta}能带极值附近、抛物线近似4.8
准动量p\mathbf pkg⋅m/s\mathrm{kg\cdot m/s}p=ℏk\mathbf p=\hbar\mathbf k准经典近似4.8
能态密度N(E)N(E)J−1\mathrm{J^{-1}}N(E)=2V(2π)3∫SEds∣∇kE∣N(E)=\dfrac{2V}{(2\pi)^3}\displaystyle\int_{S_E}\dfrac{\mathrm ds}{\lvert\nabla_{\mathbf k}E\rvert}三维、等能面非奇异4.10
自由电子能态密度N(E)N(E)J−1\mathrm{J^{-1}}N(E)=V2π2(2mℏ2)3/2EN(E)=\dfrac{V}{2\pi^2}\left(\dfrac{2m}{\hbar^2}\right)^{3/2}\sqrt E三维自由电子、E>0E>04.10
朗道能级EnE_nJ\mathrm JEn=(n+12)ℏω0+ℏ2kz22mE_n=\left(n+\dfrac12\right)\hbar\omega_0+\dfrac{\hbar^2k_z^2}{2m}自由电子、均匀磁场4.11
回旋频率ωc\omega_crad/s\mathrm{rad/s}ωc=eBmc∗\omega_c=\dfrac{eB}{m_c^*}晶体中电子/空穴、单带4.11
回旋共振质量mc∗m_c^*kg\mathrm{kg}1mc∗2=cos⁡2θmt∗2+sin⁡2θml∗2\dfrac1{m_c^{*2}}=\dfrac{\cos^2\theta}{m_t^{*2}}+\dfrac{\sin^2\theta}{m_l^{*2}}旋转椭球等能面4.11
电导率σ\sigmaS/m\mathrm{S/m}σ=ne2τm\sigma=\dfrac{ne^2\tau}{m}特鲁德模型、稳态5.1
电阻率ρ\rhoΩ⋅m\mathrm{\Omega\cdot m}ρ=mne2τ\rho=\dfrac{m}{ne^2\tau}特鲁德模型5.1
平均自由程λ\lambdam\mathrm mλ=vτ\lambda=v\tau经典自由电子气5.1
费米分布f(E)f(E)无量纲f(E)=1e(E−EF)/kBT+1f(E)=\dfrac{1}{e^{(E-E_{\mathrm F})/k_{\mathrm B}T}+1}热平衡、费米子5.2
0 K 费米能EF0E_{\mathrm F}^0J\mathrm JEF0=ℏ22m(3π2n)2/3E_{\mathrm F}^0=\dfrac{\hbar^2}{2m}(3\pi^2n)^{2/3}三维自由电子、0 K5.3
费米半径kFk_{\mathrm F}m−1\mathrm{m^{-1}}kF=(3π2n)1/3k_{\mathrm F}=(3\pi^2n)^{1/3}三维自由电子5.3
费米速度vFv_{\mathrm F}m/s\mathrm{m/s}vF=ℏkFmv_{\mathrm F}=\dfrac{\hbar k_{\mathrm F}}{m}三维自由电子5.3
费米温度TFT_{\mathrm F}K\mathrm KTF=EF0kBT_{\mathrm F}=\dfrac{E_{\mathrm F}^0}{k_{\mathrm B}}三维自由电子5.3
0 K 平均能量Eˉ\bar EJ\mathrm JEˉ=35EF0\bar E=\dfrac35E_{\mathrm F}^0三维自由电子、0 K5.3
电子热容CVeC_V^{\mathrm e}J/K\mathrm{J/K}CVe=γTC_V^{\mathrm e}=\gamma T低温、kBT≪EF0k_{\mathrm B}T\ll E_{\mathrm F}^05.4
电子热容系数γ\gammaJ/K2\mathrm{J/K^2}γ=π2NkB22EF0\gamma=\dfrac{\pi^2Nk_{\mathrm B}^2}{2E_{\mathrm F}^0}三维自由电子5.4
金属总热容CVC_VJ/K\mathrm{J/K}CV=γT+bT3C_V=\gamma T+bT^3低温金属5.4

附录 B 常用常数与符号 ​

常数/符号数值单位
约化普朗克常量 ℏ\hbar1.0546×10−341.0546\times10^{-34}J⋅s\mathrm{J\cdot s}
普朗克常量 hh6.626×10−346.626\times10^{-34}J⋅s\mathrm{J\cdot s}
玻尔兹曼常量 kBk_{\mathrm B}1.38×10−231.38\times10^{-23}J/K\mathrm{J/K}
电子电荷量 ee1.602×10−191.602\times10^{-19}C\mathrm C
电子质量 mem_e9.11×10−319.11\times10^{-31}kg\mathrm{kg}
真空介电常数 ε0\varepsilon_08.854×10−128.854\times10^{-12}F/m\mathrm{F/m}
电子伏特11eV=1.602×10−19 J\mathrm{eV}=1.602\times10^{-19}\ \mathrm J

补充:ℏ,h,e,me,ε0\hbar,h,e,m_e,\varepsilon_0 的数值为通用物理常数;课件中明确给出的有 kB=1.38×10−23 J/Kk_{\mathrm B}=1.38\times10^{-23}\ \mathrm{J/K}、ℏ=1.0546×10−34 J⋅s\hbar=1.0546\times10^{-34}\ \mathrm{J\cdot s}、m=9×10−31 kgm=9\times10^{-31}\ \mathrm{kg}(费米能数量级估算用)。

符号约定:

  • a,b,ca,b,c:晶格常数/晶胞基矢长度(m\mathrm m);
  • a1,a2,a3\mathbf a_1,\mathbf a_2,\mathbf a_3:原胞基矢(m\mathrm m);
  • b1,b2,b3\mathbf b_1,\mathbf b_2,\mathbf b_3:倒格子基矢(m−1\mathrm{m^{-1}});
  • k,q\mathbf k,\mathbf q:波矢(m−1\mathrm{m^{-1}});ω\omega:角频率(rad/s\mathrm{rad/s});ν=ω/2π\nu=\omega/2\pi:频率(Hz\mathrm{Hz});
  • nn:原胞内原子数(无量纲)或电子数密度(m−3\mathrm{m^{-3}}),按上下文区分;
  • NN:原胞数或电子总数,按上下文区分;
  • VcV_c:晶体体积(m3\mathrm{m^3});Ω\Omega:原胞体积(m3\mathrm{m^3});
  • TT:绝对温度(K\mathrm K);kBk_{\mathrm B}:玻尔兹曼常量(J/K\mathrm{J/K})。