[论文解读] Directed Percolation and the Golden Ratio

Stephan M. Dammer, Sílvio R. Dahmen+1 more|arXiv (Cornell University)|Jun 20, 2001

Theoretical and Computational Physics参考文献 4被引用 4

一句话总结

本文将杨-李零点理论应用于非平衡定向渗透模型，证明在复数渗透参数 $ p = (1 \pm \sqrt{5})/2 $（黄金分割率）及 $ p = 2 $ 处，存活概率可被精确计算。在这些取值下，存活概率在 1 与 $ (\pm\sqrt{5} - 1)/2 $ 之间呈现精确周期性振荡，揭示了非可积模型中隐藏的代数结构。

ABSTRACT

Applying the theory of Yang-Lee zeros to nonequilibrium critical phenomena, we investigate the properties of a directed bond percolation process for a complex percolation parameter p. It is shown that for the Golden Ratio p=(1(+/-)sqrt(5))/2 and for p=2 the survival probability of a cluster can be computed exactly.

研究动机与目标

• 研究非平衡系统中杨-李零点的分布，特别是如定向渗透等非可积模型中的情况。
• 探讨渗透参数 $ p $ 的特殊代数取值是否会导致尽管模型整体不可积，仍能实现精确可解。
• 将原本仅适用于平衡系统的杨-李框架，通过存活概率作为生成函数，推广至非平衡相变问题。
• 通过在复平面上特定 $ p $-值处的精确解，识别可积性或不可积性的特征信号。

提出的方法

• 在 1+1D 倾斜正方形晶格上建立定向键渗透模型，渗透概率为 $ p $，并将 $ p $ 视为复变量。
• 将存活概率 $ P(t) $ 定义为所有延伸至时间 $ t $ 的簇的加权和，权重为 $ p^n(1-p)^m $，其中 $ n $ 和 $ m $ 分别表示导电键与边界键的数量。
• 引入补全概率 $ Q(t) = 1 - P(t) $，表示未存活的簇，并根据晶格边界处的边界条件将其分解为若干子集。
• 利用簇的分解与抵消：证明当 $ p = (1 \pm \sqrt{5})/2 $ 时，两簇子集的贡献因 $ \hat{Q}_1^{(a)} + \hat{Q}_1^{(b)} = 0 $ 而代数抵消。
• 通过归纳法证明 $ Q(t) = Q(t-2) $，从而推出 $ P(t) = P(t-2) $，导致 $ P(t) $ 呈周期性振荡。
• 验证在 $ p = 2 $ 处的精确解，此时对所有 $ t > 0 $ 有 $ P(t) = 0 $，且 $ P(0) = 1 $，表明立即灭绝；在黄金分割率处，$ P(t) $ 在偶数 $ t $ 时为 1，在奇数 $ t $ 时为 $ (\pm\sqrt{5} - 1)/2 $。

实验结果

研究问题

• RQ1 尽管模型整体不可积，是否可在 $ p $ 的特殊复数值处找到定向渗透中存活概率的精确解？
• RQ2 杨-李零点方法是否能揭示代数结构（如黄金分割率）——即在非平衡系统中精确可解出现的标志？
• RQ3 簇的边界条件与抵消机制在特定 $ p $-值处实现精确解的过程中起何种作用？
• RQ4 在 $ p = (1 \pm \sqrt{5})/2 $ 处，存活概率的行为如何？这对杨-李零点的分布有何启示？
• RQ5 在定向渗透或相关模型中，是否存在其他特殊 $ p $-值，使得出现类似的精确解？

主要发现

• 在 $ p = (1 + \sqrt{5})/2 \approx 1.618 $ 处，存活概率 $ P(t) $ 在偶数 $ t $ 时振荡于 1，在奇数 $ t $ 时振荡于 $ (\sqrt{5} - 1)/2 \approx 0.618 $。
• 在 $ p = (1 - \sqrt{5})/2 \approx -0.618 $ 处，存活概率 $ P(t) $ 在偶数 $ t $ 时振荡于 1，在奇数 $ t $ 时振荡于 $ (-\sqrt{5} - 1)/2 \approx -1.618 $，但该值位于物理允许范围之外。
• 当 $ p = 2 $ 时，存活概率满足对所有 $ t > 0 $ 有 $ P(t) = 0 $，且 $ P(0) = 1 $，表明系统立即灭绝。
• 精确解源于代数抵消：当 $ p $ 取黄金分割率时，两簇子集的贡献恰好代数抵消，仅剩在 $ t-2 $ 处具有非活跃位点的构型。
• 在黄金分割率处，对所有 $ t \geq 2 $，递推关系 $ P(t) = P(t-2) $ 成立，其证明基于簇分解与边界条件约束的归纳法。
• 这些精确结果表明，$ p $ 的特殊代数取值可能作为理解非平衡临界现象中杨-李零点分布的复平面上的标志性点。