形成与自组织14 分钟

晶体怎样从无序的材料中形成,中间态如何改变最终结果?

材料聚到一起、排列有序和选定晶体结构,是相互联系却不能混为一谈的过程。双组分胶体实验显示,作用强度和变化历史能使同一类颗粒走向直接成晶、团簇后成晶或观察期内的无序聚集。本文重绘原始数据、独立复算盐扩散,并用首次到达模型说明中间态既可能加快也可能拖慢结晶;无定形局部结构的预测方法则进一步揭示了可达路径与最低能结构之间的区别。

聚到一起、排成晶体、最后留下什么,是三件事

把能够相互吸引的颗粒放进溶液,并不意味着它们会直接排成整齐的晶体。它们可能先形成稠密但无序的团块,再在内部出现晶核;也可能迅速黏成不规则聚集体,在观察期内迟迟排不好。即便最终都是同一种晶体,中间经过的路也未必相同。

因此,至少要分开三个问题:颗粒有没有聚到一起,邻居之间有没有形成规则排列,以及形成的究竟是哪种晶体结构。

本篇以2025年的双组分胶体实验为主要案例,研究相互作用和环境变化怎样改变路径,再比较无定形前驱态的结构预测。这里的“结果”包括首次出现的晶体、观察期末的结构以及所需时间;它们不自动等于无限久后的热力学平衡。12

中间态为什么有机会起作用

晶体内部规则排列有时可以降低自由能,但形成一个很小的晶核还要付出界面和结构重排的代价。刚开始,新增有序区域带来的收益未必足够抵消这些代价,所以系统即便具备结晶的驱动力,也可能等很久才跨过去。

中间态会改变这段路。例如,先形成较密的聚集体,可以让更多颗粒靠近,给局部有序化创造条件;已有的一些排列方式,也可能更容易接续某种晶体结构。然而,颗粒如果被抓得太紧,重新排列所需的移动又会变难。1

这里没有“多一步一定更快”的保证。中间态既可能提供容易继续转化的起点,也可能让物料停留在不易离开的状态。

“两步”也不自动意味着经典理论全部失效。每一步仍可能需要跨越自己的成核势垒,并可在合适条件下用经典或连续理论描述。需要检查的是模型保留了哪些结构和动力学变量,而不是只按步骤数量给整个过程贴标签。

2025年实验:作用更强,并不一路通向更好的晶体

Zang等人研究带相反电荷、表面带有聚合物刷层的胶体颗粒。电荷相互作用与刷层排斥共同形成颗粒间的有效作用;通过改变盐浓度,可以改变作用的范围和强度。2

在这套具体体系中,减少盐浓度会增加德拜长度,并增强相反电荷颗粒之间的吸引。这个方向依赖所研究的相互作用,不能直接推广到所有盐溶液中的材料。

连续成像显示,一些条件下会先出现较密、没有晶体长程排列的团簇,随后晶核在其中形成,结晶前沿继续推进。三维共聚焦成像和固定样品的电子显微图像,帮助区分了无序聚集体与已经形成的晶体。

作者公开表中的39组实验条件。每个点保留初态和48小时状态的分类,横轴为单对颗粒势阱深度的绝对值;它不是整个晶核的成核势垒。点之间没有插值成普适相界,初态观察约在制样后5分钟。
作者公开表中的39组实验条件。每个点保留初态和48小时状态的分类,横轴为单对颗粒势阱深度的绝对值;它不是整个晶核的成核势垒。点之间没有插值成普适相界,初态观察约在制样后5分钟。

公开表给出了一幅清楚的过程图景:作用较弱时,颗粒可以保持稀分散;某些较窄条件下出现作者归类的直接结晶;再改变条件,则可出现聚集后形成单晶或多晶;很强的吸引对应观察期末仍无序的聚集状态。2a

因此,真正需要配合的是相遇、结合与重排。吸引太弱,形成稳定有序区域可能很难;吸引增强,颗粒更容易聚集,但错误接触是否还能松开、颗粒能否调整位置,也成为限制。

图中的分类是有限条件和有限时间的观察,不能把“48小时仍未结晶”说成永远不能结晶。初始画面已经出现晶体,也不能独自排除在更早、更小尺度上存在过短暂中间结构;判断路径还要结合连续成像。

此外,论文中的“气相”是胶体颗粒稀疏分散的状态,背景仍是液体;“CsCl-like”等名称描述的是颗粒排列的相似性,不能据此把胶体颗粒当成真实的氯化铯离子。

势能降低了,结构可能还没排好

团簇出现时,颗粒获得更多相互作用,势能可能明显下降。但这还不等于形成晶体:很多无序接触也能降低势能。

原研究同时使用结构指标来检查邻居排列。下图左侧是作者两条模拟轨迹的势能,右侧是不同构型中的局部取向有序指标分布。22a

重绘原始模拟数据。左侧两条轨迹的体积分数、德拜长度及电位条件不同,不能据此作同条件速度排名;步数不是实验分钟。直接途径原表有一处时间标签回跳,保留数据并在那里断开连线。右侧柱份额因舍入略偏离1,绘图时各自归一化;中途构型不等于已证实的第三种热力学相。
重绘原始模拟数据。左侧两条轨迹的体积分数、德拜长度及电位条件不同,不能据此作同条件速度排名;步数不是实验分钟。直接途径原表有一处时间标签回跳,保留数据并在那里断开连线。右侧柱份额因舍入略偏离1,绘图时各自归一化;中途构型不等于已证实的第三种热力学相。

取向有序指标大致描述:从一个颗粒指向邻居的连线,是否形成某种规则的方向组合。这里的“键”是几何邻接,不必是真实化学键。它需要与具体晶体结构、邻居定义及平均方式配合,不能拿一个阈值识别所有材料的晶化。

对于本组模拟,两步路径可以先在浓集状态停留,再进一步有序化。要分辨这两段过程,势能、密度与结构读数各自提供不同信息。

密度和有序化究竟谁先变化,也没有必要预设统一答案。Russo和Tanaka在硬球模型中强调局部取向预有序的作用;Berryman等人随后对硬球成核过程的分析,则报告有序化与致密化在所研究尺度上同步。后者认为,探测指标的灵敏度会影响我们怎样描述先后顺序。34

这不否定胶体实验中实际观察到的团簇路径。它提醒我们:一个体系中能看见明确中间阶段,不等于所有结晶都先变密、再排队。 对照研究时,材料、指标、时间分辨率和空间平均尺度都要对上。

同一根毛细管里,为什么会出现不同结晶结果

2025年研究还通过持续透析,让一根毛细管中的盐逐渐向外扩散。靠近开口的地方先经历变化,较远处改变得慢;于是,同一批颗粒在不同位置经历了不同的相互作用历史。

作者观察到,靠近开口、条件变化较快的区域更容易出现小而不规则的结构和聚集体;改变较缓慢的某些区域能长出较大的规则晶体;另一些位置则出现不同结构共存。2

先把外部驱动这一环单独算清。按论文给定的条件,盐在50毫米长的通道中扩散,初始浓度为4.5 mM,开口处浓度固定为零,封闭端无通量。在横截面均匀、忽略对流的模型中:

$$ \frac{\partial c}{\partial t} =D\frac{\partial^2c}{\partial x^2}, \qquad D=1.61\times10^{-9}\ {\rm m^2/s}. $$

本库独立求解了这一边界问题,没有重新拟合扩散系数。对原表中全部正时间的浓度场值,最大绝对差约为0.019 mM,均方根差约为0.003 mM;这是与作者数值场的比较,不是实验测量误差。

独立扩散计算与公开实验位置。左:不同位置的盐浓度经历,阴影为原文两条参照浓度之间的范围。右:等浓度线与成晶前沿、生长结束位置比较,未平移实验起点。图展示96小时;原数值表延伸到168小时,前述误差比较覆盖全部正时间。原表零时近开口并非理想均匀初场,未纳入比较。等浓度线不是成核动力学预测。
独立扩散计算与公开实验位置。左:不同位置的盐浓度经历,阴影为原文两条参照浓度之间的范围。右:等浓度线与成晶前沿、生长结束位置比较,未平移实验起点。图展示96小时;原数值表延伸到168小时,前述误差比较覆盖全部正时间。原表零时近开口并非理想均匀初场,未纳入比较。等浓度线不是成核动力学预测。

在还没有明显受到封闭端影响时,扩散带来的变化时间尺度随距离平方增长。这使不同位置进入、停留和离开某个作用强度范围的时间不同。

但盐扩散方程并不知道哪里会出现晶核,也不知道晶体有多少缺陷。它只描述了驱动条件怎样变化。成核、颗粒重排、物料供应和表面作用,还要在这个变化的背景上继续发生。

因此,透析实验的价值是把“历史”变成可以计算和比较的条件。要把最终晶体完全归因于某一项历史变量,仍需检查其他同时改变的因素。

一个最小例子:中间态怎样帮助,也怎样拖慢

用一个自选的四状态模型,可以把这两种可能性算清。设无序状态为D,中间态为M,两个可形成的晶体状态为A和B。系统可以从D直接到A或B,也可以先进入M,再返回D或继续到晶体。

转移速率采用常见的越垒形式:

$$ k_{i\to j}=\nu \exp\!\left[-\frac{G^\ddagger_{ij}-G_i}{k_BT}\right]. $$

这里的状态自由能和过渡态位置全部为本库自选,尝试频率设为一个时间单位的倒数;它们不是实验中单对颗粒的势阱深度,也没有从胶体数据拟合。

我们把B设为比A更低的自由能状态,但让中间态更容易接上A的通路。模型在首次到达A或B时结束记时;这示意首次晶体选择,不计算真实空间中的成核过程、后续增长和无限久的平衡。

自选四状态模型的实际概率演化。左为首次到达任一晶体的累计概率,时间使用自选单位。右为首次形成A的最终概率。“更深中间态”只降低M的自由能,保持过渡态位置不变,因而出口速率一起减慢。
自选四状态模型的实际概率演化。左为首次到达任一晶体的累计概率,时间使用自选单位。右为首次形成A的最终概率。“更深中间态”只降低M的自由能,保持过渡态位置不变,因而出口速率一起减慢。
模型条件 平均首次到达晶体状态的时间 第1000个时间单位前已到达的比例 首次到达A的最终概率
只有直接通路 约966 约64.5% 约11.9%
加入可继续转化的中间态 约79.7 接近100% 约97.3%
中间态更深,出口过渡态位置不变 约21,391 约6.34% 约97.3%

第二种情况中,新增通路让系统很容易经过中间态到达A,整体过程更快。

第三种情况中,中间态更稳定,但离开它的代价也变大了。系统可以长时间留在那里,平均等待时间反而远长于没有这条路的时候。

还有一个值得注意的结果:第二和第三种情况的首次晶体选择概率相同,时间却相差很大。因为在这组特定参数变化下,M的各条出口同比减慢,出口之间的相对竞争关系没有改变。

这说明,中间态占了多少、离开它要多久、最后先到哪一个晶体,是不同的问题。 观察到大量团簇,不能独自判断它们是高效前驱态,还是长期滞留的物料。

本库用概率演化、线性方程和独立求积相互核对了该例子。它没有模拟真实空间中的临界晶核、颗粒尺寸或生长过程,数字不能拿来预测胶体实验所需小时数。

初始结构与所在表面,也会改变哪条路容易走

中间态的作用还不只是改变密度。已有的局部排列,可能与某些晶体更容易衔接;表面则可能改变第一层颗粒怎样排列,以及后续结构形成的代价。

在胶体实验中,作者观察到带有贯穿空通道的L3S4结构。它在实验中与表面和透析条件有关;相应模拟比较显示,改变带电表面条件,可以改变L3S4与CsCl-like结构的成核倾向。2

同一批“积木”不等于每一种产物都具有相同组成。L3S4的大小颗粒比例为3:4,而CsCl-like的排列对应另一种比例;比较这些结构时,局部物料组成和表面也必须进入账本。不能把它们全当成相同化学式的多晶型,再只凭每颗粒势能作完整自由能排序。

有些晶体还会沿已有晶体表面长出另一种结构,形成异质外延。作者结合粒子排列与模拟,提出部分界面具有可衔接的几何结构。这展示了另一种历史作用:先形成的结构,能成为后续结构的环境。

然而,不同晶体在一张最终图像里共存,并不能唯一确定它们的成核势垒相近;可用位点、物料组成、成核时间和生长速度也会影响最后看到谁。共存是需要解释的观测,不能替代完整动力学。

从无定形结构里找“晶体种子”,能预测到哪一步

Aykol等人在2025年提出的a²c方法,把这个问题推进到多种无机材料:先生成包含局部短程结构的无定形模型,从中取出许多小区域,再在周期性边界下放松这些区域,寻找能够接续的晶体结构。5

这里的想法是,无定形材料虽然没有长程周期排列,却并非什么结构信息都没有。局部配位、邻居几何和成分关系,可能让一些晶体结构比另一些更容易从它出发到达。

该方法在论文列举的多个体系中,与已有的初始结晶产物吻合。但它不是一段真实成核录像:人为截取并周期重复小区域,再沿能量梯度放松,与真实材料中晶核通过涨落形成、跨过势垒并长大,是不同的过程。

论文与随机结构搜索作比较时,使用了一个叫“加速因子”的量。按其定义,它比较的是找到已知目标后,目标在低能候选结构中的比例,不能解释成晶体长快了几倍,也不能直接当作程序运行时间缩短了几倍。

本库根据公开表重新计算了六个可比较体系的比例:

体系 随机搜索的候选数 a²c候选数 按文中定义可比较的精度比
SiO₂ 24 4 6.00
TiO₂ 25 7 3.57
GeTe 23 10 2.30
BiBO₃ 9 5 1.80
VO₂ 23 13 1.77
Al₂O₃ 15 12 1.25

这些比较采用了论文的已知目标和候选能量窗口,不能把“候选更少”脱离目标是否找到来理解。第七个Mn₂C案例中,论文注明随机搜索在所用设置下没有找到实验目标。附表仍列出候选数21与24及其比值0.875,但这一比值不能按同一精度定义解释;本篇保留原值记录,不将它列为可比较的有限倍数。5a

更重要的科学边界是,无定形前驱态本身也有历史。同样成分的薄膜,制备方式、缺陷、氧含量或厚度不同,初始局部结构就可能不同。作者明确指出,某些TiO₂条件下的板钛矿,以及部分材料转化途中短暂出现的相,需要更具体的制备与转化模拟才能处理。5

因此,这种方法有价值地利用了前驱态的信息,但没有把结晶路径压缩成“给一个化学式,就能唯一预测最后长成什么”。

怎样让证据真正区分路径

如果只有一张最终晶体图像,多种路径可能兼容它。要进一步区分,关键是把结构、时间和物料联系起来。

例如,连续观察浓集团簇内部出现有序区域,与只看到团簇和晶体同时存在,证据强度不同。要区分颗粒先离开团簇、经溶液进入晶体,还是团簇接触后直接被吸收,需要位置追踪或相应的物料流线索;要区别普通聚集与晶粒定向附着,还要核对取向及接触区重排。

测量也会改变可见范围。电子显微镜制样可能改变脆弱中间态,原位成像常只能跟踪有限对象;“这一条路径确实能发生”与“它在全部样品中占主要地位”仍需要不同统计。1

本篇能支持的理解是:中间态通过浓度、局部排列、移动能力和界面条件改变路径,但这些作用没有统一的方向。合适的中间态可以帮助晶体出现,较深或难以重排的状态也可能拖慢它;同一首次晶型选择可以对应很不一样的等待时间。

把这些区别讲清,才能理解为什么相似的起始材料和相似的最终晶体之间,仍然藏着不同的形成机制。尚未解决的,是如何在具体材料中同时测准路径、速率和选择概率,并把有限模型的条件与真实制备历史对上。

可回查:计算结果、实验及模拟状态表、原模拟势能轨迹、结构指标分布、盐浓度场比较、实验前沿位置、自选路径模型、a²c候选比较。

来源与依据

  1. S1 领域综述,概念及观测地图

    Non-classical crystallization in soft and organic materials

    Jingshan S. Du, Yuna Bae, James J. De Yoreo · 2024

    Nature Reviews Materials9,229–248;OSTI作者稿;2024版权更正已核 · 相关全文已核

    Introduction多步/CNT关系、前驱态移动能力及in-situ/ex-situ技术/选择偏差相关段

    全文取得,实际读相关段而非全部通读;路径存在不等于占主导。review-correction2024.html只更正版权行,无科学修订。

  2. S2 原始胶体成像、透析与粒子模拟

    Direct observation and control of non-classical crystallization pathways in binary colloidal systems

    Shihao Zang, Sanjib Paul, Cheuk W. Leung, Michael S. Chen, Theodore Hueckel, Glen M. Hocky, Stefano Sacanna · 2025

    Nature Communications16,3645,2025-04-17;公开全文/SI · 相关全文已核

    Results10–37、Methods64–71;SI图S1–S3的2/3/4页实际渲染读取

    核聚集/成核/生长、密度匹配和bulk/表面对照、透析与结构组成。SI提取乱码,改读落盘PNG。CsCl-like非真实CsCl,L3S4与CsCl粒子配比不同,不能全当同组成多晶型。共存不唯一确定势垒相近;Fig9d第一层能量非全自由能。代码快照879038f2ff1dd55de9c12efcc633edf4c367caa8只读README及define_quantities单位相关函数,未跑MD或通读全部合成方法。

  3. S2DATA s2原实验/模拟数值表,非独立证据

    Source Data for Zang et al. 2025

    Shihao Zang et al. · 2025

    出版MOESM11,2026-09-09取得 · 相关全文已核

    Fig3a状态、6d能量、6e有序指标、5a_color map与5a_data

    39实验/10模拟条件;两不同参数轨迹10000/20000点,direct原4187行时间回跳保留并断线,不编原因。三q6直方和.999/.998/.998各归一。盐场实际纵坐标Y0–50mm,901节点/673时点0–168h;比672正时场,最大差.0186303mM/RMS.00304285mM;原t0开口附近非理想均匀初态而排除。图96h,不平移实验时间/位置起点;cm→mm,误差列保留。未重分割图像,盐场差非实验误差。

  4. S3 原始密度与有序成核模拟

    The microscopic pathway to crystallization in supercooled liquids

    John Russo, Hajime Tanaka · 2012

    Scientific Reports2,505;全文和SI · 相关全文已核

    引言、Voronoi密度/q6/Q6/S定义、硬球路径与多晶型讨论;SI开头分类/方法

    作者强调有序前驱;不把条件分布当一般自由能证明。水/GCM与SI全证明未逐项复核,模拟未跑;和s4的解释分歧保留对象尺度。

  5. S4 原始硬球模拟再分析及竞争解释

    The Early Crystal Nucleation Process in Hard Spheres shows Synchronised Ordering and Densification

    Joshua T. Berryman, Muhammad Anwar, Sven Dorosz, Tanja Schilling · 2016

    arXiv1503.07732v4,2016-04-25;对应2016 JCP论文 · 相关全文已核

    摘要/引言;Conclusion及个体轨迹和集体涨落,提取文本390–425

    所研究尺度密度与有序同步,指标敏感性为作者提出的分歧解释。不宣称推翻所有有序前驱框架;未跑轨迹。机构PDF超时后取得预印本,实际读取版本如实注明。

  6. S5 原始无定形前驱态结构搜索研究

    Predicting emergence of crystals from amorphous precursors with deep learning potentials

    Muratahan Aykol, Amil Merchant, Simon Batzner, Jennifer N. Wei, Ekin Dogus Cubuk · 2025

    Nature Computational Science5,105–111,2025;全文及SI · 相关全文已核

    全主线Results、Fig2、Methods28–42;SI方法入口、Fig3和Notes1/2/4相关段

    周期子区域放松不是实际成核轨迹/势垒/速率。读前驱态、已知产品、精度及边界,未重跑GNN、MQMD或RSS,SI未全读。12例非普适前瞻12/12准确率;TiO2氧含量/厚度/brookite与暂态相需要具体过程。

  7. S5DATA s5候选数表,非独立验证

    Source Data Fig.2 for Aykol et al. 2025

    Muratahan Aykol et al. · 2025

    出版SourceDataFig2全七行 · 相关全文已核

    全表,配s5 Methods precision定义及Fig2图注

    六行目标均找到时精度比=RSS候选数/A2C候选数,实际重算;Mn2C的RSS目标未找到,precision=0,21/24=.875不能解释为有限精度比,保留原值。Fig1CSV已得未分析。不当运行加速、物理成核速度或逐轨迹成功概率。

最后更新 2026-09-09

来源