非球形杂质对非晶态固体延性-脆性转变临界点的影响

背景与学术渊源

起源与学术渊源

非晶态固体力学失效的研究因其在工业和日常生活中的广泛应用而具有至关重要的意义。历史上，尽管晶体固体在外部变形下的力学行为已通过缺陷得到了全面阐释，但非晶态固体缺乏长程结构有序性，这使得精确描绘其力学响应变得困难。这导致了数十年的严谨研究，但仍未对非晶态固体的屈服机制形成全面理解。

一个关键的研究领域是延性-脆性转变，即材料从逐渐变形（延性）转变为灾难性失效（脆性）。近期研究，特别是在非热准静态加载（AQS）条件下，表明这种转变受材料固有的无序强度控制，存在一个临界点来区分延性屈服和脆性屈服。然而，这种临界点的存在在学术界仍存在争议，一些人认为在热力学极限下，对于大体系而言，屈服始终是脆性的。

先前方法的根本局限在于它们无法明确且可实验地控制固有无序性来研究这种转变。例如，退火等方法虽然在分子模拟中有所应用，但存在有限尺寸效应，且不适用于非布朗胶体系统。粒子钉扎，另一种提出的方法，在分子玻璃中仍难以实现。此外，尽管微合金化（添加少量杂质）是提高材料强度的常见工业实践，但其对延性-脆性转变影响的微观机制仍不明确。缺乏一种稳健、实验可行且在微观上得到理解的调控无序性并探测这一临界转变的方法，正是促使作者们进行这项研究的痛点。他们的目标是提出一种使用非球形杂质的新方案来解决这些局限性。

直观领域术语

• 非晶态固体 (Amorphous Solids): 想象一堆随机堆放的乐高积木，而不是整齐地堆叠成一面墙。它是固态的，但没有重复的图案或有组织的结构，不像一个完美搭建的乐高房子（晶体固体）。
• 延性-脆性转变 (Ductile-to-Brittle Transition): 想象一块软黏土与一块干脆饼干。黏土（延性）在断裂前可以显著拉伸和塑形，而脆饼干（脆性）在极小的变形下就会突然断裂。这种转变描述了材料在失效模式上变得更像脆饼干而更不像黏土。
• 非热准静态加载 (Athermal Quasistatic Straining, AQS): 想象在一个温度变化不起作用的房间里，用一只极其轻柔、稳定的手缓慢地推动一块果冻。 "非热"意味着不涉及热量，而 "准静态"意味着过程非常缓慢，系统几乎始终处于平衡状态，使其能够适应变形。
• 旋转自由度 (Rotational Degrees of Freedom, rDoF): 想象在一碗弹珠中加入微小的、细长的珠子。弹珠只能相互滑动，但细长的珠子还可以翻滚和旋转，这给了它们额外的移动和重新排列的方式。这些额外的 "旋转"运动就是它们的旋转自由度。
• 剪切带 (Shear Band): 想象侧向推压一副扑克牌。不是所有的牌都平滑滑动，而是几层牌可能会突然非常快速地相互滑动，形成一个局部化的 "断层线"或带，大部分变形发生在这里，导致突然的、灾难性的失效。

符号表

符号	描述	类型
$\gamma$	剪切应变，变形的度量。	变量
$\sigma_{xy}$	剪切应力，平行于表面的单位面积力。	变量
$\chi_{dis}$	非连通磁化率，应力涨落的度量，指示脆性。	变量
$\Theta$	结构序参量，量化非晶态固体的局部结构有序性或稳定性。	变量
$S_r$	旋转弛豫函数，表征棒状杂质旋转迁移的程度。	变量
$D_{min}^2$	非仿射位移，宏观变形无法解释的粒子重排的度量。	变量
$c_s$	球形杂质的数分数。	参数
$c_d$	二聚体杂质的数分数。	参数
$c_r$	棒状杂质的数分数。	参数
$L_r$	棒状杂质的长度。	参数
$\sigma_{AA}$	母体A型粒子的直径。	参数
$\sigma_{BB}$	母体B型粒子的直径。	参数
$\sigma_s$	球形杂质的直径。	参数
$\sigma_b$	构成二聚体/棒状杂质的珠子的直径。	参数
$T$	温度。	参数
$N$	系统中粒子的总数（系统尺寸）。	参数
$\Delta \sigma_{max}$	最大塑性下降，延性-脆性转变的序参量。	变量
$\chi_d$	$\Delta \sigma_{max}$ 涨落的磁化率，指示临界性。	变量
$c_r^*$	临界棒分数，发生延性-脆性转变的杂质浓度。	参数

问题定义与约束

核心问题表述与困境

本文所解决的核心问题围绕着理解和精确表征非晶态固体的延性-脆性转变。

输入/当前状态： 非晶态固体缺乏长程结构有序性，在外部变形下表现出两种截然不同的失效模式：延性屈服（渐进的材料流动，应力响应连续）和脆性屈服（灾难性失效，通过突然的剪切带形成，应力下降不连续）。已知这种屈服的性质受材料固有无序强度的控制，并且认为在有限的临界无序强度下会发生延性行为和脆性行为之间的转变。微合金化，即添加少量杂质，是提高材料机械性能的常用工程实践，但这些改进对非晶态固体影响的微观机制知之甚少。

输出/目标状态： 本文旨在精确定义和表征一个有限无序临界点，该临界点标志着非晶态固体延性屈服和脆性屈服之间的边界。这是通过引入具有不同旋转自由度（rDoF）的非球形杂质来系统地控制固有无序性来实现的。最终目标是建立一种新颖的、实验可及的协议，用于调控非晶态固体（特别是软玻璃，如胶体系统）的力学性能，以诱导所需的延性或脆性响应。

缺失环节或数学鸿沟：
1. 争议中的临界点： 尽管先前的研究表明存在一个有限无序临界点用于延性-脆性转变，但其存在性和性质，尤其是在热力学极限下，在文献中仍存在争议。先前的方法，如热退火或弹塑性模型，得出了不确定的结果或受到有限尺寸效应的限制。
2. 微合金化的微观机制： 微合金化（特别是使用非球形杂质）影响非晶态固体屈服转变和力学稳定性的精确微观机制在很大程度上是未知的。非球形杂质引入的额外旋转自由度的作用是一个关键的未探索方面。

痛苦的权衡或困境：
核心困境在于杂质特性、旋转自由度与力学响应之间的复杂关系。虽然添加杂质通常可以提高力学强度，但要实现向脆性或延性的显著转变并非易事。球形杂质只能提供最小的屈服应变改善。非球形杂质，特别是棒状杂质，引入了旋转自由度，这可以显著 增强 延性和承载能力。然而，本文揭示了一个反直觉的权衡：减小 这种旋转自由度（通过增加杂质的长径比或人为冻结其旋转）会导致更脆的响应和超稳定般的力学行为。因此，研究人员被困在设计能够通过 rDoF 增强延性的杂质和设计能够通过限制 rDoF 诱导脆性的杂质之间，使得临界转变的精确控制变得具有挑战性。

约束与失效模式

作者们遇到了几个严峻的现实障碍，使得这个问题极其难以解决：

• 物理/结构约束：
  • 缺乏长程有序性： 与晶体固体中可以精确描绘缺陷不同，非晶态固体缺乏长程结构有序性，这使得定义和追踪缺陷本身就具有内在困难，从而使理解其力学响应变得复杂。
  • 标准退火方法不适用： 调控无序性的传统方法，如热退火或气相沉积，对于所有系统，特别是对于非布朗胶体玻璃，效果不佳或不适用。这限制了这些系统中延性-脆性转变的实验可及性。
• 计算约束：
  • 大尺寸杂质的缓慢动力学： 包含非常大的球形杂质或大直径二聚体杂质的系统表现出极慢的动力学。这使得在可行模拟时间尺度内使这些系统达到平衡在计算上具有挑战性，限制了可研究的杂质尺寸范围。
  • 棒长模拟限制： 存在一个实际的最大棒长 ($L_r = 2.5\sigma_{AA}$) 限制。这阻止了直接观察极长棒的完全旋转冻结，而极长棒自然具有消失的旋转扩散性。
  • 长棒系统退火困难： 具有较长棒的系统动力学更慢，这意味着标准的制备方案产生的退火状态较差。这种模拟困难可能导致对这些系统中观察到的脆性估计不足。
  • 自由旋转自由度计算成本高昂： 在允许其完全旋转自由度的情况下模拟具有非常大长径比杂质的系统在计算上非常密集。为规避此问题，作者们采用人为冻结 rDoF 的方法来模拟非常长、旋转不动的棒状杂质的影响。
• 数据驱动/分析约束：
  • 微观机制难以捉摸： 尽管有实验证据表明微合金化可以提高屈服应变，但其精确的微观机制仍然难以捉摸，需要详细的计算分析来揭示。
  • 参数范围有限： 研究主要集中在长径比和旋转自由度上。其他关键参数，如粒子形状（例如，椭圆形与棒状）和边界相互作用的性质（例如，构成棒状的球体重叠程度），已被承认但由于复杂性而留待未来研究。
  • 非热准静态假设： 当前分析是在非热准静态加载（AQS）条件下进行的。作者们承认，包含热效应和有限应变率的结果可能会发生变化，而这些在本文中并未得到充分探讨。

为什么选择这种方法

选择的必然性

作者选择采用分子动力学（MD）模拟，并引入具有冻结旋转自由度（rDoF）的非球形、棒状杂质，这并非仅仅是偏好，而是由现有方法的局限性和问题的特定要求所驱动的战略性必需。核心挑战在于研究非晶态固体的延性-脆性转变，并利用一种既实验可及又能精确控制固有无序性的方案来确立有限无序临界点的存在。

传统的“最先进”（SOTA）方法，如热退火，因多种原因被证明不足。虽然热退火可以改变样品的无序强度，但先前的研究 [15] 仅提供了“由于有限尺寸效应而得到的部分支持”，并且大规模弹塑性模型 [22] 的结论仍然“不确定”。更关键的是，热退火“不适用于”非布朗胶体系统，而胶体系统是研究玻璃转变的一个关键实验平台。同样，虽然粒子钉扎 [28] 提供了控制无序性的另一种途径，并且可以在胶体玻璃实验 [29] 中实现，但它“在分子玻璃中仍然具有挑战性”。作者们明确指出，“识别新的协议来研究这种转变，特别是那些实验可及的协议，可能对该领域产生重大影响。”

引入非球形杂质，特别是棒状杂质，直接解决了这一差距。通过操纵其长径比和旋转自由度，作者们发现了一种调控固有无序强度的新方法，这种方法“在分子玻璃和胶体玻璃实验中都易于实现”。此外，模拟“具有消失的旋转扩散性的非常长的棒”由于其缓慢的动力学而带来的计算困难，使得在能量释放过程中“手动冻结”rDoF成为必要。这种模拟技术使他们能够模拟这些长棒的影响，而这些影响在可用的模拟时间尺度内是难以处理的，从而使这种特定方法成为探索所需区域的唯一可行途径。

比较优势

所选择的包含非球形杂质（特别是具有可控旋转自由度的棒状杂质）的方法，与包括球形杂质和热退火在内的先前黄金标准相比，表现出压倒性的定性优势。

在比较球形杂质和非球形（二聚体）杂质时，本文强调了球形夹杂物在力学性能方面“最小”的改善。例如，对于 3DKA 系统，在 $c_s = 0.1$ 时，屈服应变 ($\gamma_y$) 仅从纯系统的 $\gamma_y = 0.09$ 增加到 $\gamma_y = 0.107$（图 2(a,b)）。相比之下，在 $c_d = 0.1$ 时，非球形二聚体杂质使屈服应变增加了“近 100%”，达到 $\gamma_y = 0.127$（图 3(a,b)）。这是一个显著的定性差异。

结构优势在于非球形粒子引入的“额外的旋转自由度（rDoF）”。这些 rDoF 提供了“额外的局部路径来耗散内部应力”，使系统能够“在非局部剪切带事件导致失效之前，承受更高的载荷并储存更多应力”。这种机制使得系统能够承受显著更大的载荷，并比缺乏这些旋转自由度的球形杂质更有效地延迟屈服转变。

此外，在模拟中 冻结 这些 rDoF 的能力，模拟非常长的棒状杂质，会导致“超稳定般的力学响应”和“极度脆性”的失效（图 5(e,f)）。这证明了一种强大的结构控制机制：rDoF 的存在增强了延性，而它们的限制或缺失则促进了脆性和超稳定性。结构序参量 $\Theta$ 随着杂质分数和棒长增加而持续降低（图 2(d), 图 4(e)），表明结构稳定性增强，这比控制较差的方法具有直接的定性优势。通过简单地调整杂质的长径比和旋转自由度，该方法提供了对材料力学响应的精确、可调控的控制，从延性到高度脆性。

与约束的对齐

所选方法与问题定义中概述的约束完美对齐，清晰地展示了“严苛要求与解决方案独特属性的结合”。

一个主要约束是需要一种“在更广泛的实验系统范围内可及的”协议来研究延性-脆性转变，包括“对于热退火不适用的非布朗胶体系统”。使用非球形杂质，特别是棒状杂质，直接解决了这个问题。作者们明确指出，他们的方法提供了一种“在分子玻璃和胶体玻璃实验中都易于实现的”新颖延性-脆性转变。这使其成为处理软玻璃的实验人员的一种高度相关且实用的方法，在软玻璃中，传统退火方法不可行。

另一个隐含的约束是需要一种能够精确控制“固有无序强度”以探测临界点的方法。通过改变非球形杂质的数量分数和长径比，以及控制它们的旋转自由度，作者们获得了对系统结构稳定性和力学响应的精细控制。结构序参量 $\Theta$ 作为这种受控无序性的定量度量，显示出与观察到的力学行为的明确相关性。

最后，模拟非常长的棒状杂质（表现出“非常缓慢的动力学”和“消失的旋转扩散性”）的计算约束，通过在应力释放过程中人为冻结 rDoF 来克服。这种巧妙的变通方法使作者们能够“模拟非常长的棒状杂质的影响”，并研究由此产生的超稳定、脆性行为，而这在计算障碍下本不可能实现。这表明了在存在计算障碍的情况下，研究现实物理现象的需求与实际解决方案的务实对齐。

替代方案的拒绝

本文基于其在解决非晶态固体延性-脆性转变问题（尤其是在对各种系统进行实验可及性方面）的局限性，隐含且明确地拒绝了几种替代方法。

1. 热退火： 作者们指出，虽然热退火可以改变无序强度，但它“由于有限尺寸效应而提供部分支持”，并且其结论仍然“不确定”[15, 22]。更重要的是，它“不适用于”非布朗胶体系统，而胶体系统是实验验证的关键目标。即使对于分子玻璃，通过退火获得“现实系统”也很困难，因为“即使是可及的最低冷却速率也会导致退火状态不佳”。像交换蒙特卡洛 [34, 35] 这样的方法可以产生超稳定玻璃，但“不适用于实验”。这使得热退火成为一种不合适的通用方法。

2. 粒子钉扎： 尽管粒子钉扎在最近的研究 [28] 中被用于控制胶体玻璃实验 [29] 中的无序性，但作者们表示它“在分子玻璃中仍然具有挑战性”。虽然这是一种有前途的技术，但它不像他们提出的微合金化策略那样，在不同材料类型之间提供相同的广泛适用性。

3. 球形杂质： 本文直接将其非球形杂质方法与球形杂质进行了比较。虽然球形杂质确实导致了“更高的屈服应变和脆性增加”，但作者们得出结论，“改善是最小的”（例如，对于 3DKA，$\gamma_y$ 从 0.09 增加到 0.107）。这是一个明确的、数据驱动的拒绝，认为球形杂质不足以有效控制延性-脆性转变。球形粒子缺乏旋转自由度限制了它们与非球形粒子相比耗散应力和增强稳定性的能力。

本文没有讨论机器学习方法，如 GANs 或扩散模型，因为其重点是材料性质的物理模拟，而不是生成模型。拒绝替代方案的根源在于它们无法为所研究的特定现象提供足够的控制、实验可及性或所需程度的力学性能增强。

FIG. 6. Ultra-stability with frozen rDoF: Kinetic as- pect. (a) The potential energy per particle e is plotted for a system with 10% dimers subjected to heating-cooling cycles at the rate of dT/dt = 10−4 (same as the rate with which the sample was prepared). The solid lines are for a system with frozen rotational DoFs; the observed hysteresis during the first cycle indicates the ultra-stable character. The sec- ond heating cycle does not show any hysteresis. The dashed lines are for the system evolved with rotational DoFs, and the absence of hysteresis is seen as expected. (b) The specific heat CV = de/dT is calculated from data in (a) by numerical differentiation and conveys the same

数学与逻辑机制

主方程

本研究中控制粒子行为以及非晶态固体力学性能的基本相互作用由 Lennard-Jones (L-J) 势描述。这种势能函数决定了系统中任意两个粒子如何相互作用，构成了所执行的分子动力学模拟的基础。它是系统弛豫到其固有状态过程中最小化的目标函数。

$$V_{\alpha\beta}(r) = 4\epsilon_{\alpha\beta} \left[ \left(\frac{\sigma_{\alpha\beta}}{r}\right)^{12} - \left(\frac{\sigma_{\alpha\beta}}{r}\right)^6 \right]$$

逐项解剖

让我们来剖析这个方程，以理解每个组成部分的作用：

• $V_{\alpha\beta}(r)$:
  1) 数学定义： 该项表示类型为 $\alpha$ 和 $\beta$ 的两个粒子之间，距离为 $r$ 时的相互作用势能。
  2) 物理/逻辑作用： 它量化了粒子对之间相互作用中储存的能量。这种能量决定了它们之间的力，进而决定了它们的运动、排列以及材料整体的结构和力学响应。
  3) 为何是减法： 方括号内的两项代表了不同的物理现象：排斥和吸引。减法表示这些是相反的力。第一项（排斥项）将粒子推开，而第二项（吸引项）将它们拉拢。

• $4\epsilon_{\alpha\beta}$:
  1) 数学定义： 相互作用能的缩放因子，其中 $\epsilon_{\alpha\beta}$ 是 $\alpha$ 和 $\beta$ 类型粒子之间相互作用的势阱深度。
  2) 物理/逻辑作用： 这个系数设定了相互作用的强度或幅度。较大的 $\epsilon_{\alpha\beta}$ 意味着更强的吸引力和排斥力，导致系统更稳健或结合更紧密。它定义了相互作用的能量尺度。
  3) 为何是乘法： 它作为整个势函数的全局乘数，确保了吸引和排斥分量整体强度的比例缩放。

• $\sigma_{\alpha\beta}$:
  1) 数学定义： $\alpha$ 和 $\beta$ 类型粒子之间势能为零的距离。它有效地代表了粒子的“尺寸”或碰撞直径。
  2) 物理/逻辑作用： 这个参数定义了相互作用的长度尺度。它决定了粒子的有效“硬核”尺寸，影响了它们的堆积密度以及吸引力显著的范围。它对于建立材料的密度和局部结构至关重要。
  3) 为何除以 $r$： 它通过特征粒子尺寸 $\sigma_{\alpha\beta}$ 来归一化粒子间距离 $r$，使得括号内的项无量纲，并允许势函数以相对分离的形式表示。

• $r$:
  1) 数学定义： 两个相互作用粒子中心之间的瞬时距离。
  2) 物理/逻辑作用： 这是基本空间变量。粒子间的势能，以及由此产生的力，直接依赖于这个分离距离。
  3) 为何是幂律（$r^{-12}$ 和 $r^{-6}$）： 这些幂律选择用于模拟原子间和分子间力的特征性质。$r^{-12}$ 项代表由于电子云重叠而产生的非常短程、强烈的排斥力，而 $r^{-6}$ 项代表较长程、较弱的吸引力（例如，范德华力）。

• $\left(\frac{\sigma_{\alpha\beta}}{r}\right)^{12}$:
  1) 数学定义： Lennard-Jones 势的排斥分量，与 $r^{-12}$ 成正比。
  2) 物理/逻辑作用： 该项模拟了强烈的、短程的排斥力，防止粒子重叠或占据相同的物理空间。它确保了粒子的“硬核”性质。
  3) 为何指数是 12： 这个陡峭的幂律是一个经验选择，它有效地近似了粒子非常接近时排斥力的快速增加，模拟了泡利不相容原理。

• $\left(\frac{\sigma_{\alpha\beta}}{r}\right)^{6}$:
  1) 数学定义： Lennard-Jones 势的吸引分量，与 $r^{-6}$ 成正比。
  2) 物理/逻辑作用： 该项模拟了较长程的吸引力（例如，范德华力），这些力负责将材料结合在一起。
  3) 为何指数是 6： 这个较缓和的幂律是诱导偶极-偶极相互作用的特征，这种相互作用在许多原子和分子系统中普遍存在。

逐步流程

Lennard-Jones 势是用于探索延性-脆性转变的分子动力学模拟的驱动引擎。以下是抽象数据点（粒子的位置）如何通过模拟的流程：

1. 系统设置： 想象一组粒子，每个粒子由其类型（$\alpha$ 或 $\beta$）和初始空间坐标 $(X_i, Y_i, Z_i)$ 定义。这些类型决定了其相互作用的具体 $\epsilon_{\alpha\beta}$ 和 $\sigma_{\alpha\beta}$ 参数。
2. 高温平衡： 系统被“加热”到高温。每个粒子的位置根据其受到的净力不断更新，该净力源自总势能（所有与其邻居的 $V_{\alpha\beta}(r)$ 相互作用之和）的梯度。粒子剧烈运动，探索能量景观，直到达到稳定的液态。
3. 冷却与淬火： 然后通过逐渐降低粒子的动能来“冷却”系统。随着温度降低，来自 $V_{\alpha\beta}(r)$ 的吸引力变得更占主导地位，导致粒子沉淀到无序的非晶态固体状态。最后，通过最小化其总势能来将系统“淬火”到其固有状态。这意味着粒子被移动到使所有 $V_{\alpha\beta}(r)$ 之和处于局部最小值的位置，从而有效地消除所有动能并达到力学稳定的构型。
4. 非热准静态加载（AQS）循环：
   • 仿射变换： 将一个小的增量剪切应变 $\delta\gamma$ 应用于整个系统。对于位于 $(X_i, Y_i, Z_i)$ 的粒子，其 x 坐标被仿射变换为 $X_i \leftarrow X_i + \delta\gamma Y_i$。此步骤会扭曲粒子相对于其邻居的位置，暂时使系统脱离力学平衡。
   • 能量最小化（弛豫）： 在仿射步骤之后，系统不再处于能量最小值。然后允许粒子弛豫。这包括迭代地调整它们的位置以最小化总势能，使用共轭梯度法。力（$V_{\alpha\beta}(r)$ 的梯度）引导每个粒子朝着新的局部能量最小值移动。这个弛豫过程模拟了材料在应力下的变形和重排。
5. 迭代变形： 步骤 4a 和 4b 重复多次小的应变增量。在每次弛豫步骤后记录系统的响应（例如，应力）。这个迭代过程允许研究人员追踪材料的应力-应变曲线，并观察塑性事件和屈服行为。

优化动力学

这里的“优化”指的是系统主要通过能量最小化来寻找力学稳定构型的过程。

• 能量景观： 系统的总势能 $U = \sum_{\text{所有对 } \alpha\beta} V_{\alpha\beta}(r)$ 定义了一个复杂的高维能量景观。该景观的特点是存在许多局部最小值（代表稳定构型或“固有状态”），它们被能量势垒（鞍点）分隔开。
• 梯度作为驱动力： “学习”或“更新”机制依赖于作用在每个粒子上的力，这些力是总势能相对于其位置的负梯度（$\mathbf{F}_i = -\nabla_i U$）。这些力决定了弛豫过程中粒子运动的方向和大小。
• 共轭梯度法： 该迭代算法用于导航能量景观并寻找局部最小值。其操作方式是：
  1. 计算梯度： 在当前的粒子位置，计算力（梯度）。
  2. 确定搜索方向： 选择一个与先前方向“共轭”的搜索方向，旨在高效地下降到最小值。这有助于避免之字形运动并加速收敛。
  3. 线搜索： 粒子沿着选定的方向移动，以在该特定方向上找到能量最小点。
  4. 更新位置： 更新粒子位置，然后重复该过程。
• 状态更新与收敛： 系统的状态（粒子位置）被迭代更新，直到所有粒子的净力低于预定义的容差。此时，系统已收敛到局部能量最小值，表明力学平衡。这个迭代过程允许系统在受到仿射应变扰动后“弛豫”并找到稳定构型。
• 塑性与屈服： 随着应变的累积，系统最终会达到其占据的局部最小值变得不稳定的点。为了弛豫，它必须克服能量势垒并过渡到新的、不同的局部最小值。这些过渡是“塑性事件”，对应于突然的、不可逆的粒子重排。这些事件的集体行为塑造了材料的屈服响应。在延性屈服中，塑性事件数量众多且空间分布，导致应力响应渐进。在脆性屈服中，塑性事件局部化为剪切带，导致灾难性的、不连续的应力下降。能量景观的形状，由 L-J 势决定，控制着这些重排的势垒和路径，从而控制了材料的力学行为。杂质的引入，特别是非球形杂质，改变了这种景观，改变了这些塑性事件的难易程度和性质。

结果、局限性与结论

实验设计与基线

作者们通过对 Lennard-Jones (L-J) 粒子进行分子动力学模拟，精心设计了他们的实验，主要关注二维和三维的 Kob-Andersen (KA) 模型玻璃形成体。核心思想是将各种类型的杂质引入二元玻璃基体中，并在非热准静态加载（AQS）条件下观察它们对屈服转变的影响。

实验设置包括通过在高温下平衡系统来制备静态状态，然后冷却和淬火到固有状态。之后，通过增量方式施加变形。与提出的机制进行测试的“受害者”或基线模型包括：
1. 纯系统： 没有杂质的玻璃系统（$c_s = 0.00$, $c_d = 0.00$, $c_r = 0.00$）。这些作为延性行为的基本参考。
2. 带球形杂质的系统： 添加了第三种粒子，其直径较大（$\sigma_s = 2.0 \sigma_{AA}$），在 $N_T = 100000$ 个粒子的系统中改变其数量分数（$c_s \in [0, 0.1]$）。
3. 带非球形二聚体杂质的系统： 引入了由两个重叠球形珠子组成的棒状粒子（直径 $\sigma_b = 2.0 \sigma_{AA}$，长径比 $L_d : \sigma_b = 2.3 : 2$），数量分数（$c_d$）不同。这些被设计成具有旋转自由度（rDoF）。
4. 带更长棒状杂质的系统： 在固定分数（$c_r = 0.1$）下添加了由多个珠子组成的更非球形的棒状杂质（直径 $\sigma_r = \sigma_{AA}$），并改变其长度（$L_r$）以控制旋转自由度。
5. 带人为冻结 rDoF 的系统： 为了分离旋转的影响，在能量最小化步骤中人为地阻止了棒状杂质的旋转，同时仍然允许仿射变换。这是一个关键的对照实验。
6. 退火状态不佳的样品： 对于延性-脆性转变研究，样品在相对较高的温度（$T=0.6$）下制备，并快速冷却（$T=10^{-1}$）以确保退火状态不佳、固有的延性状态。

力学响应通过几个关键可观测量进行了严格表征：
* 应力-应变曲线（$\sigma_{xy}$ vs. $\gamma$）： 直接观察屈服应变、剪切模量以及应力超调或不连续下降的存在。
* 非连通磁化率（$\chi_{dis}(\gamma)$）： 定义为 $\chi_{dis}(\gamma) = N_T (\langle \sigma_{xy}^2 \rangle - \langle \sigma_{xy} \rangle^2)$，其峰值幅度和尖锐度指示了屈服转变的脆性。更大、更窄的峰值表示更脆的失效。
* 结构序参量（$\Theta$）： 一种局部度量（公式 (4)），用于量化静态状态的结构稳定性与无序性。较低的 $\Theta$ 表示更好的稳定性。
* 旋转去相关函数（$S_r(\gamma)$）： 用于量化棒状杂质在力学加载下的旋转弛豫程度。更快的衰减意味着更好的独立迁移性。
* 非仿射位移（$D_{min}^2$）： 用于可视化和量化塑性事件的空间分布和剪切带的形成。
* 单位粒子势能（$e$）和比热（$C_v = de/dT$）： 在加热-冷却循环中使用，以评估玻璃的动力学稳定性和超稳定特性。

最后，为了无情地证明其数学关于临界点的论断，作者们采用了有限尺寸标度（FSS）分析。他们改变了粒子总数（$N_T = [25000, 50000, 100000, 200000]$）对于具有旋转冻结棒的系统，分析了 $\chi_{dis}$ 峰值高度、$\langle \Delta \sigma_{max} \rangle$ 及其涨落（$\chi_d$）随系统尺寸的标度行为。这允许外推到热力学极限（$N \to \infty$）并与随机场伊辛模型（RFIM）普适类进行比较。

证据证明了什么

本文提供的证据为非球形杂质和旋转自由度在调控非晶态固体力学性能和延性-脆性转变中的作用提供了一个引人入胜的叙述。

1. 球形杂质增强脆性和稳定性：

   • 应力-应变曲线（图 2a, b）清楚地表明，添加球形杂质显著增加了屈服应变和剪切模量。例如，3DKA 模型表现出更明显的应力超调。
   • $\chi_{dis}$ 图（图 2c）显示，随着球形杂质分数（$c_s$）的增加，峰值移向更高的应变并变得明显更尖锐，这是向更脆性失效转变的明确证据。
   • 结构序参量 $\Theta$（图 2d）随着 $c_s$ 的增加而减小，表明这些杂质导致玻璃基体结构稳定性增强和有序性提高。这种效果比简单降低冷却速率所达到的效果更显著。

2. 具有 rDoF 的非球形二聚体提供卓越的力学增强：

   • 引入轻微非球形的二聚体杂质（图 3a, b）与球形杂质相比，在屈服应变和剪切模量方面产生了更显著的增加。对于 3DKA 模型，$c_d = 0.1$ 的二聚体分数导致屈服应变增加了 40%，而球形杂质仅增加了 18%。
   • 应力贡献的微观分析（图 3d）提供了不可否认的证据，表明棒-球相互作用在更高的应变下继续维持应力，而球-球相互作用则较早饱和。这表明含有棒状结构的区域在结构上更稳定，并且可以承受更高的载荷。
   • 这种卓越的增强直接归因于非球形粒子的额外旋转自由度（rDoF）。这些 rDoF 提供了耗散内部应力的局部路径，使系统能够在灾难性失效之前承受更高的载荷并储存更多应力。

3. 限制 rDoF（更长的棒状杂质）导致脆性：

   • 随着棒状杂质长径比的增加（即更长的棒状杂质），它们的旋转自由度减小，如旋转去相关函数 $S_r(\gamma)$ 的更快衰减所示（图 4b）。
   • rDoF 的这种减少导致屈服点向较低值移动，并且力学响应变得越来越脆性（图 4a）。
   • $\chi_{dis}$ 峰值（图 4d）随着棒长的增加而变得更窄、更大，证实了脆性相的出现。
   • 至关重要的是，非仿射位移图（$D_{min}^2$）（图 4f, g, h）显示了一个清晰的转变：较短的棒状杂质允许塑性事件在空间上分散（延性行为），而较长的棒状杂质则导致形成局部化的、贯穿整个系统的剪切带（脆性失效）。这是核心机制（rDoF）直接影响失效模式的有力证据。

4. 人为冻结的 rDoF 诱导超稳定性和极端脆性：

   • 最令人信服的证据来自人为冻结杂质的 rDoF。这种操作极大地改变了力学响应，导致了具有不连续应力下降的极端脆性失效（图 5a, c, e）。具有自由 rDoF 的系统观察到的显著屈服应变增强完全消失，再次证实了旋转的关键作用。
   • $\chi_{dis}$ 峰值（图 5b, d, f）变得异常尖锐和高，表明了高度脆性的响应。
   • 动力学稳定性分析（图 6a, b）显示，对于具有冻结 rDoF 的系统，在加热-冷却循环中势能存在明显的滞后现象，这是超稳定玻璃的标志。相比之下，具有自由 rDoF 的系统则没有这种滞后现象。这证明了冻结 rDoF 同时导致了力学和动力学超稳定性。

5. 有限无序临界点控制延性-脆性转变：

   • 通过改变退火状态不佳样品中旋转冻结棒的比例（$c_r$），作者们展示了一个从延性到脆性屈服的系统性转变（图 7a, b）。
   • 有限尺寸标度分析为热力学极限下的有限无序临界点提供了明确、不可否认的证据：
     • 对于掺杂系统，$\chi_{dis}$ 峰值幅度随系统尺寸呈幂律 $N^{1.1 \pm 0.03}$ 增长，而对于纯系统，它保持较小且无响应（图 8b, c）。
     • 对于延性系统（$c_r=0$），最大塑性下降 $\langle \Delta \sigma_{max} \rangle$ 随 $N^{-0.4 \pm 0.02}$ 消失，而对于脆性系统（$c_r=0.1$），它保持有限并随 $N^{0.2 \pm 0.04}$ 增长（图 8d, e）。
     • $\langle \Delta \sigma_{max} \rangle$ 的涨落，记为 $\chi_d$，呈现非单调峰值，该峰值随系统尺寸变尖锐和增大，呈 $N^{0.36 \pm 0.027}$ 发散，而其半高全宽（FWHM）随 $N^{-0.46 \pm 0.036}$ 消失（图 8f, g）。
     • 使用这些幂律的曲线拟合（图 8h）证实了转变的临界性质。
     • 将 $c_r^*(N)$ 外推到热力学极限（$N \to \infty$）得到一个有限的临界棒分数为 $c_r^* \sim 0.045 \pm 0.0003$（图 9b）。这一结果反驳了在热力学极限下不存在延性-脆性转变的观点，为真正的无序驱动的临界转变提供了可靠证据。

局限性与未来方向

尽管本文提出了一个出色且彻底的分析，但与任何科学探索一样，它也在一定的约束下进行，并开启了许多未来探索的途径。

一个显著的局限性是对较高浓度旋转受限杂质下微观机制的不完全理解。作者们承认，剪切诱导的向列序对力学响应的影响，特别是在较高杂质浓度下，仍然不清楚。尽管作者们表示，对于各向同性的初始取向和低旋转扩散或冻结 rDoF 的系统，这些影响似乎很小，但对这一方面的深入研究可能会揭示更细微的相互作用。

另一个约束是主要关注非热准静态加载（AQS）条件。作者们自己指出，热效应和有限应变率可能会改变从延性到脆性的临界转变，特别是对于旋转冻结的棒状杂质。尽管他们预计即使在有限应变率和温度下，微合金化系统也会表现出增强的力学强度，但包含这些变量的全面研究对于获得完整图景至关重要。

此外，研究强调了粒子形状（椭圆形与棒状）和边界相互作用的性质（例如，构成棒状的球体重叠程度）是控制力学响应的关键参数，但这些影响留待未来研究。这表明当前模型存在简化，可以进行扩展。计算限制也阻止了模拟具有消失旋转扩散性的非常长的棒状杂质，这使得采用“冻结 rDoF”方法作为代理成为必要。虽然有效，但直接模拟这些系统可以提供额外的见解。

展望未来，本文的研究结果为讨论和研究方向提供了令人兴奋的起点：

1. 探索杂质设计的完整参数空间： 未来的工作可以系统地研究各种杂质形状（例如，椭球体、分支聚合物）和与主体基体不同的相互作用势的影响。这将提供对杂质特性如何影响力学性能和延性-脆性转变的更全面理解。
2. 整合热效应和应变率： 将当前的 AQS 分析扩展到有限温度和应变率至关重要。这将使研究人员能够绘制出在更现实条件下的延性-脆性转变相图，可能揭示新的临界点或交叉点。将热退火方法（如交换蒙特卡洛）与棒状杂质掺杂相结合，可以为固有无序性、杂质诱导的无序性以及动力学效应之间的相互作用提供宝贵的见解。
3. 在胶体和分子系统中的实验验证： 本文强调了这种新颖延性-脆性转变的实验可及性，特别是在胶体系统中。未来的研究应侧重于使用先进成像技术直接进行实验验证，以观察塑性事件、剪切带形成，并在掺有非球形杂质的胶体玻璃中测量力学性能。这将弥合模拟与现实世界材料之间的差距。对于传统退火困难的软玻璃而言，这种微合金化策略可能是一个突破。
4. 超稳定性和剪切带的微观机制： 对于较高浓度旋转受限杂质下驱动转变的微观机制，有必要进行更深入的研究。这包括理解剪切诱导的向列序如何影响力学响应，以及长棒施加的静态关联长度与过冷区域超稳定玻璃的点-设定（PTS）关联长度的比较。这可能涉及先进的理论建模和高分辨率模拟。
5. 在材料设计中的应用： 研究结果提供了一种简单、可控的微合金化策略，用于调控非晶态固体的力学性能。未来的工作可以探索如何将这些原理应用于设计具有定制延性或脆性的新材料，以满足特定的工业应用需求，从坚固的复合材料到自修复材料。通过冻结 rDoF 诱导超稳定性和脆性失效的能力，可能在需要高力学稳定性的应用中特别相关。
6. 普适类探索： 进一步详细的有限尺寸标度分析，可能使用更多系统尺寸和不同的长径比，可以巩固与 RFIM 普适类的联系，并探索在不同条件下是否会出现其他普适类。这将加深我们对无序系统中非平衡相变的基本理解。

FIG. 1. Components of the system. (a, b) Particle type ‘A’ and ‘B’ of the parent KA system. (c) Spherical impurity with twice the diameter of the particle ‘A’. It has a variable number fraction of cs in the system. (d) Slightly aspherical dimers composed of larger particles with σb = 2.0, added at number fraction cd in the system. (e) Rods with larger aspect ratio, made by attaching A-type particles, studied across dif- ferent aspect ratios at a fixed number fraction of cr = 0.1 FIG. 4. Mechanical properties with aspherical impurities: Effect of the length of rod impurity. (a) Stress-strain curves for 3dKA systems doped with 10% rods (Fig. 1(e)) of different lengths, Lr. The yield point shifts back as Lr increases, while the response becomes increasingly brittle due to restricted rotations. (b) Rotational de-correlation function for rods of different lengths with mechanical loading; the red points indicate the equal net D2 min . (c) Single ensemble stress-strain plots for systems with Lr = 2.5 (Lr = 1.3) show abrupt-brittle (continuous-ductile) yielding. (d) Corresponding susceptibility plots with increasing magnitude and sharpness indicate the sample’s emerging brittle behavior. (e) The structural order parameter decreases with increasing Lr, indicating greater structural stability. (f) D2 min, averaged in spatial-strips perpendicular to the shear band (at dx = 0) for systems with rods of varying Lr. Strain values are selected to yield a similar area under the curve, ensuring equal plasticity. Larger spread for smaller Lr and large displacement peak for larger Lr advocates the emerging brittle behavior. (g, h) D2 min maps obtained at equal net displacements for Lr = 1.3 and Lr = 2.5, respectively, illustrating the contrast between ductile and brittle responses FIG. 8. Impurity driven ductile-to-brittle transition, FSS study. (a) System size dependence of stress-strain curves for the 3dKA system prepared at a high temperature of T = 0.6. The pure system, shown in lighter colors, exhibits a ductile macroscopic response with no system size dependence, while systems with 10% rod impurities (aspect ratio 1.9 : 1) show progressively brittle behavior as system size increases. (b) Susceptibility plots for both the pure and doped systems at various system sizes. (c) System-size dependence of peak height for different impurity concentrations. For cr ≤0.02, the χp dis saturates, while it follows the indicated power law for cr ≥0.05. (d) System size dependence of ⟨σmax⟩with increasing cr. The response shifts from vanishing to not-vanishing in the thermodynamic limit on increasing cr. (e) The increase in ⟨σmax⟩with increasing cr becomes sharper with increasing system size. (f) Ensemble-level fluctuations of ⟨σmax⟩reveal a nonmonotonic peak structure, marking the transition point. As system size increases, the peak narrows, and its height increases, highlighting the critical nature of the transition. (g) The peak height diverges as χp d ∼N 0.36±0.027, while the FWHM vanishes as FWHM ∼N −0.46±0.036. (h) The obtained collapse of data points in panel (f) using the determined power laws

与其他领域的同构性

结构骨架

本文的核心机制描述了嵌入式各向异性元素的局部动力学自由度如何决定无序基体中的全局临界转变行为。