固体物理复习：从Miller指数到色心一条线打通

简介这是四川大学固体物理课程的期末复习大纲PPT讲义面向物理、材料等相关专业本科生及考研备考者用于考前快速梳理晶体结构、衍射、晶格振动等核心考点。内容覆盖晶面晶向标定与Miller指数、七种典型晶体结构、晶体点群对称操作、Bragg定律与Laue方程、倒格子及布里渊区、离子/共价/金属/分子晶体结合方式并延伸至单双原子链振动、声子、固体比热模型和常见缺陷类型。资源为1个pptx文件压缩包约133KB虽体量精简但知识点密度高便于直接用于期末背诵与查漏补缺。目前已有95人浏览学习。通过该讲义可系统掌握sc/bcc/fcc/hcp等结构的倒格子关系理解结构因子与原子散射因子的计算以及Madelung势、Lennard-Jones势等考点是考前快速过一遍固体物理重点的高效资料。1. 固体物理期末复习从 Miller 指数到色心一条线打通四川大学固体物理期末复习大纲这份 PPT 一共 32 页把晶体结构、衍射、倒格子、晶格振动、晶体结合、固体缺陷六个板块压缩成了一条可执行的复习链。做器件物理、材料模拟或者半导体工艺的人回头看这门课最有价值的不是背下某个结构而是理解正空间和倒空间两套描述语言如何通过衍射条件彼此切换。这篇博文把大纲里每个知识点的推导逻辑、关键公式的物理边界、以及考试和面试中常见的丢分点逐一拆开配合可复现的 Python 计算脚本和参数表格让复习可以直接落在纸上而不是停在课件翻页上。2. 晶体结构和 Miller 指数先把实空间的格点排布算清楚2.1 Miller 指数的确定方法与高频易错点Miller 指数hkl标定的标准流程是取晶面在三根基矢上的截距取倒数再约化成互质整数。这里有一个高频丢分点——截距为负时要在对应指数上方加横线例如截距为 a、-b、2c 时倒数为 1、-1、1/2约化后得到(2̄12)——注意负号写在数字上方写反整个晶面就反了。另一个容易混淆的地方是截距为零的情形。截距为零意味着晶面平行于该基矢方向对应倒数为无穷大所以指数中会出现 0。例如 (100) 面是垂直于 a₁ 且平行于 a₂、a₃ 的平面而不是经过原点的平面经过原点的平面需要先平移再标定。很多学生在标定 fcc 的 (111) 面时出错原因是没有先画出晶胞内的等效原子位置而是直接用立方晶胞的角点算截距。# 用原子坐标判定某个(hkl)面实际切割了哪些原子 import numpy as np # fcc 晶胞内的四个等效原子位置分数坐标 fcc_atoms [(0,0,0), (0.5,0.5,0), (0.5,0,0.5), (0,0.5,0.5)] def plane_hits(h, k, l, atoms, tol1e-6): hits [] for idx, (x, y, z) in enumerate(atoms): # 平面方程 h*x k*y l*z normalized val h*x k*y l*z if abs(val - round(val)) tol: hits.append(idx) return hits print(fcc(111) 经过的原子:, plane_hits(1,1,1,fcc_atoms)) print(fcc(100) 经过的原子:, plane_hits(1,0,0,fcc_atoms))这段脚本的思路是任意晶面 h x k y l z NN 为整数在晶胞内切割的是那些分数坐标代入后恰好凑成整数的格点。fcc(111) 会经过全部四个原子而 fcc(100) 只经过 (0,0,0) 和 (0,0.5,0.5) 两个原子。原因是 fcc 的面心原子并不都在 (100) 族平面上这个判断直接决定了你在画晶面投影图时会不会漏原子。2.2 五种典型结构的原子坐标、配位数与填充率大纲中列出 sc、bcc、fcc、hcp、diamond、zinc sulfide 六种结构复习时建议用一张表横向对比而不是逐个去背。表中的晶胞原子数、配位数、填充率是按几何关系推出来的推一遍比记三遍牢靠。结构晶胞原子数配位数最近邻距离填充率说明sc16aπ/6 ≈ 0.524简单立方Po 是代表bcc28√3a/2√3π/8 ≈ 0.680体心立方Fe 室温相fcc412√2a/2√2π/6 ≈ 0.740面心立方密堆积hcp2基元12a理想 c/a√2π/6 ≈ 0.740六方密堆Mg、Zndiamond8基元为 2 个 C4√3a/4√3π/16 ≈ 0.340共价键低配位zinc sulfide闪锌矿型4Zn 4S4√3a/4约 0.340与 diamond 同骨架填充率计算以 fcc 为例面对角线上原子相切最近邻距离 √2a/2因此原子半径 r √2a/4四个原子的体积 4 × (4/3)πr³除以晶胞体积 a³得到 √2π/6。diamond 结构的填充率只有 0.34不到 fcc 的一半这个数值直接解释了为什么钻石的密度低于同族元素按密堆积估算的预期值。注意hcp 是布拉菲格子里的简单六方格子加两原子基元不是独立的布拉菲点阵类型。diamond 的结构也不属于金刚石布拉菲点阵——它的基元是两个全同原子空间群是 Fd3̄m。考卷上问布拉菲格子种类有多少时答案只有 14 种sc、fcc、bcc、简单六方都在其中hcp 和 diamond 不在。2.3 对称操作与晶体学限制定理平移对称操作 T(Rn) u₁a₁ u₂a₂ u₃a₃ 是晶体区别于非晶的本质旋转对称操作则要求转角只能是 2π/nn 1、2、3、4、6。为什么不能有 5 次轴因为晶体的平移对称性要求旋转操作必须保持格点间距不变。设想把一个格点绕 5 次轴旋转产生的 5 个等距邻居无法同时满足与相邻格点距离等于原格矢长度的条件只有 n 1、2、3、4、6 能在平面内密铺。这个限制定理在准晶发现前被认为是绝对成立的而准晶的 5 次衍射花样之所以被当成异常长达十年根源就在这里。3. 倒格子、布里渊区与衍射条件正空间看不懂的问题换到倒空间就清楚了3.1 Bragg 定律与 Laue 方程的等价性大纲同时列出 Bragg 定律 2d sinθ nλ 和 Laue 方程 Δk nGh这两个条件本质上是同一个衍射条件的正空间与倒空间表示。Bragg 条件用人眼可见的面间距 d 说话Laue 方程用入射和散射波矢差 Δk 说话。两个条件的等价性可以由倒格矢的几何性质直接推出满足 Δk G 时Δk 的模长等于 2|k|sinθ而 |k| 2π/λ|G| 2π/d代入即得 2d sinθ λ。这里 n 的物理意义是衍射级次 n 对应的是倒格矢 G nG₀G₀ 是垂直于衍射面的最短倒格矢。3.2 倒格子基矢定义与关键性质推导倒格子基矢的定义是 b₁ 2π/V (a₂ × a₃)b₂ 和 b₃ 按轮换顺序取。直接写出的性质有三条需要能够独立推导而不是背诵第一aᵢ·bⱼ 2πδᵢⱼ第二倒格矢 G_hkl hb₁ kb₂ lb₃ 垂直于正空间的 (hkl) 晶面族第三晶面间距 d_hkl 2π/|G_hkl|。垂直关系证明的关键步骤是构造 (hkl) 面上的两条非共线矢量(a₁/h - a₂/k) 和 (a₂/k - a₃/l)让它们分别与 G_hkl 做点积利用 aᵢ·bⱼ 2πδᵢⱼ 得到结果为零即 G_hkl 同时垂直于这两条矢量因而垂直于整个平面。面间距公式则是实数空间周期函数的倒格矢展开的直接推论。3.3 sc、bcc、fcc 互为正倒格子正空间倒空间倒空间晶格常数sc晶格常数 asc晶格常数 2π/a2π/abcc晶格常数 afcc晶格常数 4π/a4π/afcc晶格常数 abcc晶格常数 4π/a4π/abcc 的倒格子是 fcc 这一结论可以由基矢定义直接验证。把 bcc 的初基基矢 a₁ a/2(ŷ ẑ - x̂)、a₂ a/2(ẑ x̂ - ŷ)、a₃ a/2(x̂ ŷ - ẑ) 代入倒基矢定义得到的 b₁ 2π/a(ŷ ẑ)这正是 fcc 的正空间基矢结构。fcc 的倒格子是 bcc 同理可证。这一对关系决定了 Al、Cu 这类 fcc 金属的 X 射线衍射花样中出现的衍射峰对应的是 bcc 倒空间的格点序列。3.4 结构因子与系统消光规律基元内各原子的散射波之间会发生干涉结构因子 S 描述这个干涉效果。给定 (hkl) 反射位置矢量 rⱼ 处的原子散射因子为 fⱼS(hkl) Σⱼ fⱼ exp[i 2π(hxⱼ kyⱼ lzⱼ)]总和遍历基元内全部原子。bcc 结构的基元包含 (0,0,0) 和 (1/2,1/2,1/2) 两个原子代入后 S ∝ 1 exp[iπ(hkl)]。hkl 为奇数时 S 0对应系统消光(100)、(111) 反射在 bcc 中消失。fcc 结构的基元包含四个原子代入后 S ∝ 1 exp[iπ(hk)] exp[iπ(hl)] exp[iπ(kl)]h、k、l 奇偶混合时消光(100)、(110) 反射不会出现。写一段 Python 脚本验证消光规律直接模拟出 fcc 铝粉的粉末衍射峰位序列import numpy as np def structure_factor(h, k, l, basis): S 00j for (x, y, z) in basis: S np.exp(2j*np.pi*(h*x k*y l*z)) return abs(S) fcc_basis [(0,0,0), (0.5,0.5,0), (0.5,0,0.5), (0,0.5,0.5)] hkl_list [(1,0,0),(1,1,0),(1,1,1),(2,0,0),(2,2,0),(3,1,1)] for hkl in hkl_list: print(hkl, |S| , round(structure_factor(*hkl, fcc_basis), 4))运行结果是 (100) 和 (110) 的 |S| 为零(111)、(200)、(220)、(311) 非零。这解释了为什么实际测 fcc 铝的衍射谱时第一强峰是 (111) 而不是 (100)——不是因为你没对准晶体取向而是这个反射在结构上被消光了这是初学衍射的人最容易误判的一个点。4. 晶格振动与热学性质色散关系、声学支/光学支和德拜模型4.1 单原子链色散关系的推导与长波极限单原子链是最简单的振动模型大纲里的色散关系是 ω(q) 2(β/m)^(1/2)|sin(qa/2)|其中 β 是最近邻力常数m 是原子质量a 是晶格常数。推导过程为写出第 n 个原子的运动方程 m·d²uₙ/dt² β(uₙ₊₁ uₙ₋₁ - 2uₙ)代入行波解 uₙ Aexp[i(qna - ωt)]约掉公共因子后得到 ω² (4β/m)sin²(qa/2)开方取正根。长波极限 q → 0 时 sin(qa/2) ≈ qa/2色散关系退化为 ω √(β/m)·|qa|即声速 c √(β/m)·a。这个线性关系意味着长波极限下晶格可以看作连续弹性介质弹性波速由力常数和原子质量决定。注意这里角频率对 q 的依赖是线性的不是二次的这一点在判断声子态密度低频行为时至关重要。4.2 双原子链光学支与声学支的判定双原子链每元胞含质量为 M重原子和 m轻原子的两个原子二元耦合方程的解给出两支色散ω±² β(Mm)/(Mm) ± β√[(Mm)²/(Mm)² - 4sin²(qa)/(Mm)]下标减号对应声学支q → 0 时 ω ≈ 0色散线性重原子和轻原子同向振动代表晶格的刚性平移。下标加号对应光学支q → 0 时 ω √[2β(Mm)/(Mm)]两类原子反向振动质心不动。如果是离子晶体反向振动产生净电偶极矩可与红外光耦合——这正是光学支名称的由来。在布里渊区边界 q π/a 处两支频率分别为 ω₊ √(2β/m)轻原子振动和 ω₋ √(2β/M)重原子振动。判定一个模式是光学支还是声学支考试时看 q→0 时频率是否为零即可。更严格的说法是看原子相对运动方向同相为声学支反相为光学支。4.3 Born-Karman 边界条件与模式计数有限长度 N 个原子的链需要边界条件确定允许的 q 值。周期性边界条件要求位移 uₙ₊ₙ uₙ代入行波解得到 exp(iqNa) 1于是 q 2πl/(Na)l 取整数。在第一布里渊区 -π/a q ≤ π/a 内l -N/2 1 到 N/2共 N 个允许波矢。这个 N 恰好等于元胞数因此每个振动支贡献 N 个模式。推广到三维若晶体有 N 个元胞每个元胞含 m 个原子总振动模式数 N·m 个自由度这是晶格比热计算的基础。许可的波矢数永远是 N正比于实空间体积振动支数为 m·维数其中声学支占维数支光学支占 m·维数 - 维数 支。4.4 态密度与德拜比热的 T³ 定律模式对热容的贡献需要用态密度 g(ω) 加权积分。一比热 C_V ∫ g(ω) ħω · ∂n(ω,T)/∂T dω其中 n 1/[exp(ħω/kBT) - 1] 是玻色-爱因斯坦分布。德拜模型的关键近似是用线性色散 ω vq 和一个截止频率 ω_D 替代真实色散关系G(ω) 9N ω²/ω_D³。低温极限 kBT ≪ ħω_D 时积分上限可以延拓到无穷得到 C_V (12π⁴/5) NkB (T/Θ_D)³。德拜温度 Θ_D ħω_D/kB 是材料的内禀参数铜的 Θ_D 约 343 K钻石约 2230 K。这个 T³ 依赖关系在液氦温度到几十开尔文区间内与实验高度一致是检验德拜模型是否适用的经验判据。4.5 N 过程与 U 过程倒逆散射决定热阻声子-声子散射中动量守恒有两种形式正常过程N 过程满足 q₁ q₂ q₃q₃ 落在第一布里渊区内散射后声子总动量守恒不直接带来热阻倒逆过程U 过程满足 q₁ q₂ q₃ G其中 G 是倒格矢散射前后准动量相差一个倒格矢动量的流向发生了反转是晶格热阻的主要来源。用数值脚本直观展示 U 过程如何发生# 一维链上演示倒逆过程 # 设 q10.6*pi/a, q20.7*pi/a超出第一布里渊区边界 pi/a q_total 0.6 0.7 # 单位为 pi/a if q_total 1.0: G_shift 2.0 # 最近倒格矢单位为 pi/a q3 q_total - G_shift print(fN过程失败发生U过程: q3 {q3} * pi/a) else: print(fN过程q3 {q_total} * pi/a)这段脚本模拟了两个高温声子碰撞后合成波矢越过第一布里渊区边界被倒格矢折回的情形。q₃ 为负值声子运动方向反转净动量传递被破坏这部分动量直接转化为热阻。温度越高激发的高频声子越多U 过程占比越大——这就是晶格热导率随温度升高而下降的本质原因。4.6 频率计算脚本把色散曲线画出来import numpy as np def phonon_dispersion_1d(beta, m, MNone, a1.0, N500): 计算单/双原子链色散关系。M为None时为单原子链。 q np.linspace(0, np.pi/a, N) if M is None: omega 2*np.sqrt(beta/m)*np.abs(np.sin(q*a/2)) return q, omega mu 1/m 1/M omega_sq beta*mu omega_sq beta*np.sqrt(mu**2 - 4*np.sin(q*a/2)**2/(m*M)) omega_minus np.sqrt(beta*mu - beta*np.sqrt(mu**2 - 4*np.sin(q*a/2)**2/(m*M))) return q, omega_sq, omega_minus q, omega phonon_dispersion_1d(beta1.0, m1.0) print(布里渊区边界频率:, omega[-1], ω_max 2√(β/m))脚本输出单原子链在布里渊区边界 q π/a 处的频率为 2√(β/m)。物理含义是相邻原子反相振动恢复力最大因此频率最高q 0 处所有原子同相平动频率为零。提示高频声子的波长已经小到与晶格常数相当连续介质近似不再成立色散曲线从线性弯成平顶。这也是德拜模型在高温端必然失效的原因——它把 ω ∝ q 的线性关系用到了超出适用范围的地方。5. 晶体的结合Madelung 势、Lennard-Jones 势与结合能判据5.1 五类结合方式的特征对照类型键的物理来源典型晶体特征熔点/硬度离子键库仑吸引 泡利排斥NaCl、CsCl绝缘、透明、高熔点、低膨胀高共价键电子对共享方向性金刚石、Si、Ge低配位、高硬度、弱导电很高金属键价电子离域正离子实浸没在电子海中Cu、Fe、Na密堆积、高延展、金属光泽中分子键范德瓦尔斯瞬时偶极-诱导偶极惰性气体晶体密堆积、低熔点、透明、易压缩低氢键电负性原子与 H 间强偶极作用冰、DNA 碱基对弱相互作用、方向性低石墨是典型的混合键晶体层内 C-C 以 sp² 共价键连接层间是范德瓦尔斯力垂直层方向的电子离域又带来类似金属键的性质因此它导电、润滑、耐高温三种特性集于一身。5.2 双粒子势 U(r) 的通用分析工具大纲中给出一般形式的相互作用势 u(r) -α/rⁿ β/rᵐα 0, β 0以及平衡条件 ∂U/∂r|_{r₀} 0。处理这类问题有三步固定的计算流程。第一步由平衡条件求 r₀nα/r₀ⁿ⁺¹ mβ/r₀ᵐ⁺¹解出 r₀ (mβ/nα)^(1/(m-n))。第二步求结合能 U(r₀) 并化简为只含 α 和 r₀ 的形式。第三步由稳定条件判断哪个指数更大——U(r) 在 r₀ 处取极小值要求 U(r₀) 0代入第一、二步结果后得到 m n。m n 的物理含义是排斥项必须比吸引项更短程、更陡峭。排斥项起源于电子云重叠时的泡利排斥随距离下降极快吸引项库仑力或范德瓦尔斯力是长程力衰减慢。如果 m n势能曲线没有极小值原子会被压碎或飞散晶体无法结合。5.3 离子晶体的 Madelung 势计算离子晶体的每对离子能量写为 U(r) N/2[-αe²/(4πε₀r) B/rⁿ]其中 α 是 Madelung 常数定义为对晶格中所有正负离子对的库仑势求和α Σ±1/pᵢⱼpᵢⱼ 是以最近邻距离为单位的距离求和排除 i j。直接逐项求和收敛极慢——三维离子晶体中长程库仑项按 1/r 衰减条件收敛需要对离子分组后再求和。NaCl 结构的 α 1.7476CsCl 结构 α 1.7627一维无限链的 α 2ln2 ≈ 1.386。考试中常见的计算题方式是给出最近邻距离 r₀ 和 Madelung 常数求结合能或压缩率。先把吸引项积分求和得到 -αe²/4πε₀r₀再由平衡条件确定排斥项参数 B两步即可。5.4 Lennard-Jones 势的参数标定与平衡距离惰性气体晶体用 LJ 势描述U(r) 4ε[(σ/r)¹² - (σ/r)⁶]其中 ε 是势阱深度σ 是 U(r) 0 时的距离r 为原子间距。求平衡距离dU/dr 0得到 (σ/r₀)⁶ 1/2即 r₀ 2^(1/6)σ。势阱深度 U(r₀) -ε。由于 vdW 键极弱ε 通常只有几十 meV惰性气体晶体的熔点极低Ar 的熔点在 84 K 附近。5.5 脚本NaCl Madelung 常数的直接求和收敛演示import numpy as np def madelung_1d(N500): 一维交替离子链的 Madelung 常数 s 0.0 for i in range(1, N1): s (-1)**(i1) / i return 2 * s def madelung_nacl_naive(L): 三维 NaCl 结构直接求和收敛慢 s 0.0 count 0 for i in range(-L, L1): for j in range(-L, L1): for k in range(-L, L1): if i j k 0: continue r np.sqrt(i*i j*j k*k) # NaCl: 奇数坐标之和为阴离子偶数为阳离子 s (-1)**(ijk) / r count 1 return s print(1D α , madelung_1d(), 解析值:, 2*np.log(2)) print(3D 直接求和 L10:, madelung_nacl_naive(10))一维链的求和收敛很快三维直接求和的截断误差却大到无法接受L 10 时结果与 1.7476 仍偏差明显。这是 Ewald 求和方法存在的意义——它把实空间慢收敛的 1/r 求和转化为实空间快速衰减的误差函数项加上倒空间快速衰减的高斯项两部分都指数收敛。复习时理解到直接求和为什么不行、Ewald 为什么能行这个层面才算真正掌握了 Madelung 常数的计算逻辑。6. 固体缺陷、色心与大纲复习闭环6.1 用一组判断条件区分缺陷类型大纲第 19 页集中出现 Schottky 缺陷、Frenkel 缺陷、F 心、FA 心、VK 心、R 心、M 心。记忆这些概念可以用一个统一的判断框架缺陷由什么带走电荷空位、填隙、杂质还是电子以及缺陷之间的空间关系是什么。Schottky 缺陷正负离子空位成对出现晶体体积膨胀密度下降。核心特征是只产生空位不产生填隙。Frenkel 缺陷一个离子离开正常格位变成填隙离子同时留下空位。核心特征是空位和填隙浓度相等密度几乎不变。F 心负离子空位俘获一个电子。碱卤晶体受辐照后变色就是 F 心吸收了可见光。F 心本身是色心的总称其他色心都是在 F 心基础上的变体M 心是两个相邻 F 心R 心是三个相邻 F 心VK 心是相邻卤素离子空位对。色心的能级位于禁带中吸收边红移效应可以用来说明碱卤晶体颜色的变化。扩散机制部分大纲提到空位机制、间隙机制和复合机制。判断用哪种机制看原子尺寸置换型杂质原子走空位机制小尺寸填隙原子如碳在铁中走间隙机制复合机制则适用于间隙原子先变为替换原子再扩散的复合过程。Fick 第一定律 J -D∇c 只适用于稳态扩散非稳态问题必须用第二定律 ∂c/∂t D∇²c。考试最常见的错误是在稳态问题里用第二定律或者在扩散系数随浓度变化明显的体系里强行用常数 D。6.2 面缺陷小角晶界与堆垛层错小角晶界可以看作一列平行刃型位错相邻位错间距 D b/θb 为伯格斯矢量模长θ 为两侧晶粒的取向差。θ 从 1° 变到 5° 时位错间距从约 57b 缩小到约 11b界面能从低角度晶界逐渐过渡到大角度晶界的特征。堆垛层错则发生在密堆积结构的堆垛顺序被打断时。fcc 的理想堆垛是 ABCABCABC如果在某处变成 ABABAB就插入了一段 hcp 堆垛。层错能很高的材料如铝不容易出现层错层错能低的材料如黄铜容易出现这会直接改变材料的加工硬化行为。6.3 把 32 页大纲变成可自测的复习闭环这份大纲结构清晰适合用写-推-默三遍法闭环。第一遍合上 PPT在一张白纸上画出知识树晶体结构、衍射、倒格子、振动、结合、缺陷六个根节点每个根节点下写公式名和它解决的问题。第二遍打开大纲对着每个公式独立推导Miller 指数标定三步骤、倒格子基矢定义式代入 bcc/fcc 验证互倒关系、单原子链色散推导、LJ 势求平衡距离、U(r) 平衡条件推 m n每个推导控制在 10 分钟以内。第三遍用下面这张检查表逐一打勾打勾过程中如果某个条目需要翻课件才能确认就标记为薄弱点回到对应章节重新看。知识点验收标准Miller 指数能标定负截距晶面并写出互质指数不混淆 (100) 与 (ī00)结构填充率能推导 fcc 和 diamond 的填充率说出 hcp 为何不属于布拉菲格子倒格子与布区能证明 bcc 与 fcc 互为正倒格子会用 Wigner-Seitz 元胞构造第一布里渊区衍射消光能判断 bcc 的(100)(111) 和 fcc 的(100)(110) 消光原因色散关系能推导单原子链的 ω~q 关系说出光学支与声学支在 q→0 的差异比热模型能推导德拜 T³ 定律并从 Θ_D 区分金刚石与铜结合能能用平衡条件解 r₀ 并证明 m n 的稳定性来源缺陷能现场举出 F 心、Schottky、Frenkel 各一个实例并说明电荷携带方式固体物理面试的常见追问方式是为什么这里用正空间解释不了。答题的核心能力是把问题翻译成倒空间语言衍射问题用 Δk G振动问题用 q 的布里渊区周期场中的电子问题用能带。这份大纲的 32 页正是围绕这一主线组织的把每个公式自己推一遍胜过把课件翻三遍。本文还有配套的精品资源点击获取