Learning objectives / 学习目标
完成本课后，你应能：区分五类成键的主要电子机制和方向性；从 \(U'(r_0)=0\)、\(U''(r_0)>0\) 判断稳定间距；把势阱曲率连接到一维弹性模量的尺度；用平移有序、晶粒、玻璃化时间尺度与弱相互作用分类固体；在所有结论中声明势模型、温压条件与各向同性近似。
Prerequisites / 先修知识
- Coulomb potential、Pauli exclusion、orbital overlap 与能量最低原理。
- 导数、Taylor expansion、small strain \(\epsilon=\Delta L/L\)。
- Identical Particles：交换对称与 Pauli 不是经典排斥力。
Scope and conventions / 范围与约定
- 研究对象为中性体材料；默认 SI，原子尺度常辅助用 eV、nm，但结果明确换算。
- \(U(r)\) 是固定其余自由度后的一条有效 pair-coordinate 能量曲线；多体金属键、角向共价键不能普遍还原为 pair potential。
- 应力取拉伸为正；小应变、低温、无塑性时才使用线性弹性。
- “材料类别”描述主导机制，不排除混合成键；soft solid 是力学/时间尺度类别，不是第五种晶体学周期性。
Notation and modeling contract / 符号与建模契约
| Symbol | Meaning | SI unit |
|---|---|---|
| \(r\) | representative separation | m |
| \(r_0\) | stable equilibrium separation | m |
| \(U(r)\) | effective interaction energy | J |
| \(E_{\mathrm{coh}}\) | cohesive energy magnitude per atom/formula unit | J |
| \(k_{\mathrm{eff}}=U''(r_0)\) | local bond stiffness | N m\(^{-1}\) |
| \(\epsilon\) | engineering strain | 1 |
| \(\sigma\) | normal stress | Pa |
| \(Y\) | Young modulus in stated direction | Pa |
1. Ionic, covalent, metallic, molecular, and hydrogen bonding / 离子、共价、金属、分子与氢键
Knowledge ID: P09-U01-L01-K01
判断键型时应先问电子自由度如何重新组织，而不是先背“硬、脆、导电”等宏观标签。离子键的库仑项具有长程性，单个离子对的能量不足以决定晶体稳定，必须结合整个配位网络和电中性；共价键的轨道重叠带有角向选择，因此键长与键角共同进入能量；金属键中占据态延伸跨越许多晶胞，局部配位变化会改变能带而非只改变一根“键”。分子晶体与氢键体系又同时受分子取向、极化率和热运动影响。实际材料常为混合机制，例如带部分共价性的离子晶体，因此分类的正确用法是指出主导能标、方向性和电子局域程度，再由实验检验，而不是从一个性质反推唯一键型。
成键不是“两个原子之间多出一根实体棒”，而是 nuclei、electrons 与量子统计共同决定的总能量降低。Ionic bonding 由电荷转移后的长程 Coulomb attraction 主导；covalent bonding 由定向 orbital overlap 与共享电子密度主导；metallic bonding 由离域价电子对离子实的整体稳定作用主导；molecular solids 主要靠 dispersion、dipole interactions；hydrogen bond 是带强方向性的特殊非共价相互作用，常写作 donor–H···acceptor。
图 1 只画局部离子对；真实 ionic crystal 的能量必须对整个 lattice 求和，并有短程排斥稳定间距。
图 2 强调 covalent bond 的方向性；bond angle 改变也会付出能量，因此单一 \(U(r)\) 不足以描述网络固体。
图 3 的电子点不是固定在某两个离子之间；离域性解释金属键缺少固定 bond direction，并为导电性留下可能，但真正输运需能带与散射理论。
图 4 比较较弱且较各向平均的 dispersion 与有 donor/acceptor 几何要求的 hydrogen bond。强弱排序不是普适表：总结合还取决于配位数与结构。
2. Cohesive energy and equilibrium spacing / 结合能与平衡间距
Knowledge ID: P09-U01-L01-K02
实验上还应把结合能与晶格常数放在同一温压条件下比较；若样品含表面、缺陷或不同相，测得平均值可能偏离理想体相。特别要避免把势能 minimum 的深度误作单次断键阈值，真实断裂会选择低能路径并产生新表面。
结合能是一张能量账本，参考态必须与拆解过程一致：把离子晶体拆成中性原子、气态离子或分子会得到不同数值。平衡间距则是总能量的一阶导数为零，并非吸引力“消失”；在该点吸引与排斥贡献恰好平衡。若外压存在，应最小化焓或适当热力学势，而不是仍把真空中的 \(U(r)\) minimum 当作唯一答案。有限温度时自由能还包含振动熵，平均间距可偏离零温势能 minimum。由此可见，势阱深度、minimum 位置和曲率分别回答拆解成本、几何尺度和局部刚度，是三个独立信息；任何经验势都应同时用这三类观测量校准，并检查远距离极限与短距离发散是否合理。
取无限分离原子能量为零，稳定固体的最低能量为负；cohesive energy magnitude 定义为拆成相隔无穷远中性组分所需的正能量。一个教学势为
\[U(r)=-\frac{A}{r^m}+\frac{B}{r^n},\qquad n>m>0.\]
平衡条件与稳定条件分别为
\[U'(r_0)=0,\qquad U''(r_0)>0,\]
故 \(r_0=[nB/(mA)]^{1/(n-m)}\)。若只有 attraction，体系会塌缩；若只有 repulsion，不会形成 bound solid。
图 5 同时编码零力点、势阱深度与两侧不对称；热振动在非谐势中平均位置会偏移，为热膨胀埋下伏笔。
3. Elastic response from interatomic potentials / 由原子间势理解弹性响应
Knowledge ID: P09-U01-L01-K03
从势曲率预测模量时，还需说明加载是定温还是绝热、横向是否自由以及晶体取向。若计算得到负曲率或能量随某个剪切方向下降，说明假定结构发生机械不稳定，而不是得到“负的普通弹性模量”。
从原子势到宏观弹性必须经过“位移场如何分配”的桥梁。若所有键随外加应变仿射伸长，单位体积能量的二阶变化可以直接给 elastic constants；若基元内原子还能内部松弛，实际刚度会比冻结内部坐标的结果小。三维晶体的应变和应力都是二阶张量，立方晶体也通常需要多个独立常数，单一 \(Y\) 只描述指定加载方向和横向边界条件。势的三阶、四阶项控制非线性、声子散射和热膨胀；缺陷滑移则会产生不可逆塑性，已超出局部 minimum 的弹簧图景。计算后应检查能量密度是否为正、对称性是否允许相应耦合，以及小应变极限能否回到 Hooke law。
令 \(x=r-r_0\)，在小位移下
\[U(r)\simeq U(r_0)+\frac12U''(r_0)x^2, \qquad F\simeq-k_{\rm eff}x.\]
对间距 \(a\)、截面积 \(A_c\) 的一维链，单键力 \(F=k_{\mathrm{eff}}a\epsilon\)，故 \(\sigma=F/A_c\) 给出尺度关系 \(Y\sim k_{\mathrm{eff}}a/A_c\)。三维实际模量要由单位体积能量对 strain tensor 二阶求导；晶体通常各向异性。
图 6 比较相同平衡点、不同曲率的势阱；深势阱不必等于陡曲率，故 cohesive energy 与 modulus 不能互换。
图 7 是小应变 chain model；超过 harmonic region 后会出现非线性、缺陷滑移或断裂，不能继续外推 Hooke law。
4. Crystalline, polycrystalline, amorphous, and soft solids / 单晶、多晶、非晶与软固体
Knowledge ID: P09-U01-L01-K04
结构分类要区分空间相关长度和动力学时间尺度。单晶的平移相位与取向能跨越宏观样品，多晶只在每个 grain 内保持这种关系，grain boundary 处原子配位和局部应力不同；非晶体虽无无限周期，却保留由键长、键角规定的短程甚至中程关联。软固体的“固态”常来自网络连通、缠结、界面或玻璃化冻结，其模量比原子晶体低很多，且随频率变化。判断样品时应同时使用 diffraction 的 sharp/diffuse features、显微取向图和流变时间扫描；只凭外观看起来坚硬或透明无法区分。温度、制备速率、老化和观察时间都可能改变类别表现，因此分类应附测量条件。
Single crystal 具有跨样品的取向一致 long-range translational order；polycrystal 由不同取向 grains 拼接，grain boundaries 打断平移相位；amorphous solid 有近邻 short-range order 但无长程周期；soft solids（polymer networks、gels、foams 等）由较弱能标、复杂内部自由度和长 relaxation time 支撑形状，其 response 显著依赖时间、频率和温度。
图 8 用位置关联示意三类结构；非晶并非“完全随机气体”，第一、第二近邻壳层仍可清楚存在。
图 9 表明 soft solid 的分类依赖 \(t_{\rm obs}/\tau\)；在短时间像固体、长时间可能流动。
图 10 是本课建模闭环：成键标签、势参数和结构类别必须相互校验，不能由单一宏观性质反推唯一微观键型。
Worked Examples / 完整例题
Worked Example 1 — 幂律势的平衡间距
已知与任务： 给定势参数，求正的稳定平衡距离并判断 stationary point 类型。
研究对象与系统边界： 只取两个有效粒子的 separation coordinate；外界和其余邻居冻结，能量零点在无限分离。
参考系与正方向： 质心系中 \(r\) 从一个粒子指向另一个，取增大间距为正；最终只报告标量距离。
模型与理由： 采用题给 attraction–repulsion 幂律，因为任务是训练 minimum test；它不是完整晶格势。
方程与计算： 先求 \(U'=0\)，再代入 \(U''\)；保留 SI 参数到最后，避免把 eV 与 J 混用。
单位与有效数字： \(B/A\) 的量纲为 \(\mathrm m^8\)，八次方根为 m；输入三位有效数字，结果取三位。
合理性、极限与解释： \(B\to0\) 时 \(r_0\to0\) 显示失去短程稳定；正曲率确认 minimum，所得亚纳米尺度符合原子间距。
已知： \(U=-A/r+B/r^9\)，且 \(A=2.00\times10^{-28}\ \mathrm{J,m}\)、\(B=1.00\times10^{-106}\ \mathrm{J,m^9}\)。对象/边界： isolated effective pair；忽略多体。坐标： \(r>0\)。模型： stationary minimum。
\[U'=\frac{A}{r^2}-\frac{9B}{r^{10}}=0, \qquad r_0=\left(\frac{9B}{A}\right)^{1/8}=2.15\times10^{-10}\ \mathrm m.\]
\(U''(r_0)=8A/r_0^3>0\)，故为稳定点。量纲为 m；原子尺度合理。结果不代表真实材料的完整 Madelung sum。
Worked Example 2 — 结合能与温度尺度
已知与任务： 将每原子的 \(3.20\ \mathrm{eV}\) cohesive energy 换算为每摩尔能量及形式温度尺度。
研究对象与系统边界： 一摩尔相同原子从指定固体参考态拆为无限远中性基态原子；不追踪实际升温或相变路径。
参考系与正方向： 能量账本不依赖惯性系；把外界对系统输入的拆解功取正。
模型与理由： 使用 \(N_A\) 个等价原子且忽略表面项，适合宏观 bulk 极限。
方程与计算： 依次乘 elementary charge conversion 与 Avogadro constant，再除 \(k_B\) 得单粒子等效温标。
单位与有效数字： eV 转 J 后乘 mol\(^{-1}\) 得 J mol\(^{-1}\)；保留三位有效数字。
合理性、极限与解释： 数值远高于典型熔化焓，因为 melting 不要求完全 atomization；温度尺度不是预测熔点。
已知： 每原子 cohesive energy \(3.20\ \mathrm{eV}\)。对象： 一摩尔原子固体。模型： 独立拆解能量账本，不含相变路径。
\[E_{\rm mol}=3.20(1.602\times10^{-19})(6.022\times10^{23}) =3.09\times10^5\ \mathrm{J,mol^{-1}}.\]
等效 \(E/k_B=3.71\times10^4\ \mathrm K\)。这不是熔点：熔化涉及 entropy、缺陷与两相自由能，而非把所有键独立切断。
Worked Example 3 — 势曲率估算一维模量
已知与任务： 由单键曲率、平衡间距和分配截面积估算轴向 Young modulus。
研究对象与系统边界： 无限、均匀的一维链及分配给它的横截面；忽略横向收缩、缺陷和基元内部松弛。
参考系与正方向： 实验室系中链沿 \(+x\)；拉伸应变与拉应力为正。
模型与理由： 小应变使每键延长 \(a\epsilon\)，并用 harmonic force \(F=k_{\mathrm{eff}}a\epsilon\)；这是 affine-chain approximation。
方程与计算： 用 \(\sigma=F/A_c\) 与 \(Y=\sigma/\epsilon\) 消去 \(F,\epsilon\)，再统一 nm 到 m。
单位与有效数字： \((\mathrm{N,m^{-1}})\mathrm m/\mathrm{m^2}=\mathrm{Pa}\)；数据三位，答案三位。
合理性、极限与解释： \(k_{\mathrm{eff}}\to0\) 时模量趋零；大数值反映单链很硬，不能未经体积平均直接当作多晶 bulk 值。
已知： \(k_{\rm eff}=120\ \mathrm{N,m^{-1}}\)、\(a=0.250\ \mathrm{nm}\)、每链截面积 \(A_c=0.0600\ \mathrm{nm^2}\)。对象： 轴向 chain；小应变。正方向： 拉伸为正。
\[Y=\frac{k_{\rm eff}a}{A_c} =\frac{120(2.50\times10^{-10})}{6.00\times10^{-20}} =5.00\times10^{11}\ \mathrm{Pa}.\]
单位 \(\mathrm{N,m^{-2}}\) 正确；数量级对应 stiff covalent chain，但真实体材料需三维连接与方向平均。
Worked Example 4 — 观察时间决定软物质响应
已知与任务： 对单一 Maxwell relaxation time，计算两种观察时间的 Deborah number 并判断固态或流动态响应。
研究对象与系统边界： 线性粘弹单元在小振幅测试下；温度固定、无老化、无多个 relaxation modes。
参考系与正方向： 样品固定在流变仪坐标中；本题只比较时间尺度，不涉及方向性应力。
模型与理由： \(De=\tau/t_{\mathrm{obs}}\) 直接比较内部松弛与外部实验时钟，适合 regime classification。
方程与计算： 分别代入 \(20.0/0.20\) 与 \(20.0/200\)，再与 1 比较。
单位与有效数字： 秒相除得到无量纲；输入至少两位有效数字，报告三位便于比较。
合理性、极限与解释： \(t_{\mathrm{obs}}\to0\) 时 \(De\to\infty\)、结构来不及松弛；长时间极限 \(De\to0\)，材料有充分时间流动。
已知： Maxwell relaxation time \(\tau=20.0\ \mathrm s\)；测试时间分别 \(0.20\ \mathrm s\)、\(200\ \mathrm s\)。对象： 线性 viscoelastic element；小扰动。定义 Deborah number \(De=\tau/t_{\rm obs}\)。
短测量 \(De=100\)，以 elastic-like 为主；长测量 \(De=0.100\)，显著 relaxation/flow。\(De\) 无量纲；结论只适用于单一 relaxation-time 近似。
Lesson synthesis / 全课综合
电子结构决定 attraction 的来源与方向性，Pauli/核排斥阻止塌缩；二者形成势阱，minimum 给间距、depth 给拆解能、curvature 给小振动和弹性尺度。许多原子如何排列以及观察多久，又决定单晶、多晶、非晶或软固体的宏观表现。可执行流程是：定 constituents → 声明主导键型 → 写总能量 → 求 minimum/curvature → 检查结构有序与时间尺度 → 再预测性质。
Formula-and-decision table / 公式与决策表
| Task | Relation | Use only when | Model boundary or convention |
|---|---|---|---|
| stable spacing | \(U'(r_0)=0\), \(U''(r_0)>0\) | chosen coordinate captures relaxation | stationary point alone may be maximum |
| cohesive energy | \(E_{\rm coh}=U(\infty)-U(r_0)\) | fragments/reference state declared | differs from melting enthalpy |
| harmonic bond | \(k_{\rm eff}=U''(r_0)\) | small \(r-r_0\) | ignores anharmonicity |
| one-dimensional modulus scale | \(Y\sim k_{\rm eff}a/A_c\) | affine chain and assigned area | not full elastic tensor |
| soft-solid regime | \(De=\tau/t_{\rm obs}\) | characteristic \(\tau\) meaningful | real spectra may have many \(\tau\) |
Model boundaries / 模型边界
- Pair potentials are controlled approximations; metallic and covalent energies are generally many-body/orientation dependent.
- Cohesive-energy reference must specify atoms, ions, molecules and electronic states.
- Linear elasticity excludes plasticity, fracture, defects moving between minima and large rotations.
- “Crystalline” concerns long-range order, not perfection; crystals contain thermal motion and defects.
- “Amorphous” does not mean no correlations; “soft” does not mean liquid or structurally disordered.
Misconception clinic / 误区诊所
- “键就是一根弹簧。” 弹簧只是 minimum 附近的二阶近似；键源自总量子能量。
- “cohesive energy 大就一定 modulus 大。” 势阱深度与底部曲率是不同导数信息。
- “金属电子完全自由。” 它们离域但仍受周期势、Pauli filling 与 scattering。
- “非晶没有结构。” 它缺 long-range periodicity，却有 short-range order。
- “polycrystal 等于许多孤立 single crystals。” Grain boundaries 改变输运、强度与散射。
Practice and retrieval / 练习与提取
- [K01] 分别给 ionic 与 covalent bonding 的主导电子图景。
- [K01] 为什么 hydrogen bond 不能简单等同于 dispersion？
- [K02] 对 \(U=-A/r^2+B/r^6\) 求 \(r_0\)。
- [K02] 若 minimum 为 \(-4.0\ \mathrm{eV}\)，cohesive-energy magnitude 是多少？
- [K03] \(U''(r_0)\) 的 SI 单位是什么？
- [K03] 同一 \(a,A_c\) 下曲率加倍，chain-model \(Y\) 如何变？
- [K04] single crystal 与 polycrystal 的取向区别是什么？
- [K04] 为什么 glass 仍会出现宽的近邻峰？
- [K04] \(\tau=5\ \mathrm s\)、\(t_{\rm obs}=0.01\ \mathrm s\)，求 \(De\) 并判断响应。
- [综合] 指出“熔点等于 cohesive energy 除以 \(k_B\)”的两处问题。
Answer 1
Ionic：显著 charge transfer 后的 lattice Coulomb energy；covalent：方向性 orbital overlap 与共享密度降低能量。
Answer 2
Hydrogen bond 有 donor–H···acceptor 的电性、极化与方向几何；dispersion 普遍存在且较非定向。
Answer 3
\(U'=2A/r^3-6B/r^7=0\)，所以 \(r_0=(3B/A)^{1/4}\)，并须检验 \(U''>0\)。
Answer 4
\(4.0\ \mathrm{eV}\)；magnitude 为正，最低能量相对分离态为 \(-4.0\ \mathrm{eV}\)。
Answer 5
\(\mathrm{J,m^{-2}}=\mathrm{N,m^{-1}}\)。
Answer 6
在该近似中加倍。
Answer 7
前者在样品尺度保持同一晶格取向；后者由多个不同取向 grains 构成。
Answer 8
局部配位与键长仍受势能约束，故有 short-range order。
Answer 9
\(De=500\)，观察比 relaxation 快，呈 elastic-like。
Answer 10
熔化比较两相 Gibbs free energy 并含 entropy；熔化不需把材料拆成无限远孤立原子。
Exit ticket / 离堂检验
- 画一条有效势并同时标出 \(r_0\)、\(E_{\rm coh}\)、\(U''(r_0)\) 的物理意义。
- 给一个材料写三层判据：主导键型、空间有序、\(t_{\rm obs}/\tau\)；说明为何三层不能互相替代。
References and source notes / 参考资料与来源说明
以下直接来源于 2026-08-30 实际打开核验；正文重组，TikZ 均由方程原创。
- OpenStax, University Physics Vol. 3, 9.3 Bonding in Crystalline Solids. https://openstax.org/books/university-physics-volume-3/pages/9-3-bonding-in-crystalline-solids
- OpenStax, Chemistry, 10.5 The Solid State of Matter. https://openstax.org/books/chemistry/pages/10-5-the-solid-state-of-matter
- MIT OCW 3.012, Fundamentals of Materials Science: Bonding and crystal structures. https://ocw.mit.edu/courses/3-012-fundamentals-of-materials-science-fall-2005/resources/lec16b/
- MIT OCW 1.978, Problem Set 2, Problem 4: Calculation of Elastic Coefficients from Energy Density Function. https://ocw.mit.edu/courses/1-978-from-nano-to-macro-introduction-to-atomistic-modeling-techniques-january-iap-2007/181227c33a5e1029ceb33659685af9e3_ps2.pdf
- IUCr, An introduction to the scope, potential and applications of X-ray analysis. https://www.iucr.org/education/pamphlets/2/full-text
Source-to-content map
| Source | Used for | Knowledge ID / Section | Model boundary or convention |
|---|---|---|---|
| OpenStax University Physics 9.3 | bonding categories, equilibrium spacing | K01–K02 | introductory effective-potential treatment |
| OpenStax Chemistry 10.5 | ionic/molecular/metallic/network solids, amorphous distinction | K01, K04 | chemistry classification, mixed bonding allowed |
| MIT 3.012 | lattice/bonding and common material structures | K01, K04 | undergraduate materials convention |
| MIT 1.978 Problem Set 2, Problem 4 | energy curvature and elastic coefficients | K03 | continuum coefficient derived from stated atomistic model |
| IUCr pamphlet | crystalline order and lattice context | K04 | crystallographic terminology |
复习与拓展
- 下一课把“周期排列”精确拆成 lattice + basis，并引入 symmetry 与 Miller planes。
- Unit 9.2 将 \(U''(r_0)\) 推广为 coupled dynamical matrix，得到 phonons。
- Unit 9.3 将离域电子放入周期势，说明为何键型与导电分类并非一一对应。