Dynamical renormalization group analysis of $O(n)$ model in steady shear flow

通过将强各向异性纳入动力学重正化群分析，本研究确定了$O(n)$模型在稳态剪切流下的一个新稳定高斯不动点，揭示了剪切流稳定了二维中的长程有序，并改变了守恒与非守恒序参量的上临界维度，从而违背了平衡态的霍亨伯格-梅尔明-瓦格纳定理。

要点

想象一个拥挤的舞池，每个人都试图同步移动。在一个平静的房间（平衡态）里，如果音乐停止，舞者们可能会僵在原地；或者如果他们试图形成某种图案，仅凭人群相互碰撞的巨大数量，他们就可能被挤散。在物理学中，这就像一种材料试图决定自己是处于有序状态（如磁铁）还是无序状态（如气体）。

现在，想象有人开始将整个舞池向侧面推动，形成稳定的“剪切流”。舞者们不再只是随机碰撞；他们被扫向一个特定的方向。本文问道：这种持续的推动如何改变舞者们的自我组织方式？

作者 Harukuni Ikeda 和 Hiroyoshi Nakano 使用了一种名为“重整化群”的复杂数学工具（将其想象为一台显微镜，通过放大和缩小来观察不同尺度下图案的变化）对此进行了研究。他们观察了两类舞者：

1. 模型 A：可以自由移动并轻松更换位置的舞者（非守恒）。
2. 模型 B：被困在网格中且只能与邻居交换位置的舞者（守恒）。

以下是主要发现，以简单的方式解释：

1. 推动的“魔力”

在正常、平静的房间里，关于舞池必须多小才会导致舞者无法形成大规模有序图案，存在着严格的规则。

• 旧规则：在一个二维房间（如平面地板）中，如果舞者试图打破对称性（例如选择特定的朝向），霍亨伯格 - 梅尔明 - 瓦格纳定理指出这是不可能的。随机的推挤太强了，图案会被破坏。你至少需要一个三维房间才能实现这一点。
• 新发现：作者发现，当你施加那种稳定的“推动”（剪切流）时，规则完全改变了。这种推动实际上稳定了图案。即使在平坦的二维房间里，舞者们现在也能形成完美的长程有序。“推动”抑制了通常破坏派对的混乱推挤。

2. “新常态”（不动点）

在物理学中，系统通常会稳定在一个“不动点”——即无论你将视野放大还是缩小，游戏规则停止变化的状态。

• 没有推动时：系统试图稳定在一个“高斯不动点”（一种标准的、可预测的状态），但推动使这种状态变得不稳定。这就像试图在有人摇晃桌子时，将铅笔平衡在笔尖上。
• 有推动时：作者发现了一个新的、稳定的不动点。由于推动非常强烈，系统找到了一种新的平衡方式。这种新状态是“高斯”的（简单且可预测），但其行为与平静状态截然不同。

3. 维度的收缩

本文引入了两个关键数值：

• 上临界维度 ($d_{up}$)：指“简单”规则（平均场理论）开始完美起作用的房间大小。

  • 之前：你需要一个四维房间，简单规则才能生效。
  • 之后：有了推动，简单规则甚至在二维房间（对于模型 A）甚至零维房间（对于模型 B，这意味着它们在任何地方都适用）中也能生效。
  • 类比：这就像推动让舞者们如此协调，以至于即使他们身处狭小拥挤的空间，其行为也仿佛处于一个更大、更简单的世界。

• 下临界维度 ($d_{low}$)：指有序成为可能的最小房间尺寸。

  • 之前：你需要一个大于二维的房间才能拥有有序。
  • 之后：有了推动，即使房间小于二维（数学上表明 $d_{low} < 2$），有序也是可能的。
  • 类比：推动在组织人群方面如此有效，以至于即使是在狭窄到他们无法正常站立的走廊里，他们也能排成一行。

4. “拉伸”效应

最有趣的视觉变化是舞者的移动方式。

• 在平静的房间里：如果你观察舞者之间的距离，各个方向都是一样的（各向同性）。
• 在推动中：舞者们被拉伸。沿着推动的方向，他们变得非常长而细；垂直于推动的方向，他们保持短小。
• 结果：“关联”（一个舞者的动作预测另一个舞者动作的程度）发生了变化。在推动的方向上，连接变得较弱，并遵循一种奇怪的分数幂律（例如 $1/|q|^{2/3}$ 而不是通常的 $1/|q|^2$）。这就像舞者们手拉手形成了一条长长的、被拉伸的链条，而不是一个紧密的圆圈。

5. 为什么之前的实验感到困惑

作者提到，过去的计算机模拟得出了令人困惑的结果。有人说序参量（群体的组织程度）是 0.37，有人说是 0.48，而“简单”理论预测的是 0.5。

• 解释：作者认为，“拉伸”（各向异性）如此极端，以至于标准的计算机模拟规模不够大，无法看到真正的图案。
• 类比：想象试图给一条非常长而细的蛇拍照。如果你的相机画幅是正方形的，你可能会切掉尾巴或头部，使其看起来像一条短而粗的蠕虫。要看到整条蛇，你需要一个宽度是高度 100 倍的相机。作者认为，之前的模拟在“蛇状”系统上使用了“方形相机”，导致了错误的测量。

总结

本文声称，稳定的剪切流就像一个强大的组织者。它打破了旧有的物理规则，即“你无法在二维中拥有有序”。相反，这种流动创造了一种新的、稳定的状态，在其中有序更容易实现，规则变得更简单（平均场），并且系统沿流动方向发生剧烈的拉伸。作者认为，这解释了为什么一些实验观察到“平均场”行为，而另一些实验感到困惑——他们只是没有考虑到这种极端的拉伸。

技术摘要

问题陈述
$O(n)$ 模型在稳态剪切流下的临界行为仍是非平衡统计力学中的一个未决问题。虽然平衡态临界现象已通过重整化群（RG）方法得到充分理解，但剪切流的引入破坏了旋转对称性，并将系统驱动至非平衡态。先前的理论工作，如德热纳（De Gennes）的平均场分析以及大贯（Onuki）和川崎（Kawasaki）对模型 H 的 RG 研究，表明剪切流会抑制沿流动方向的临界涨落，并可能改变临界维度。然而，关于剪切伊辛（Ising）模型和$O(n)$ 模型的数值模拟，在临界指数（特别是$\beta$）以及二维（$d=2$）中长程有序的存在性方面得出了相互矛盾的结果。一个关键的理论挑战在于，RG 分析中缺乏一个能够正确纳入剪切流诱导的强各向异性的稳定不动点，从而导致理论预测与数值观测之间存在差异。

方法论
作者采用了一种专为各向异性系统设计的动力学重整化群（RG）分析。他们将此框架应用于稳态剪切流下的$O(n)$ 模型，针对两种不同情形：

1. 模型 A：非守恒序参量动力学。
2. 模型 B：守恒序参量动力学。

核心方法论的创新在于采用了各向异性标度假设。与早期假设各向同性标度或仅通过重整化输运系数处理各向异性的分析不同，本工作明确假设平行于流动方向（$x_1$）和垂直于流动方向（$x_\perp$）的坐标具有不同的标度指数。标度变换定义为：
$$x_1 = l^\zeta x'_1, \quad x_\perp = l x'_\perp, \quad t = l^z t', \quad \phi_a = l^\chi \phi'_a$$
其中$\zeta$、$z$和$\chi$是独立的临界指数。作者通过要求特定物理量（如剪切率和噪声强度）的标度不变性，推导了剪切率（$\dot{\gamma}$）、输运系数、噪声强度及相互作用项的 RG 流方程，以识别稳定的高斯不动点。

主要贡献与结果
本文的主要贡献是识别出了一个新的稳定高斯不动点，该不动点仅在稳态剪切流条件下存在。与在任何维度下都会被剪切流破坏的平衡态高斯不动点不同，该不动点对于剪切率是稳定的。

分析得出了以下具体结果：

• 临界维度：

  • 上临界维度（$d_{up}$）：本文发现模型 A 的$d_{up} = 2$，模型 B 的$d_{up} = 0$。这意味着在物理维度$d=2$和$d=3$下，两个模型均表现出平均场临界指数。
  • 下临界维度（$d_{low}$）：分析显示模型 A 的$d_{low} = 0$，模型 B 的$d_{low} = -2$。关键在于，这表明即使在$d=2$下也能发生连续对称性破缺（从而产生长程有序），而在平衡态系统中，霍亨伯格 - 梅尔明 - 瓦格纳（Hohenberg-Mermin-Wagner）定理禁止此类有序的存在。

• 临界指数：

  • 对于模型 A，指数为$\zeta=3$，$z=2$，且$\chi = -d/2$。序参量标度指数对于所有$d \geq 2$均为负值，证实了长程有序的稳定性。
  • 对于模型 B，指数为$\zeta=5$，$z=4$，且$\chi = -(d+2)/2$。
  • 在这两种情况下，序参量的临界指数均为$\beta = 1/2$，与平均场理论一致。

• 关联函数：

  • 关联函数表现出独特的各向异性标度。在傅里叶空间中，沿流动方向（$q_1$）的关联函数对于模型 A 标度为$C(q_1) \sim |q_1|^{-2/3}$，对于模型 B 标度为$C(q_1) \sim |q_1|^{-2/5}$。
  • 这些分数幂（小于平衡态值 2）表明，与平衡态相比，沿流动方向的临界涨落受到了更强的抑制。这种抑制在模型 B 中比在模型 A 中更为显著。

意义与主张
本文主张，围绕剪切模型中平均场特征的“争议”之所以得以解决，源于对标度假设中各向异性的正确处理。通过识别具有这些特定各向异性指数的稳定高斯不动点，作者为以下现象提供了理论基础：

1. 在$d=2$和$d=3$中观察到平均场指数（$\beta=1/2$）。
2. 在$d=2$下连续对称性破缺中长程有序的存在，这在平衡态中是不可能的。
3. 解释了为何先前的数值模拟可能难以收敛到这些值：巨大的各向异性指数（特别是模型 B 中的$\zeta=5$）意味着模拟需要极端的纵横比（$L_\parallel \sim L_\perp^5$）才能捕捉到正确的标度，而这一条件在标准的有限尺寸标度分析中往往无法满足。

作者将这项工作定位为协调理论 RG 预测与数值结果的必要步骤，并指出早期分析中“缺失”的稳定不动点是由于在标度假设中忽视了强各向异性所致。他们明确指出，模型 A 和 B 的结果与大贯和川崎推导的模型 H 的无限剪切率极限不同，并将此归因于不同的理论形式（显式各向异性标度与各向异性系数的各向同性标度），同时建议未来的工作应致力于弥合这些框架。