以下内容是对文献 Phys. Rev. Lett. 131, 108301 (2023) - Extreme Spontaneous Deformations of Active Crystals (aps.org) 的机翻 + 简单修改。

翻译可能存在一些不足之处，还望见谅。

摘要

我们展示了由活性粒子构成的二维晶体可以在不熔化的情况下经历极其大的自发变形。通过使用主要通过成对斥力相互作用的粒子，我们发现这种活性晶体保持长程 (long-range) 键序 (bond order) 和代数衰减的 (algebraically decaying) 位置序 (positional order)，但指数 不受 (平衡) KTHNY 理论给出的 的上界。我们使用线性弹性理论 (linear elastic theory) 解释了我们的发现，并展示了量化大尺度变形和键序波动 (bond-order fluctuations) 的两个明确定义的有效温度的存在。这些现象的根源在于粒子固有轴的唯一时间持久性，因此应该在许多不同的情况下观察到它们。

正文

二维晶相通常由位置序的缓慢代数衰变和长程键序的存在来定义。在热平衡中，位置相关函数 (positional correlation function) 的衰变指数 随温度而增加。在许多情况下，著名的 Kosterlitz-Thouless Halperin-Nelson Young (KTHNY) 理论 1-2 很好地描述了，二维晶体的熔化是一个两步过程，在晶体和液相之间产生具有短程位置序和准长程键序的中间相——通常是六角相 (hexatic phase) (参见，例如 3，了解经典综述)。第一步被定义为当热激活的位错对可以解离时。这通常发生在孤立位错的熵和弹性能平衡时。KTHNY 理论告诉我们，当随着温度升高的 达到 到 之间的值时，(对于六方晶体) 这一定会发生。这种结合可以被视为对二维晶体的可变形性施加了限制。

由活性粒子构成的晶体通过在局部时间尺度和长度尺度上持续注入机械功来维持非平衡状态，因此不受与其平衡对应物相同的约束。特别是基于这样的一般性原因，人们不应该期望 KTHNY 理论适用于二维活性晶体的熔化。因此，迄今为止，大多数对这类晶体 (不具有手性或定向序) 熔化的简化版本的研究都得出了 KTHNY 理论仍然适用的结论，或者假设它仍然适用于决定熔化何时以及如何发生。4-5

在这篇文章中，我们展示了二维活性晶体可以在不熔化的情况下经历极端的自发变形：它们可以保持真正的长程键序，即使位置序相关性衰减得非常快，也可以抵抗位错对的解除结合，尽管仍然是代数衰减。我们认为衰变指数 实际上可以任意大，并通过两个明确定义的有效温度来解释我们的发现，这两个温度分别量化了弹性变形和键序波动。

许多研究中的活性晶体由旋转单元或个体自组织的涡旋组成，在二维中形成六角格子 6-7。这种手性活性晶体被研究其特定属性，如边缘模式和奇异弹性。8-9

还考虑了具有粒子固有极性的强排列的活性晶体。模型证明了在相对较小的系统中传播晶体 10-11 的可能性。Dauchot 等人的实验研究了极性排列粒子的小晶体 12,13。在理论方面，Maitra 和 Ramaswamy 最近的工作 14 对定向活性固体进行了一般预测，但迄今为止尚未进行测试。

在这里，我们研究了由自推进粒子制成的更简单的活性晶体的熔化，这些粒子受到成对排斥力和旋转噪声的作用，其极性没有或弱排列。先前的工作已经得出结论或假设 KTHNY 理论仍然相关。它们通常依赖于对位置序相关函数衰减的数值估计，这是一项困难的任务，因为通常只考虑了相对较小的系统尺寸，并且通常考虑的函数通常是固有振荡 (intrinsically oscillatory) 的。下面我们绕过这个问题，这样可以获得更清晰的结果。

我们首先考虑在具有周期性边界条件的二维域中， 个 过阻尼 的 活性布朗粒子 (active Brownian particles, ABPs) 通过软谐波成对力和局部铁磁排列相互作用。它们的位置 和极性 服从：

其中 表示 方向的单位矢量，， 表示由 指向 的单位矢量，求和的范围是对所有与粒子 距离在 之内的粒子， 是均值为零的高斯噪声 (Gaussian white noise)

数值细节由 15 给出。

我们在不失一般性的情况下固定 ，主要在旋转噪声强度 (rotational noise strength) 条件下，改变剩余的两个参数，自推进速度 (self-propulsion speed) 和对齐强度 (alignment strength) 。为了设置阶段，我们首先在这个 平面上显示我们的模型的相图。足够强的排列 产生全局极性序 ()。在这个过程中，我们专注于排斥相互作用主导排列的互补区域，其中 ，极性相关性以指数形式迅速衰减。这里主要感兴趣的晶相是在较小的 或 值下发现的。在较大的 或 值下发现了液相。根据 KTHNY 的两步方案，这两个相被一个六角序区域分开。

我们的主要兴趣是晶体排列的稳定性，我们考虑最初的完美六方构型，晶格步长 ，产生平均数密度 。下面， 是近似正方形域在 中的线性大小。

晶体相中的位置序通过两点相关函数 16 的大尺度行为进行估计

其中， 是粒子 在初始完美晶格上的位置， 是其位移矢量，而 是完美晶格的倒格矢中的一个。

取向键序 (对于三角晶格的晶体来说是六角的) 的关联函数 可以类似地定义，用 代替 中的 e 指数部分，其中局部序 (local order) ， 表示粒子 的 Voronoi 邻居， 是 的方向。

关联函数 适用于最多有界位错对的构型。它们具有无固有振荡的巨大优势，因此它们的渐近行为可以很容易地测量。

正如预期的那样，在晶体区域中观察到准长程位置序，，如图 1b 所示，绘制了在增加系统尺寸时获得的结果 17。值得注意的是，在这种情况下测得的 值远大于 KTHNY 理论设定的上界 ，但存在真正的长程键序 (如图 1c 所示)，并且没有发现未结合的位错 (几乎没有发现任何结合的位错对)。在不同的对准强度 和自推进速度 值下重复这些测量，我们发现 随着 和 的增加而增加 (图 1d 和 图 1e)。有趣的是，我们观察到 (在常数 下)。

我们现在提出了一般性的论点，主张在活性晶体中观察到 的可能性。在热平衡中，线性弹性理论 18,19 给出

其中 和 都是拉梅弹性系数。熔化大致发生在熵和形成位错所需的弹性能量平衡时。这一由 KTHNY 20 精确描述的声明，给出了 熔化温度

通过消除 将这两个方程结合起来，又因为对于三角形晶格有 ，可得熔化时的

这个表达式有一个上界 。

当然，在不平衡的情况下，情况可能会大不相同。我们在这里主张并稍后展示，可以在某种意义上仍然调用两个平衡关系 (4) 和 (5)，但在两个不同的噪声水平或有效温度下：参与 (4) 的 表征晶体结构的大尺度弹性变形，参与 (5) 的 推动位错运动。假设持久极性的弱极性排列平滑了局部不规则性，可以预期 ，当 (4) 和 (5) 结合时，在熔化时产生大于 的 值。

这些想法可以通过考虑极性场 来变得更加具体。我们观察到的 表明， 可以被视为有效的空间和时间相关噪声。我们已经计算了 的空间功率谱 (图 2a)。小波数下的平台定义了对应于 的有效温度，而较低的噪声水平 可归因于 (近晶格间距) 行为 21。这些光谱可以通过将其除以 (未显示) 来叠加，而对应于 的低 平台值随着 的增加而增加 (图 2a 上的插图)。请注意，这与我们的观察结果一致，即 随 (线性) 和 增加 (图 1d 和图 1e)。

我们的大多数发现都可以在线性弹性理论中解释。控制位移场 的方程可以表示为 22,23

在平衡时， 是定义热温度 的白噪声，但在我们的活性晶体中，可以预期

其中白噪声 可以从粗粒化方程 (1) 和 (2) 中获得 。系数 和 是正的： 是一个阻尼项，用于解释极性的短程相关性，拉普拉斯算子是由弱对齐相互作用产生的。

的空间光谱可以根据线性方程 (6) (7) 计算 24。它在短波极限给出

与平衡情况相比，它给出有效温度 25。对于 ，线性理论与数值估计有很好的一致性：固定 ，改变 ，通过 和系数 重新缩放波数产生了极好的坍塌，揭示了在低 时 的行为，在高 时 的衰减 (图 2b，计算过程在 26)。乘积 提供了方程 (8) 中 因子的估计。因此，使用方程 (4)，衰变指数 可以重新表示为 。正如预期的那样，这些估计与通过直接计算 的衰变而获得的估计非常匹配 (图 2c)。

值得注意的是，基于光谱的 值甚至比直接的 值延伸得更远，图 2c 所示的情况下达到 。当 时，直接评估 会失效，因为此时具有幂律发散的布拉格峰的概念消失了 27。然而，即使晶体的畸变超过晶格间距，变形场 的行为仍然平滑，局部保持了晶体序，而六角键序明显是长程的 (图 3a 和图 3b)。正如预期的那样，结构因子在傅里叶空间呈现出六重对称性 (证明了晶体序)，但它由光滑的、非发散的、有噪音的延伸峰组成 (见图 3c)。中心区域当 时变平 ，且高 区域中的其他峰值呈现出 的缩放。所有这些都与较小 值下的行为有质的不同，其中 ，并且观察到具有幂律散度 的峰值 (图 3d)。因此，通过 获得的 的估计仍然是良好定义的。与直接测量相反，它们几乎可以延伸到接近晶体-六角晶相转变。

我们现在考虑我们的活性晶体的键序性质。在固定的参数值下，可以用来自诸如图 1c 中的数据估计渐近的非零值，即 。在平衡状态下，线性理论预测， 随温度呈指数下降 28。在这里，我们确实观察到，对于各种 值，确实有关系 (图 2d)。这表明存在“键序温度” 与重新标度的温度 与图 2d 中的线的 的预因子成比例。然而，为了成为有意义的温度，必须调整 的绝对值，使其在系统受到白噪声影响时，在极限 下，与平衡时的 一致。由于 ，我们定义了重新缩放的温度 ，并全局调整 ，使其与极限 下的 一致。然后，我们的数据显示，即使在没有对齐强度 的情况下，随着 的减少， 是如何偏离 的。类似地，，并且当 增加时， 越来越大 (图第 2e)。 可以容易地比 小得多，这表明键序的波动可以保持相当平缓，而位置序的波动很强。由于键无序有助于触发位错对的解离，因此可以认为 与 密切相关，这证实了之前介绍的情况。

我们最终证明，我们的结果并不是软相互作用势所特有的，并证实了——这已经在图 2e 中显示了——对于观察到极易变形的活性晶体，对齐 (alignment) 不是必须的。我们重新审视了 Paliwal 和 Dijkstra 最近研究的 ABPs 晶体 29，这些粒子被宣布表现出“无缺陷六角相”，这是从比 幂律衰减更快的位置序相关性得出的结论。现在粒子仅通过成对的硬核斥力相互作用。角度遵循等式 (2)，其中 ，位置遵循

其中 是在截断距离 处截断的排斥型 Weeks-Chandler-Andersen 势。使用类似于报道的“无缺陷六角序”的参数值，但考虑到更大的系统，我们发现了一个真正的晶相，具有明显的位置序代数衰变—— 值很容易大于 ——而键取向序是真正的长程的 (图 4) 30。改变 和 ，我们还是发现 仅取决于 ，而 不然 (图 4 的插图)，这表明在这种情况下，。

总之，我们已经证明，二维活性晶体中的波动可以在不熔化的情况下引起极大的变形：键序保持在长范围内，位置序代数衰减，但指数 不受 KTHNY 理论给出的 上界的限制。对于所研究的简单情况，线性弹性理论使我们能够根据定义明确的有效温度 和 来合理化我们的发现，分别量化大规模变形和键序波动。我们特别发现，预期控制熔化的 可以比 小得多，从而在保持长程键序的同时允许极端的自发变形。我们的结果的起源于这些极性的唯一时间持久性，这些结果是在活性粒子的固有轴弱排列和不弱排列的情况下获得的 31。因此我们预计，类似的现象也会出现在其他类型的活性晶体中，例如由手性或自旋单元制成的活性晶体。进一步推广，可以看出关键因素是一些有效的时间相关局部噪声的存在。根据这一想法，初步结果表明，即使是被动晶体 (passive crystals)，在受到活性浴 (active bath) 的影响时也可能表现出大幅的自发变形。

参考文献

  1. J. M. Kosterlitz and D. J. Thouless, J. Phys. C 5, L124 (1972).↩︎

  2. A. P. Young, Phys. Rev. B 19, 1855 (1979).↩︎

  3. K. J. Strandburg, Rev. Mod. Phys. 60, 161 (1988).↩︎

  4. J. Bialk ́ e, T. Speck, and H. Löwen, Phys. Rev. Lett. 108, 168301 (2012).↩︎

  5. P. Digregorio, D. Levis, L. F. Cugliandolo, G. Gonnella, and I. Pagonabarraga, Soft Matter 18, 566 (2022).↩︎

  6. I. H. Riedel, K. Kruse, and J. Howard, Science 309, 300 (2005).↩︎

  7. N. Oppenheimer, D. B. Stein, M. Y. B. Zion, and M. J. Shelley, Nat. Commun. 13, 804 (2022).↩︎

  8. B. C. van Zuiden, J. Paulose, W. T. M. Irvine, D. Bartolo, and V. Vitelli, Proc. Natl. Acad. Sci. U.S.A. 113, 12919 (2016).↩︎

  9. L. Braverman, C. Scheibner, B. VanSaders, and V. Vitelli, Phys. Rev. Lett. 127, 268001 (2021).↩︎

  10. G. Gr ́ egoire, H. Chat ́ e, and Y. Tu, Physica (Amsterdam) 181D, 157 (2003).↩︎

  11. C. Huang, L. Chen, and X. Xing, Phys. Rev. E 104, 064605 (2021).↩︎

  12. G. Briand, M. Schindler, and O. Dauchot, Phys. Rev. Lett. 120, 208001 (2018).↩︎

  13. M. N. van der Linden, L. C. Alexander, D. G. A. L. Aarts, and O. Dauchot, Phys. Rev. Lett. 123, 098001 (2019).↩︎

  14. A. Maitra and S. Ramaswamy, Phys. Rev. Lett. 123, 238001 (2019).↩︎

  15. See Supplemental Material at http://link.aps.org/ supplemental/10.1103/PhysRevLett.131.108301 for numerical details of particle simulations, phase diagram, comparison with other correlations functions, and some results of linearised theories.↩︎

  16. M. Engel, J. A. Anderson, S. C. Glotzer, M. Isobe, E. P. Bernard, and W. Krauth, Phys. Rev. E 87, 042134 (2013).↩︎

  17. Measurements performed using the more traditional, but oscillatory correlation function g(r), shown in [^48], give essentially the same results, albeit with not as clean an algebraic decay as observed with ˆ gðrÞ.↩︎

  18. D. R. Nelson and B. I. Halperin, Phys. Rev. B 19, 2457 (1979).↩︎

  19. M. Kardar, Statistical Physics of Fields (Cambridge University Press, Cambridge, England, 2007).↩︎

  20. K. J. Strandburg, Rev. Mod. Phys. 60, 161 (1988).↩︎

  21. For κ ¼ 0 (no alignment), unsurprisingly, the spatial spectrum is flat. This does not necessarily mean that TS ¼ TD, since the s field is still time correlated.↩︎

  22. C. Huang, L. Chen, and X. Xing, Phys. Rev. E 104, 064605 (2021).↩︎

  23. P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter Physics (Cambridge University Press, Cambridge, England, 1995).↩︎

  24. See Supplemental Material at http://link.aps.org/ supplemental/10.1103/PhysRevLett.131.108301 for numerical details of particle simulations, phase diagram, comparison with other correlations functions, and some results of linearised theories.↩︎

  25. In the κ ¼ 0 nonaligning limit, where the s spectrum is flat, one has a ¼ Dr, and thus TS ¼ 1 2 s20=Dr, as expected for active Brownian particles; see, e.g., [^57].↩︎

  26. See Supplemental Material at http://link.aps.org/ supplemental/10.1103/PhysRevLett.131.108301 for numerical details of particle simulations, phase diagram, comparison with other correlations functions, and some results of linearised theories.↩︎

  27. See Supplemental Material at http://link.aps.org/ supplemental/10.1103/PhysRevLett.131.108301 for numerical details of particle simulations, phase diagram, comparison with other correlations functions, and some results of linearised theories.↩︎

  28. M. Kardar, Statistical Physics of Fields (Cambridge University Press, Cambridge, England, 2007).↩︎

  29. S. Paliwal and M. Dijkstra, Phys. Rev. Res. 2, 012013(R) (2020).↩︎

  30. Note that this is only observed beyond a crossover scale of the order of 100 lattice steps, a sizable numerical difficulty, which comes in addition to having to use small timesteps.↩︎

  31. The reader might wonder why most of our results are for the rather “exotic” case of weak alignment. The main reason is practical: one can use much larger time steps in this case than with the hardcore WCA potential, saving about 1 order of magnitude of computation time.↩︎