晶体里的缺陷，什么时候是在帮助整体释放应力？

曲面上出现缺陷，究竟解释了什么

同样的胶体颗粒在平面上偏好整齐的六配位排列；在逐渐收腰的液体界面上形成晶体时，晶格中先出现相邻的五配位与七配位，再出现由这些位错排成的线。缺陷增多，但大尺度晶格未必因此更难维持：它们可能正在承担原本分布在整片晶体中的几何冲突。

Irvine、Vitelli与Chaikin在2010年的原始实验中，用共聚焦显微镜重建油—甘油界面上的胶体晶体。颗粒相互排斥，在平面上偏好六配位；液体界面的形状由毛细作用和固定边界控制，颗粒层不能任意把基底展平。拉开两端时，近圆柱界面逐渐变成负高斯曲率的腰部，出现无缺陷、孤立位错、中性位错线，再到带净荷的疤痕及孤立七配位的序列。1

这篇论文称中性位错线为 pleat。它在晶格内部由位错排列造成取向调整，不是薄片向空间鼓起的皱纹。两者都能改变应力，却使用不同自由度：前篇改变曲面形状，本篇在基本固定的曲面上改变晶格连接。

要判断缺陷是否帮助释放应力，至少要分别回答：几何和拓扑强制多少净荷；允许另外增加多少中性缺陷；这些缺陷的位置与取向能否降低能量；实验中它们是否确由这一过程产生。单看一张五、七配位图，不能走完这条因果链。

拓扑数的是净荷，边界也要入账

令$z_i$为粒子的邻居数。在三角晶格内部，定义无量纲配位荷$Q_i=6-z_i$：五配位为$+1$，七配位为$-1$；相应位向错的角荷为$q_i=(\pi/3)Q_i$。相邻的一五一七组成位错，净荷为零，但它仍有非零Burgers矢量，表示绕其一圈后晶格平移的闭合差。12

对没有边界的三角网，$3F=2E$、$\sum_i z_i=2E$，与Euler恒等式$V-E+F=\chi$合用，得到

闭球的$\chi=2$，所以净配位荷为12。如果所有缺陷都是五配位，就需要12个；存在七配位时，五配位也要相应增加。定理约束的是两者之差，不是总缺陷数只能为12。

有边界时，三角形的边计数变成$3F=2E-E_b$。对每条简单边界环，边界边数等于边界顶点数，因而

边界粒子的基准是4，而不是把所有粒子都照内部的6来数。这是原2010补充材料§4给出的开放边界修正。2

本库构造并逐条核了三组网：4

边翻转让原来的两个六配位变为五配位，另两个六配位变为七配位；于是全球从12个五配位变成14个五配位加2个七配位，净荷仍为12。这里只证明计数允许增加中性缺陷，没有算这一边翻转的实际能量收益。

开放六角片的六个角各承担一个边界单位荷，内部可以全是六配位。因此，“圆盘拓扑需要总荷6”不能改说成“圆盘内部必须有六个五配位”。在开放晶帽上，内部与边界怎样分配还要由几何、边界条件和能量决定。圆柱腰部保持环带拓扑，$\chi=0$，即使曲率改变，也不要求非零的内部加边界总荷。

配位荷、实际角亏和弹性应力不是一个量

把六配位周围的六个三角形都当成等边三角形，角度恰好凑成$2\pi$；有五个则少$\pi/3$。但实际晶体的三角形会被拉伸、压缩，角度不再相同。

在分片平面的三角网中，内部顶点的实际角亏是

边界的转角则是$\pi-\sum_f\theta_{if}$。角亏由实际角度决定，配位荷由连接决定，二者仅在相应等角条件下逐点相等。

本库把上面的平面六角片保持所有连接不变，将顶点抬到单位球面的一个帽状区域，再按直边三角面连接。两个网内部配位荷都为0，边界配位荷都为6；但分片平面角亏的内部总和由约0变为0.817464弧度，边界转角由$2\pi$变为5.465722弧度。两者之和仍为$2\pi$。这不是光滑球帽的精确积分，也没有宣称变形后仍无应变；它直接展示了同一连接计数容得下不同的实际几何。

光滑曲面的对应关系是Gauss–Bonnet式

分段光滑边界还要加入角点转角。闭球没有边界项，积分为$4\pi$；开放晶帽的曲率积分则可以连续改变，由边界项补上。不能仅凭曲率积分变化，就把一个整数的内部缺陷数硬改成同样连续的数。

原弹性理论用位向错密度$s(x)=\sum_\alpha q_\alpha\delta(x,x_\alpha)$与高斯曲率$K(x)$的差作为不相容源。在其连续弹性近似中，Airy应力函数满足

能量可写成由这一源、双调和Green函数和边界条件决定的二次型。这里$Y$是二维面内模量，单位为力/长度；不是前篇以Pa计的三维Young模量。2

加入一对净荷为零的位向错并不改变拓扑账，却改变$s(x)$的空间分布。若它抵消了曲率源中一部分代价很大的空间分量，分布应变能就可降低；但缺陷自身的弹性场、相互作用和核心仍然要付费。局部核心附近的应力可能变大，同时整片的弹性能下降。“释放应力”在这里应落实为给定条件下的总弹性能比较，不能等同于每一点都更平静。

走完一个原始模型：为何取向和位置决定收益

2010补充材料§3给出了近圆柱腰部的一对位错模型。设高度为$H$，颈部环向曲率半径为$R_\phi$、纵向曲率半径为$R_z>0$，采用

其高斯曲率为负；本库直接积分得到这一曲面近似的精确关系

后面的弹性计算进一步要求$H/R_z\ll1$，用自由轴向牵引的近圆柱解。无缺陷的曲率应力可以从

得到。取七配位朝向颈部的位错，其Burgers矢量大小为$b$，轴向Peach–Koehler力为

在颈部附近，合适取向的位错被推向外；越过$H/\sqrt{12}$后受到向内的力。只计几何相互作用时，稳定位置因此在$z=H/\sqrt{12}$。把五、七取向对调，几何力和几何势都变号。2

原PDF式16在$\chi_A^G$系数中多印了一个$b$；它与无缺陷源方程、式22的位向错势及式17的力不一致，且量纲不同。本篇用由这些关系一致确定的上式，再独立求四阶导数和应力核对；不把直接照抄那一项称为原公式复现。这一更正及原件定位留在来源记录。

再引入一对Burgers矢量相反的位错，分别位于$\pm z$。它们每一个都由一五一七组成，所以这一“位错对”共含两五两七，净位向错荷仍为0。以无缺陷曲面晶体为参照，原式23给出

即

$a$是晶格间距，$E_c$是每个位错的核心能。三项分别是对数相互作用代价、与背景曲率应力耦合的收益、两份核心代价。前两项不是各自独立的一对正负位向错的全部微观能量；这是原文已约化的两个位错模型。

固定$H=R_\phi=1$、$a=b=0.001$、$E_c=0.05Yb^2$，用$YH^2$作能量单位。这些是透明示例输入，没有拟合实验。只改变

在原文几何位置$z=H/\sqrt{12}$处，实际算得：4

弱曲率时，释放的应变能不足以补偿新增位错；曲率增强后，同样两个位错变得有利。这个差别不来自拓扑要求变化，也不来自把位错“免费”加入。

本库没有把几何项极小位置直接当成全对能量极小。对完整$\Delta E$求导，正取向驻点满足

在$\Gamma=0.05$时，全对能量最低的分离位置为$z/H=0.278605$，比几何项单独给出的0.288675稍靠近颈部；$\Delta E=-1.97265\times10^{-7}YH^2$。在相同位置反转取向，几何项成为正代价；该反向对在本例整个允许间距内都没有低于无缺陷态。

原无缺陷弱圆柱解的能量为

在这个例子中$F_0=1.09083\times10^{-5}YH^2$，加入最优位错对后仍为正，降低约1.81%。负的$\Delta E$表示相对释放，不表示绝对能量变负、全部应力消失，或单对位错已经是所有可能缺陷排列中的全局最优。

尺寸、核心代价和形成时间怎样改变答案

把位错放在几何项极小位置，令$\Delta E=0$，可由原式23直接得到

本库直接用有量纲的式23计算，不依赖原式25另印的$a'$定义。再允许完整能量选择位置，固定$E_c/(Yb^2)=0.05$、$b=a$时，得到：

这里固定$R_\phi/H=1$。较大的$H/a$让释放分布应变的收益更能补偿局部缺陷代价；提高核心能会抬高阈值。本库另实际比较了0.02和0.2两档核心能。它们是原弱圆柱近似内的位错对能量判据，不是完整实验相图。24

两个边界必须留住。其一，对数弹性不适用于重叠的位错核心。计算只在$2z\ge a$作区间比较；弱曲率的最小值常落在这一截止边缘，不能把它解释成准确的束缚对尺寸。上述过零阈值发生在$z/a>270$，远离核心，但模型仍未处理真实核心能、所有边界像力和精确曲面相互作用。

其二，把$E_c$设为零时，截止位置的对数项恰好为零，几何项对任意正曲率已经为负。首次数值求根因此不存在有限起始阈值；本库留快照后把这一情况保留为核心近似失效的诊断，没有捏造一个有限根，也不据此声称真实晶体随任意小曲率成核。

即使存在能量更低的分离位错态，系统也可能跨不过创建核心、交换邻居或从边界进入的能垒。本例的间距扫描有局部峰，但它没有模拟微观位错产生，不能用这一个峰计算实验等待时间。热涨落、黏滞阻力、拉伸速度和已有缺陷都可能改变实际出现时刻。

更多位错能排成pleat，调节晶格取向来分担曲率。原文用无缺陷角度应变能约$Y\theta^2H^2$，与稀疏低角晶界代价约$-YbH\theta\ln\theta$比较，得到$\theta_c\sim(a/H)\ln(H/a)$的量级估计。它解释尺寸越大越容易以位错线释放应变的趋势，省略的几何和核心因素不能拿来给每个实验精确报阈值；单对复算也没有完整求解多条pleat的数量、长度和竞争。

怎样区分释放曲率应力、增加晶格行数和驱动流动

原2010实验最有用的证据不只是缺陷变多，还包括它们的五、七取向和出现位置。七配位朝向负曲率较强的颈部，符合几何力的符号；随负曲率增大，极化位错增多并排成中性线，与背景应变收益变强的预测一致。12

但拉开毛细桥也增加高度，并非只改曲率。原补充视频说明§2明确指出：同时出现另一类未按曲率极化的位错对，它们滑移并增加竖直晶格行数，以适应拉伸；释放曲率应力的轴向五七排列不承担这种增行作用。本文核的是原说明和理论，未重新跟踪视频颗粒。把同一张图中所有位错都归于曲率屏蔽，会漏掉这条实际存在的替代机制。

有判别力的比较因而是：在明确边界、密度和模量下测位错的取向、位置及应变场；分别控制曲率与拉伸行数需求；改变核心代价或相对尺寸，看分离位错态的有利范围是否改变；再改变等待时间或拉伸速度，辨认动力学推迟。仅有净荷恒等式、缺陷与曲率相关，或一个末态图，都不足以独自决定这些机制的权重。

2026年的旋转颗粒研究提供另一个边界对照。其对象是平面圆形容器中的振动驱动颗粒，没有2010实验的曲面高斯曲率。边界附近的疤痕影响摩擦力矩传递；调节顺、逆时针颗粒比例后，边缘与内部流动分离，可出现自剪切、径向密度变化和再入熔化。整体面积分数保持不变，并不代表局部内部密度保持不变。3

原补充材料又给出限制：在实验面积分数0.75和0.79时未观察到边缘与内部的反向流；0.72下较小的500、1000粒子系统仍有疤痕和角速度突降，却没有反向自剪切。更关键的是，原作者的最小模拟只在暂态复现反向流，长期被边缘拖成同向，实验则观察到稳态。缺陷排列相近，不保证同一流动状态，更不能把非平衡稳态用2010的被动自由能极小化解释。

本篇取得两篇主文及补充材料，核相关原理论、实验条件和上述模拟限制；2026数据存档仅取得文件目录元数据，未下载Origin数据文件或重做模拟。被动位错对的实际复算，与原实验、主动流动模型分别记账。

目前能回答什么，哪些还留着

较稳的答案是：拓扑和边界约束净配位荷；曲率产生的分布应变允许中性位错通过合适取向与位置降低能量；新增缺陷是否有利，要用释放量补偿相互作用和核心代价。净荷相同的两片晶体，可以有不同缺陷总数、不同应力和不同可达状态。

本库的独立计算给出三条可回查的分界。第一，精确网格计数和角亏计算显示，边界配位荷、内部配位荷和实际曲率必须分别记录。第二，原位错对模型中，同一拓扑下改变曲率、取向、尺寸和核心代价就能改变能量收益；全对优化位置与几何项单独极小也不同。第三，出现更低能态只回答是否有利，未回答能否在实验时间内形成。

本文还没有得到真实胶体的核心能、完整曲面Green函数下所有缺陷排列的最低能、有限温度缺陷分布或成核时间。原主动颗粒数据未再分析，最小模拟与实验稳态的差异也未解决。若允许基底弯曲、改变颗粒相互作用、让缺陷从边界流入流出，释放应力的资源与比较条件都要重新列明。

来源与复核

<a id="source-s1"></a> [1] Irvine, Vitelli & Chaikin (2010). Pleats in crystals on curved surfaces, Nature 468, 947–950。已取得作者站发表PDF，核缺陷定义、图1—5、晶体序列与Methods。毛细基底固定、pleat中性、曲率及边界条件分别使用；未重做原颗粒跟踪。

<a id="source-s2"></a> [2] 同研究原补充材料。 17页SI，已取得并核§2基底几何、§3式1—32和注8—9、§4开放边界Euler式、补充视频说明。式16的额外$b$由原PDF页面目读及相邻方程核；本篇直接执行式23，不沿用有量纲不明的单独截止符号。

<a id="source-s3"></a> [3] Tiwari et al. (2026). Reentrant melting of scarred odd crystals by self-shear, Nature Communications 17, 1802，2026-01-16。主文及24页SI已取得，核实验/模拟方法和SI§3、5—9、13的密度、尺寸与暂态限制。数据存档仅核目录元数据，未取得逐点Origin数据或运行作者模拟。

<a id="source-s4"></a> [4] 本库实际计算 (2026)。 输入与结果、实际网格和能量曲线输入。research/calculations中的curved_crystal_stress_relief.py、_audit.py、_figures.py实际执行；网格边/面/边界环独立计数、Airy多项式求导、力积分、密扫加标量极小和阈值交叉核通过。自建网格示例与原理论重算不称实验、原作者代码或主动流动复现。