先说清，究竟是哪一种大小停了
把一批能结合的元件放在一起，最直观的预期是不断长大。但若元件各自偏好的连接角度不能在整个结构中同时满足，继续增长会积累应变；这可能让有限组装体胜过更大的结构。
这种机制确实能在明确模型中成立，却有一个容易漏掉的细节：一个方向受限，不意味着所有方向都停。 Hackney、Amey与Grason的2023年原始格点研究观察到有限宽度、可弯曲并支化的虫状组装体；它们的长度和粒子数有分布，并随浓度变化。原研究中的“自限”首先是宽度选择。1
此外，一团东西看起来不再长，可能是自由能选择了有限尺寸，也可能是附近原料耗尽、输运太慢或重排过不了能垒。即使计算出一个紧凑团块的最佳半径，若允许它展成长条、引入缺陷或与别的团块连通，答案还可能改变。
本文沿着这几种可能性逐项核查：先走完原模型中的形状和应变，再计算平衡尺寸分布；随后改变总量、开放更低自由能体相，最后单独改变速率。不同对照回答不同问题，不把一条能量曲线当作完整组装系综。
局部偏好的角度，为何不能全部满足
原研究在方格上安排可占据的元件。$\eta_i=0,1$表示某位置是否有元件，$\theta_i$是其内部取向；邻居结合能为
$J>0$控制偏离偏好角度的代价，$A_{ij}$是相邻元件偏好的角差。每个小方格绕一圈的偏好差为
若这个和不是$2\pi$的整数倍，实际取向绕一圈仍须闭合，所有边就无法同时达到最喜欢的角差。这是模型中明确的几何不相容。令$\Sigma=J+K>0$为达到最佳角差时每条已结合边的净收益；小$\Sigma/J$使结合可以在热涨落尺度上可逆，同时内部取向仍相对刚硬。1
这里每个元件只有一个内部角度。它不是任意三维刚性颗粒的完整拼装模型；原论文通过连续极限联系曲面取向序等物理对象，实际材料的伸长、弯曲和接触不能无条件省掉。
在无缺陷、角度偏离较小的团块内，取向场的弹性能为
自由边界满足法向相容条件$\hat n\cdot(\nabla\theta-\mathbf A)=0$。均匀错配使典型角度应变随横向尺度$W$增长，量级为$fW$；弹性能密度便按$Jf^2W^2$增长。这一累积可以抵消减少边界所获得的结合收益。
保持紧凑，与允许长成条带，会得出不同结论
先取一个$X\times Y$的无缺陷矩形，以晶格间距$a=1$为长度单位。矩形有$XY$个元件、$2XY-X-Y$条内部结合边，所以最佳结合部分每粒子的能量恰为
如果使用通常的物理周长$P=2(X+Y)$，对应线张力是$\sigma=\Sigma/2$。2025作者稿明确给出一般单位下的$\sigma=\Sigma/(2a)$；这与缺边计数一致。不能把这两个$\Sigma$和$\sigma$直接当成相同线张力，否则宽度系数会差一因子。3
2023发表论文附录C给出的矩形弹性能，可写为
其中可交换两轴，令$X\le Y$，且
这是给定矩形内取向场放松后的结果；它没有开放缺陷。本文采用连续尺度形式，未把它在$X$或$Y$接近1时的外推当作微观精确解。12
若强制$X=Y$，边界代价随尺寸下降，弹性代价随尺寸上升，每粒子能量出现有限最小值。现在固定粒子数$n=XY$，允许改变长宽比：超过一定$n$后，正方形失去优势，最省的做法变成长条，把应变梯度主要放在窄方向。
定义$\ell_d=[\Sigma/(Jf^2)]^{1/3}$。本库实际扫描形状分支，并用取向场流函数的独立双正弦展开核弹性级数。正方形的长宽比稳定性在
附近改变；原文报告约0.837，差约0.77%，本文保留自己的级数稳定性结果，没有将近似数硬调一致。固定矩形连续模型中，增长到长条极限后
从而选择有限宽度
此时继续延长主要增加相同的每粒子能量，已不再承担紧凑团块的超广延增长。几何错配选择了宽度，没有单靠这个能量固定总粒子数。这也解释了为何研究“停止生长”必须先指定所测的是半径、宽度、长度还是质量。
有限温度下，怎样得到尺寸分布
一个团块的每粒子能量，还要接到许多团块与单体之间的化学平衡。以刚性条带作透明近似，写每条$n$聚体的内部自由能为
$\Delta$是两端合计的额外代价；把一长条切成两条会多付一份端部代价，但会多获得一个平移自由度。平移熵因而倾向于多个有限条带，与端部收益相竞争。
本库从上面的矩形级数直接取长条极限，保留端部弹性放松，得到
括号约为0.684876。原附录C6另印$\Delta=\Sigma/W_*$，与其矩形能量的长条展开不一致；本篇使用实际展开，并以有限长条总能量减去$n\epsilon_\infty$交叉核，不称照抄该端部公式的复现。端部能量的具体系数仍限于此矩形近似。
设$\rho_n$为每单位系统面积中的$n$聚体数。考虑团块间相互作用与排斥被忽略的理想混合物，其自由能为
$g_n$记录内部简并，$\Phi$是总粒子浓度；长度和面积均已用$a$归一。对自由能在固定总质量下极小，得到
$\phi_n$是质量分布，与按团块数统计的$\rho_n$不同。本篇$\mu$有能量单位，因此单体满足$\beta\mu=\ln\rho_1$。原文有量纲与无量纲记号混用处，这里固定一种约定。1
对两种等能矩形方向取$g_n=2$，则
任何有限浓度下$z$都小于1，因而尾部可归一、平均长度有限。但随着浓度上升，$z$靠近1，长度尾部变长。有限平均长度来自整个混合物的平衡，不是某条带在一个粒子数上永远锁死。
为把这个过程实际算完，取原研究一组参数$J=1$、$\Sigma=0.09$、$f=0.016$、$\beta J=40$。本库得到$W_*=4.74542$、$\epsilon_\infty=-0.151552J$、$\Delta=0.292502J$。单体单独取$F_1=0$、$g_1=1$；条带从$\lceil n_c\rceil=43$开始，暂不计中间的短寡聚物。
对条带几何级数作无限求和、按总质量求化学势，得到：5
| 总浓度$\Phi$ | 单体浓度 | 按条带数平均的粒子数 | 平均长度，晶格单位 | 宽度 |
|---|---|---|---|---|
| 0.005 | 0.002220 | 63.19 | 13.32 | 4.745 |
| 0.02 | 0.002271 | 81.50 | 17.17 | 4.745 |
| 0.05 | 0.002290 | 101.17 | 21.32 | 4.745 |
| 0.20 | 0.002309 | 154.86 | 32.63 | 4.745 |
这些不是原Monte Carlo数据的拟合。模型省略短寡聚物、格点宽度修正、弯曲和支化的构象熵，以及团块间排斥。原模拟的尾部更接近支化聚合物统计，与本篇刚性线条的$n z^n$质量尾不同；原论文也报告连续条带理论与格点聚集浓度有明显偏移。不能把本例单体饱和值0.002330当成真实材料或原模拟的精确聚集阈值。高浓度结果尤其只是理想模型的延拓。1
小盒子里的一大团，不一定是宏观凝聚
刚才的无限混合物适用于总量足够大的极限。若整个系统只有100个元件，而典型条带就需要几十至上百个，末态自然容易被一个大团主导。
本库对同一个理想单体—条带模型进一步固定总粒子数$N$，令面积$V=N/\Phi$，使用精确正则配分函数递推
平衡团块密度为$\rho_n=q_n Z_{N-n}/Z_N$。为了避免数值上溢，实际计算把所有$F_n$减去同一线性项$n\epsilon_\infty$；固定$N$时，这只改变配分函数的共同因子，不改变任何尺寸概率。独立枚举小系统的全部整数分拆与递推一致，物料守恒也逐项核过。5
保持$\Phi=0.05$不变，只增大系统：
| 总量$N$ | 按条带数平均的粒子数 | 处于$n>N/2$团块中的质量比例 |
|---|---|---|
| 100 | 83.22 | 85.13% |
| 500 | 97.43 | 5.09% |
| 2000 | 100.25 | $7.08\times10^{-8}$ |
| 无限理想混合物 | 101.17 | 宏观质量比例趋于0 |
典型条带大小趋向有限值，但小系统中完全可能看起来“几乎都在一个大团里”。这不是新的体相转变，而是总量尚不足以同时容纳许多典型组装体。判断宏观凝聚要看扩大系统后，一大团是否持续承担非零的总质量比例，而不能只看某一张图。
上述计算仍是无排斥的理想混合物，不等于原二维格点系统的正则模拟。它的作用是把有限总量这一替代解释实际走完，显示固定浓度和固定总量是两个不同控制。
开放缺陷或连通后，有限条带还会不会胜出
紧凑结构积累应变，可以借形状变细来缓解；它也可以借缺陷改变取向场的绕转。原2023模拟的低错配体相不是无缺陷的大正方形，而是含带荷孔洞的“涡旋海绵”：孔洞附近的取向绕转屏蔽背景错配，避免把无缺陷的超广延代价一路付到无限大。1
一个有限形状族的最低能，不保证它比这个体相更稳定。为明确这一点，本库在理想条带混合物之外，加入一个单独指定的宏观体相自由能$g_B$。这是敏感性对照，未从原海绵的构形积分或模拟求得$g_B$。
若$g_B<\epsilon_\infty$，有限相的化学势最多升到$\mu=g_B$；此时$z_B=e^{\beta(g_B-\epsilon_\infty)}<1$，有限团块和单体能容纳的总浓度有上限
超过它的质量进入宏观相。以同样条带参数为例，若假设$g_B$仅比$\epsilon_\infty$低$0.001J$/粒子，则$\Phi_{\mathrm{sat}}=0.00735136$；总浓度0.05时，约85.30%的质量进入宏观相。有限条带宽度的原计算仍然成立，但已不能独自保证整个系综以有限条带为主。
反之，若体相自由能不低于$\epsilon_\infty$，这个理想条带系综的容量在$\mu\to\epsilon_\infty^-$时发散，不会因这一指定体相而在有限浓度凝聚。这两种结局由实际允许状态的自由能竞争决定，不是由“有错配”三个字决定。
入库时登记的2025预印本，现已核为同年12月15日发表的Physical Review E 112, 065419。本文取得固定作者v1全文与附录，发表版全文请求返回401，因此不声称逐式核过最终发表版。34
该研究区分两条变化路径：一条在较强错配下增加浓度，有限宽度结构先连成跨系统网络，再逐渐填成孔洞体相；另一条在较稀浓度下减小错配或降低温度，转向相分离的含缺陷体相。前者的连通变化，不能直接等同于后者的凝聚自由能平衡。
作者稿还比较了条带的平移与弯曲构象熵。较多有限条带分摊较少元件的平移自由度，能够比一个宏观团块获得更大的每粒子熵；这可以使有限条带在有限温度下比单纯零温能量预期更稳定。其相界均场估算使用了匹配模拟的系数，并不是无参数的精确定理。原稿圆形孔洞胞的部分缩放记号另有不一致，本文未用其字面表达报告精确相界，也未运行原作者的大格点模拟。
展平、脱离基底、改变层数或进一步组装成束，是其他材料可能开放的逃逸路径。这个方格内部角度模型没有自动包含它们。允许哪个资源改变，就需要把那个状态加入同一比较，而不能把无缺陷、固定形状的极小延伸成全部材料的自限保证。
怎样区分平衡自限、供料上限与动力学卡住
原论文对平衡的支持包括可逆取向和粒子移动、团块移动、组装自相关及有限系统尺寸比较。含缺陷体相尤其容易陷入亚稳态，因此作者使用非局部团块移动加速采样；这些Monte Carlo步是平衡采样工具，不是胶体实验的真实时间。12
为了把“平台不等于自限”做成可检查的对照，本库另建一个单团块生灭模型。粒子数$n$在允许区间内每次增减1，给定有效势$G(n)$，正反速率满足
因此平衡概率$P_{\mathrm{eq}}(n)\propto e^{-\beta G(n)}$。把某条通道的正反速率同时乘以同一个小因子，不改变这个概率，却可以任意拉长等待。
取自选的有限势$G(n)=0.8\sqrt n+0.002n^2-0.25n$、$\beta=4$。从$n=1$和$n=150$两个相反初态出发，最后都到平均47.8486的同一有限分布。这个模型里有真正随尺寸增长的惩罚，不靠总量上限维持47.85；允许区间取到200，概率峰远离上端。
再改成偏向大团的$G(n)=0.8\sqrt n-0.08n$。正常速率下，它趋向由200上界截断的大团分布，平均195.582。只把$70\leftrightarrow71$通道的两向速率同时乘$10^{-6}$，平衡完全不变，但到约$10^5$的模型时间仍停留在约33的平台；延长到$10^9$后，实际计算趋向同一个195.582分布。
若保留正常速率，却把总可用粒子数限制在70，同一势的平衡平均是32.6052。它与慢速率过程的中途平台很接近，但原因是可用状态被截断；增加供料上限就会改变答案。总量为70也不意味着所有粒子必须在同一个70聚体里。
这些是严格声明输入的单团块例子，未描述多团块相遇、物料池反馈或原XY模型的微观动力学。其用处是给出反例：相似的有限平台可以对应不同的平衡和不同的速率机制。 数值中极小初态概率使直接谱展开失稳，过长的直接积分也不适合长期验收；本库留快照后改用短时稀疏积分再接衰减模态，另用矩阵指数、详细平衡和物料守恒核对，未把计算失败当作物理停滞。5
在实际研究中，应同时看三个变化：固定浓度扩大系统，检验典型尺寸与宏观质量比例；改变浓度或供料，检查长的是宽度、长度、团块数量还是仅耗尽附近物料；从不同初态、延长等待或加快可逆重排，检查是否趋向同一分布。温度变化还会同时改变熵、结合可逆性与速率，不能只解释成“更容易越障”。
这篇能给出的答案
明确几何错配可以通过累积应变与边界结合竞争选择有限尺度。但在原格点研究里，最直接被选择的是宽度；长度的有限均值还依赖浓度和熵。开放缺陷体相、改变形状或增加系统连通后，有限条带未必仍是整个系综的稳定状态。
本文实际增加了四项对照：完整矩形形状放松揭示紧凑最优尺寸的逃逸；理想条带无限分布与有限总量的精确统计区分宽度、长度和小盒子效应；显式体相自由能检验有限形状极小是否足够；守详细平衡的长期演化显示相同平台如何来自供料或慢速率。这些是有边界的模型结论，未宣称完整原相图或新材料实验复现。
原二维模型的逐构形数据和执行代码尚未取得，支化、弯曲和排斥下的完整自由能分布未重算。2025最终发表版全文未取得，真实材料的制造偏差、基底变形和次级聚集也未在同一模型中量化。要给特定组装体报可靠尺寸，还需把这些允许状态、热涨落和观测时间接回其实际材料参数。
来源与复核
<a id="source-s1"></a> [1] Hackney, Amey & Grason (2023). Dispersed, Condensed, and Self-Limiting States of Geometrically Frustrated Assembly, Physical Review X 13, 041010。发表全文已取得，核模型、宽度/质量/形状统计、含缺陷凝聚、附录A采样与自相关、附录C矩形与理想聚集。原大格点Monte Carlo未执行，原分布未数字化成拟合数据。
<a id="source-s2"></a> [2] 同研究补充材料。 出版方SI与固定作者v2的归档SI均取得，SHA256一致。核参数表S1—S11及附注1—2的矩形解、格点截止与分布近似；表是模拟参数，不是逐点轨迹数据。图S1说明为单快照，未当独立样本平均。
<a id="source-s3"></a> [3] Hackney & Grason，固定作者2508.21688v1 (2025-08-29)。 Geometrically Frustrated Assembly at Finite Temperature: Phase Transitions from Self-Limiting to Bulk States。作者全文及附录已取得，核§III线张力/条带、§IV浓度贯通、§V平移/构象熵及均场适用范围。内部角度绕转的涡旋与T128三角晶格的配位缺陷分别定义；未称字面相同。
<a id="source-s4"></a> [4] 同研究发表信息。 Physical Review E 112, 065419，2025-12-15，出版方题录/摘要已核并保留网页。发表版全文请求401，未称已读；本篇方法核验使用上面固定作者稿，不把它继续称为未发表研究。
<a id="source-s5"></a> [5] 本库实际计算 (2026)。 完整输入和结果、分布和时间曲线输入。research/calculations中的self_limiting_frustrated_assembly.py、_audit.py、_figures.py实际执行。矩形单/双级数、形状密扫、端部展开、无限分布、正则递推与整数分拆、长短时生成元/矩阵指数、详细平衡和供料/体相收支均核。自建对照与原作者MC、理论、发表版本分别记账。