形成与自组织15 分钟

局部颗粒怎样拼出有长程秩序、却不周期重复的准晶?

准晶形成不仅需要局部图元拼得上,还要让它们在生长和重排中维持整体关联。本文重算2025年胶体坐标与2024年振动颗粒的取向数据,并按2026年机械组装研究的同类力学条件,独立计算低压圆盘的周期近似结构和19粒图元。局部有序指标很高不独证准晶,指标很低也可能来自不同边族抵消;形成路径、理想模板、局部机械稳定与有限温度的自由能选择须分开。胶体有效双井的来源与外驱体系的稳定机制仍有未决部分。

局部规律可以不断出现,整张图却没有重复单元

普通晶体可以用一个有限的基本单元,沿若干方向不断平移复制。准晶也有规律,局部图元甚至会反复出现,但在本篇讨论的理想二维无限结构里,没有一个非零的固定平移能让整张图完全重合。

因此,“不周期重复”不是“什么图案都不重复”,也不是“颗粒位置完全随机”。准晶中仍有跨较大距离的排列关联,能够产生规则的衍射峰。它把位置有序与平移周期这两件事分开了。1

形成这样的结构,不要求每颗粒子先看见一张无限大的图纸。局部接触可以限制下一层怎样放,已有边界又把这些限制向外传递;后续重排则让不兼容的接触改变。但这条过程链有条件:某种局部图元能形成,未必能扩展成准晶,也可能接成周期晶体或停在混合结构中。

本篇用三个尺度相接的案例理解这一点:能直接观察的胶体重排、毫米颗粒的振动组装,以及2026年提出的机械重排路线。

三角形与方形,都能围住同一个顶点

先看十二重准晶常见的一类几何:用等边三角形和正方形铺成平面。

六个三角形可以围成360°;四个正方形也可以。还有一种组合是三个三角形和两个正方形:

$$ 3\times60^\circ+2\times90^\circ=360^\circ. $$

不过,相同数量还可以有不同排列次序。围着一个顶点,三角—三角—方形—三角—方形,通常记作$\sigma$环境;三角—三角—三角—方形—方形,则是H环境。六个三角形围成的邻接称Z环境,四个方形称S环境。2

这些名字描述局部连接方式。图像里的多边形可以稍有形变,不必每个角都刚好等于理想角度。更关键的是,$\sigma$或H环境自身也能周期铺开,所以找到它们只是确定了可用的拼接零件,还没有确定整体属于准晶。

2025年的胶体快照:结构怎样一片一片改接

Gao等2025年报告,在同时存在电场和磁场的体系里,单组分微球能形成二维十二重准晶结构,并观察到从其他排列转变的过程。2

本篇依据其补充材料的共同邻居分类方法,重新分析了带明确时刻的公开坐标快照。下图中,绿色是Z环境,粉色和黄色分别是$\sigma$和H环境;其他颜色标出方格、两类拓扑缺陷或暂未归类的位置。3

公开胶体坐标的重新分类。该组转变由边界一侧推进,时间使用源表标记。邻居距离阈值为颗粒直径的1.25倍;灰色是本篇按凸包边缘距离排除的边界颗粒。颜色是连接拓扑分类,不是逐像素复制原图。
公开胶体坐标的重新分类。该组转变由边界一侧推进,时间使用源表标记。邻居距离阈值为颗粒直径的1.25倍;灰色是本篇按凸包边缘距离排除的边界颗粒。颜色是连接拓扑分类,不是逐像素复制原图。

在本篇保留的内区,17.5秒时约96%的颗粒是六角邻接的Z环境;到71.7秒时降到约8%,$\sigma$与H环境合计则从约2%增到约92%。新的连接方式从一侧进入,在前沿附近伴随缺陷,再逐渐扩展。

另一组快照显示了从内部出现新结构的路径。它的早期缺陷较多,之后混合邻接增加,观察终点仍留下较多未完全整理的部分。

两组公开快照各自的环境份额。左侧为边界推进,右侧为内部成核;它们是不同条件下的观察,不能据时间轴直接给两条路线排速度名次。每个快照独立分类,没有把颗粒数当作独立实验次数。
两组公开快照各自的环境份额。左侧为边界推进,右侧为内部成核;它们是不同条件下的观察,不能据时间轴直接给两条路线排速度名次。每个快照独立分类,没有把颗粒数当作独立实验次数。

这把生长过程具体化了:颗粒原先的邻居关系需要拆开、改接,前沿才能前进。局部接上新的图元之后,也还需要继续整理缺陷与取向。

这里的定量分析使用坐标快照。部分另附比例表存在时间列长度不齐、负比例等问题,本篇保留了这些空位,没有自行补时间或把负数改成零;边界定义和邻居阈值也另作了敏感性检查。

有效作用能解释多少,还要回到具体物理来源

胶体中的外场能够改变颗粒之间的磁作用、电作用以及与液体有关的作用。它们影响局部邻距偏好、颗粒能否重新调整,以及不同结构之间的转变。

Gao等提出一个重要假设:具有两个有利距离的有效双井相互作用,可以帮助稳定他们观察到的准晶。不过,论文明确保留了第二个势阱的物理来源问题。2

本篇核到的正式补充材料也讨论了多种可能的场致液体运动。作者的公开模拟代码主要展示小团簇的改形,例如从菱形或链状结构变为方形、环及部分局部图元;这些结果没有把“整个十二重准晶如何由已确认作用产生”的问题全部关闭。4

因此,较稳的结论是颗粒确实能重排出相应结构,并且外场条件能改变过程;将全部效果归结为一个已经独立测定、完整有效的二体势,仍需更多证据。外驱条件下使用“有效温度”等描述,也需要保留它与真实热浴温度的区别。

本篇取得了正式补充材料、公开摘要与图注、坐标数据和代码;未取得2025年正式主文全文。没有用这些公开部分冒充对全部主文论证的核查。

毫米颗粒的案例:三维接触改写了二维拼接条件

Plati等2024年研究了另一种体系:大小两类球形颗粒放在振动底板上,观察到八重准晶序。颗粒的运动主要由外部振动供能、碰撞和摩擦耗散维持,不能直接当成平衡热运动。5

这里有一个看似很小、实际很重要的几何区别:大小球都接触同一底板,但球心高度不同。从上方看,两球的水平中心距,可以小于把它们当作真正二维圆盘时的接触距离。

设两球直径为$D_L,D_S$。它们接触时,三维球心距为$(D_L+D_S)/2$,高度差为$(D_L-D_S)/2$,所以水平中心距是:

$$ d_{LS}=\sqrt{D_LD_S}. $$

当小球直径是大球的一半时,这个距离为大球直径的0.707倍;二维可加圆盘对应0.75倍。它们在三维中没有互相穿透,只是二维投影不能照搬同平面圆盘的接触规则。

这种几何改变了小球能进入什么空隙、哪些大小颗粒组合能够共同存在。实验里形成了包含不同颗粒组合的两类方形和一类等腰三角形图块,再通过重排逐渐出现八重取向关联。

左图为球心高度差导致的投影接触关系。右图是几类理想三角—方形接触铺砌的面积核算,用于说明局部图元不独自决定最密排列;右图不是对左侧二元球体系的自由能排序。
左图为球心高度差导致的投影接触关系。右图是几类理想三角—方形接触铺砌的面积核算,用于说明局部图元不独自决定最密排列;右图不是对左侧二元球体系的自由能排序。

从几何上看,单组分理想三角密排仍比这些三角—方形混合铺砌更紧密。现实体系为什么选准晶,必须继续考虑组分、作用和可访问的排列方式,不能只停在“这样装得密”一句话。

图块先出现,整体对齐可以慢得多

振动颗粒研究的一项价值,是可以长时间直接观察重排。作者区分了局部图块形成与图块之间的取向整理。公开源表覆盖0—190小时;其中与最终主取向对齐的图块面积份额,从约0.147增加到约0.580,但中间多次下降和回升。56

这个“对齐”以最终主方向为参照,是一种回看式统计。曲线不表示系统从第一刻就已经选定那个方向,而是在描述局部结构怎样逐步接成更大的关联区域。

振动颗粒公开数据重绘。左图仅使用明确的对齐面积列,未把源表中重复的两列解释成独立类别。右图分别计算长、短边方向的八阶角度指标;初、中、末三帧的边数量不同,不是三千多个独立实验。
振动颗粒公开数据重绘。左图仅使用明确的对齐面积列,未把源表中重复的两列解释成独立类别。右图分别计算长、短边方向的八阶角度指标;初、中、末三帧的边数量不同,不是三千多个独立实验。

右图还说明,低有序指标也可能来自统计对象混合。末期长边和短边各自的八阶指标约为0.497和0.638,把两类直接混合却只有0.010。原因是两类边的主方向相差约22.4°;在八阶复数表示中,这几乎是半圈,两个有序贡献会相互抵消。

所以,在这个体系里,长边与短边需要按它们各自的几何角色解释。不能见到混合平均接近零,就认定所有方向都随机了。

至于为什么这个受驱动、持续耗散的体系在某些条件下接近平衡模型的结构,原研究也没有给出完整的一般理论。相似的最终图样不意味着能量注入、碰撞耗散和不同颗粒的运动统计都可以忽略。

高“十二重指标”,为什么也不能独证准晶

对一个颗粒,把它指向邻居的连线角度记作$\theta_j$,常用局部指标是:

$$ \psi_m=\frac1{n_b}\sum_{j=1}^{n_b}e^{im\theta_j}. $$

它检查邻居方向是否与某种$m$阶对称相容。但12既是4的倍数,也是6的倍数,所以完美方格和完美六角邻接都能给出$|\psi_{12}|=1$。

本篇的实际计算进一步显示,一个周期性的$\sigma$结构,也可以有几乎为1的十二阶局部指标。

| 构型或邻接 | 局部$|\psi_4|$ | 局部$|\psi_6|$ | 局部$|\psi_{12}|$ | |---|---:|---:|---:| | 理想方格 | 1 | 0 | 1 | | 理想六角 | 0 | 1 | 1 | | 本篇算出的周期$\sigma$结构 | 约0.20 | 约0.20 | 约1.00 | | 胶体准晶片区内区平均 | 约0.25 | 约0.37 | 约0.94 |

这不是说周期结构“比准晶更准晶”,而是这把尺子只测到了一部分信息。

还要区分两种平均:先取每个颗粒的有序度大小再平均,会保留各处的局部规则性;先平均复数再取大小,则会受到不同区域方向抵消的影响。两种信息都有用,但都不能单独判定全局是否存在重复晶胞。

上排为420颗胶体坐标,下排为本篇算出的周期近似结构复制成2×2片区。右侧统一用标准结构因子S的log(1+S)显示,屏蔽中央直射峰;复制用于清楚显示周期性,不是额外的大规模组装实验。
上排为420颗胶体坐标,下排为本篇算出的周期近似结构复制成2×2片区。右侧统一用标准结构因子S的log(1+S)显示,屏蔽中央直射峰;复制用于清楚显示周期性,不是额外的大规模组装实验。

因此,研究要把局部图元、区域取向、衍射峰位置与精细结构、位置关联及周期候选放在一起看。多个普通晶畴的叠加,也可能留下看似高重数的峰;2025年的一项软颗粒理论研究就展示了这种具体竞争结构。9

有限视野能够排除许多简单解释,并建立明确尺度内的准晶序证据,但不能仅靠一张小图,证明无限结构绝无任何平移周期。周期边界模拟也需要借助不同尺度与竞争构型来理解其块体极限。

2026年的机械路线:先看局部图元能怎样长出来

Jiang等2026年提出,利用方格排列的力学不稳定性,可以让圆盘通过合作重排形成复杂结构。7

方格的水平和竖直接触,并没有同样充分地约束所有方向。某些协同位移可以在开始时很少拉长主要接触;加入压缩条件后,这些柔软方向可以成为降低焓的路径。局部方块转动,邻近空隙被挤成三角形,三角—方形混合邻接便出现了。

本篇按原文同类的短程谐和排斥势与小正压条件,独立做了两个计算。粒径和相互作用能量单位设为1,相互接近时的势能为$U(r)=(1-r)^2/2$,分开到$r\geq1$时为零;在周期方盒中最小化:

$$ H=\sum U(r_{ij})+pL^2, \qquad p=10^{-4}. $$

如此小的压力下,圆盘只产生很小重叠,接近硬颗粒极限。这里使用本篇的优化与曲率检查,未重跑作者的完整算法和全部大体系。

第一组从100颗圆盘的方格初态出发,轻微扰动后,得到全部由$\sigma$邻接构成的周期结构。第二组从81颗圆盘的方格中移去一对指定位置,剩余79颗重排后,出现了中心一颗、内圈六颗、外圈十二颗的19粒图元。

本篇独立力学计算。左侧是周期近似结构;右侧的星号与虚线标示局部19粒图元及两层几何壳,周围仍有其他结构。两组都不是完整的非周期准晶。
本篇独立力学计算。左侧是周期近似结构;右侧的星号与虚线标示局部19粒图元及两层几何壳,周围仍有其他结构。两组都不是完整的非周期准晶。

曲率检查给出了进一步区分:未扰动的100颗方格在此正压下有18个使焓下降的扰动方向;重排后的两组结构,去掉整体平移的两个中性方向,其余曲率均为正。这里固定盒子的方形形状,允许颗粒移动及盒长等比例变化;曲率结论只覆盖这些扰动,没有检验任意改变盒形的剪切。它支持这一约束下的局部机械稳定,但没有证明它们是全部可能结构中自由能最低的状态。

19粒图元也需要看完整:虽然外圈有十二颗,内圈只有六颗,整个图元只有六重对称。两种相差30°的图元能进一步组合,但“做出一个十二边形外圈”与“形成整体十二重准晶”仍然不同。

原研究对硬颗粒近似的结果也有清楚边界:增加缺位后,有些图元会被六角排列吞并,没有形成全局十二重准晶。它不是“单一尺寸硬圆盘终于普遍长成准晶”的证明。

局部能生成以后,还需要把它们维持和连接起来

为了得到更大尺度准晶序,2026年研究转向了已知适合准晶形成的双尺度排斥相互作用。这里的两个尺度来自核心和肩部,不是两个吸引势阱。

在这种模型中,适当分布的缺位可以触发多个局部图元形成,原先的结构也能在更大范围重新连接;之后的热运动或论文所报告的小幅剪切,可以帮助不够规整的区域继续调整。7 本篇没有独立复算这部分大体系的优化或剪切过程。

研究中最接近理想高层铺砌的方案,还使用了已知理想图样来安排缺位位置。本篇核对公开文件:349个模板中心对应移去698颗,原6400颗方格剩下5702颗,调整后的数密度为0.94。8

这条路线的贡献,在于把已知布局转成可由机械重排实现的结构。全局位置约束已有一部分来自模板,局部颗粒怎样落位与放松则由相互作用决定。它与随机布置缺位后形成随机铺砌的路线,回答了不同问题。

数学上的理想铺砌提供图纸,局部力学提供执行这张图纸的一种可能过程;是否不依赖这种图纸也能稳定生长,还要另外验证。

能形成、能停住、热力学上优选,是三种判断

局部图元形成,说明有一条可达路径;小扰动后保持,说明局部机械稳定;在指定温度等条件下自由能更低,才涉及平衡选择。受驱动体系的长期维持,还要解释持续输入和耗散。

Li等2025年研究广义Hertz软排斥颗粒时,找到了一批比此前模拟中的准晶更低能的周期候选结构。模拟中较高温度的准晶状态,在进一步冷却时可以转为复杂周期晶体或其共存结构。9

这并不意味着有限温度下的准晶没有稳定价值。自由能不仅包含势能,也涉及不同排列、振动等所提供的熵。多个能量很接近的拼接方式,可能在有限温度下共同影响选择。1

不过,有限温度模拟的观察与低温能量比较,并没有替本篇完成所有构型的自由能计算;该Hertz模型也不同于前述胶体、振动颗粒和双尺度肩势,不能直接用其排序替其他体系判定稳定性。

某些改变准晶允许拼接关系的重排,被称为相位子(phason)重排。它与仅把原有排列稍微弹性挪动不同,可以改变局部拼法及对应结构统计。但真实材料中的划分依赖模型和尺度,一次普通颗粒跳动不能自动称为相位子翻转。1

因此,重排既可能帮助已有图元连接、修复缺陷,也可能把结构送往一个周期竞争相。知道“动得起来”以后,还要看相互作用和当前环境把这些可行路径偏向哪里。

局部规则怎样积累成整体约束

这几组研究共同提供了一条具体过程链:几何与相互作用先限定哪些邻接容易形成;成核位置和已有结构影响它们怎样继续拼接;集体重排让不同片区调整方向、改换接触,最后形成在观察尺度内能够维持的排列。

其中有多种可行路线。不同粒径和高度可以改变有效接触;多个作用尺度可以支持混合图元;机械不稳定性可以打开快速重排路径;模板、边界与有限温度的构型选择,又会改变最后留下什么。

局部图元是整体秩序的材料,生长和重排才把它们接成结构。 准晶的有序性并不需要周期复制,但它的形成也不会由一个局部对称标签自动保证。较可靠的研究会同时检查拼得上、接得远、改得动,以及在什么条件下能够留下来。

可回查:胶体快照统计、逐颗坐标与分类、邻居阈值敏感性、振动颗粒对齐面积、长短边角度统计、100颗周期近似结构、79颗含图元结构、指标比较、投影接触几何、计算汇总与限制。

来源与依据

  1. S1 领域综述

    Quasicrystal Structure Prediction: A Review

    Michael Widom; Marek Mihalkovič · 2024

    Israel Journal of Chemistry64,e202300122;2023-10-16在线,2024卷期;作者公开发表PDF · 相关全文已核

    Introduction/结构预测与周期近似、有限模拟、自由能/振动熵/相位子段落

    区分数学铺砌、结构预测与稳定性;不把高维记账当额外物理空间,也不把所有重排叫phason。本文未复算金属DFT/自由能。

  2. S2 原始研究公开部分及正式补充信息

    Direct observation of colloidal quasicrystallization — public article material and Supplementary Information

    Yan Gao; Brennan Sprinkle; David W. M. Marr; Ning Wu · 2025

    Nature Physics21,966–973;2025-03-31;原始SI随文PDF · 摘要已核

    正式摘要/图注;SI§1结构因子、共同邻居、取向与phason;§4数值示例、§5电动流讨论及影片说明

    按主文读取范围标abstract,SI相关部分已完整取得核对。正式主文全文未得;第二势阱来源未决原摘要明确。SI原S与I平方、相位角和全局局部模的口径不盲用,本库明示自身定义;不把电动流概括为已确认纯二体势。

  3. S3 公开处理坐标/表格

    Direct Observation of Colloidal Quasicrystallization — Figure 2 coordinates and summary sheet

    Yan Gao and coauthors · 2025

    Figshare四工作簿之一;作者已识别颗粒的位置,不是原视频 · 相关全文已核

    coordinate for figure2b四时刻,标准化直径坐标;fig2c仅检查数据完整性

    按SI共同邻居法与cut1.25D重算四快照,另取凸包内1.25D内区;不是原整视野百分比复现。fig2c有2869比例行但仅718时间行、690个负值,不自行补时或钳零,不用其作时间轨迹。

  4. S4 作者模型范围说明

    Colloidal Quasicrystal numerical examples — README

    Brennan Sprinkle · 2025

    · 相关全文已核

    Square_to_Diamond、Ladder_to_Ring、Rhombus_to_Sigma三组示例及分阶段/种子说明

    发表链接QuasiCrystal大小写与当前Quasicrystal路径不同,已核目录定点。只核局部团簇示例范围,未执行外部代码或完整大尺度准晶模拟;不以能做σ团簇声称完整稳定势已建立。

  5. S5 原始颗粒实验与碰撞模型

    Quasi-crystalline order in vibrating granular matter

    A. Plati; R. Maire; E. Fayen; F. Boulogne; F. Restagno; F. Smallenburg; G. Foffi · 2024

    Nature Physics20,465–471;2024-01-19。实际正文arXiv2307.01643v2(2024-11-27),不宣称与正式版逐字相同 · 相关全文已核

    原始作者稿几何/非可加接触、碰撞规则、Fig1–3与长时重排讨论;正式数据/代码入口

    毫米、二组分、八重,非胶体十二重;大/小球高度不同导致平面投影非可加。非平衡供能/耗散不等同热平衡;未重跑5.3e10碰撞或DEM,SI未取得不称全复现。

  6. S6 正式公开源表

    Plati et al. — Source Data Figure 2

    A. Plati and coauthors · 2024

    · 相关全文已核

    Fig2a对齐面积;Bond orientation histogram初/中/末长短边;结构因子说明只核口径

    447个时间点0–190h,对齐面积.14659→.57980。misaligned/defect两列完全重复且三列和可>1,未用作三类分账。分别重算R8,末期长.49742短.63809、混合.01038,主方向差22.4148°;不把边数当独立实验。

  7. S7 原始力学/组装模拟研究

    Rapid and high-quality formation of dodecagonal quasicrystals and their approximants using a purely mechanical approach

    Zhehua Jiang; Jianhua Zhang; Mengyuan Zhan; Jiaqi Si; Junchao Huang; Hua Tong; Ning Xu · 2026

    National Science Review13,nwag244;2026-04-24在线;PMC13278518 · 相关全文已核

    Results各阶段、Methods两种势/能量焓最小化、Fig1/2图元与周期近似、Fig3–6 TLS/模板及Discussion

    硬粒极限近似为低压谐和排斥,不是真正无限硬势。硬粒结果止于周期近似和图元,完整DDQC使用已适合准晶的TLS势;理想高层结构按图样放缺位。剪切与大体系仅核主文报告,未独立复算;未取得SI包。

  8. S8 公开模板与初态说明

    Simply-form-DDQC — public initial-configuration documentation

    Jiang et al. research group · 2026

    · 相关全文已核

    01为Fig6a理想hyper-particle中心;02为移去指定缺位的方格初态

    两数据文件头分别349/5702,初始6400,恰好移去698=2*349,rho.94。支持模板信息已进入初态,不把机械落位等同从全随机发现理想全局布局。

  9. S9 原始周期候选能量与退火模拟

    Hierarchical ground-state crystals underlying Hertzian quasicrystals

    Yao Li; Yiwei Wang; Yingke Geng; Wenyu Liu; Fangfu Ye; Jeff Z. Y. Chen · 2025

    Nature Communications;2025-11-25;PMC12717090 XML · 相关全文已核

    Introduction广义Hertz势;Results能量枝/相共存/退火与Fig3衍射歧义;解释与限制

    较低能周期候选挑战把特定模拟QC认作基态;有限温度熵稳定解释未由本库自由能复算证实。本文不宣称已穷尽全部构型;Hertz模型不替胶体或TLS判稳定,非同源数据重复。

  10. S10 公开处理坐标

    Gao et al. — Figure 1 particle coordinates

    Yan Gao and coauthors · 2025

    · 相关全文已核

    420颗、直径标准化列;CNA/psi与惯用S(k)重新计算

    本篇内区316颗,σ215/H73/Z24/S2/D+1/D-1;mean|psi12|=.94039。边界与cut口径明确,未从有限S图单独证明无限非周期性。

  11. S11 公开处理坐标

    Gao et al. — Figure 3 snapshots and summary sheets

    Yan Gao and coauthors · 2025

    · 相关全文已核

    fig3_a六快照,1.8/2.1/3.5/7.2/10/14.5;其余表仅核可用性

    主要用坐标独立CNA;标记比例表有负值,未静默修正。十个快照不是十个独立实验,亦无粒子身份配对。

  12. S12 公开模型输入

    Jiang et al. — ideal template centers

    Jiang et al. research group · 2026

    · 相关全文已核

    349中心、盒长、二维位置

    与README和移去缺位后的5702粒初态同核;本篇未用它运行完整6400粒TLS过程。Git blob身份见下载记录。

  13. S13 公开模型输入

    Jiang et al. — square lattice with selected vacancies removed

    Jiang et al. research group · 2026

    · 相关全文已核

    5702粒坐标、相同盒长、密度与缺位数

    不冒充最终已经松弛的准晶坐标;这是含模板信息的初态。

最后更新 2026-09-12

来源