Occupation-condensation transition of a sublinearly vertex-reinforced random walk on regular tree

技术摘要：正则树上亚线性顶点增强随机游走的占据凝聚相变

问题陈述
本文研究了在根为 $b$ 进制树上的亚线性顶点增强随机游走（VRRW）的长期行为。在该过程中，步向相邻节点 $u$ 的概率与权重 $w(n_u) = 1 + \beta n_u^a$ 成正比，其中 $n_u$ 是之前的访问次数，$\beta \ge 0$ 是增强强度，$a \in (0, 1)$ 控制记忆增长率。

核心挑战在于理解两种力量之间的竞争：树的几何熵驱动（由于指数级增长将游走推向外部）与增强机制（将游走拉回先前访问过的位点）之间的竞争。虽然在 $\mathbb{Z}$ 上的线性增强（$a=1$）会导致几乎处处在有限集合上的局部化，且树上的边增强表现出复现/逃逸的相变，但树上的亚线性机制（$0 < a < 1$）在很大程度上仍是一个开放问题。作者旨在确定游走是局部化（范围有界）还是经历另一种类型的相变，即占据度（occupation measure）发生凝聚而范围继续增长。

方法论
作者结合了基于可逆马尔可夫链和电网络理论的解析论证，进行了广泛的数值模拟。

1. 模拟框架：

   • 在分支数为 $b \in \{2, 3, 4, 5\}$ 的根树上对游走进行模拟。
   • 运行时间延长至 $t = 3 \times 10^7$ 步。
   • 追踪两个主要的观测量：
     • 扩散（Spreading）： 范围 $R(t)$（访问过的顶点数）和范围指数 $\nu_R = \frac{d \log E[R(t)]}{d \log t}$。
     • 浓度（Concentration）： 凝聚分数 $m_1(t) = \max_v n_v(t)/t$，代表在访问次数最多的顶点处停留的时间比例。
   • 统计分析包括自助法（bootstrap）置信区间（高达 120 次实现）以及用于表征相变的 Binder 累积量和易受性（susceptibilities）计算。

2. 解析方法：

   • 冻结环境近似： 作者在权重固定的“冻结”环境中分析该链。他们建立该链是可逆的，其边电导度为 $c_{uv} = w_u w_v$，平稳测度为 $\mu_v \propto w_v W_v$，其中 $W_v = \sum_{u \sim v} w_u$。
   • 核心轮廓（Core Profile）： 他们基于可逆测度，推导了“凝聚核心”（具有正极限占据度的顶点）的非线性递归关系。
   • 前沿分析（Frontier Analysis）： 利用树上偏置游走的的分支数判据，通过分析访问集增长前沿处的漂移来估计相变点。

主要贡献与结果

1. 发现凝聚相变：
   本文识别出了在有限临界增强强度 $\beta_c(a, b)$ 处存在的尖锐相变。

   • 低于 $\beta_c$ 时： 占据度扩散；访问次数最多的顶点所占的时间比例趋于零，且范围呈线性增长（$\nu_R \to 1$）。
   • 高于 $\beta_c$ 时： 占据度发生凝聚。单个顶点占据了 $O(1)$ 的时间比例（在观测时间内保持稳定），而范围继续增长，但极其缓慢。其增长更符合 $\log t$ 而非任何幂律描述（$\nu_R \approx 0.1$）。
   • 关键区别： 作者明确区分了这种“凝聚”与“有限范围局部化”。他们发现没有证据表明范围变得有界；游走会继续发现新的顶点，尽管是以对数速度缓慢进行的。

2. 刻画相变线：

   • 相变通过四个独立的估计量定位：凝聚序参数 $m_1$、易受性峰值、Binder 累积量交叉以及范围指数 $\nu_R = 1/2$。
   • 相变线 $\beta_c(a, b)$ 在高达 $t = 3 \times 10^7$ 的观测时间内未表现出系统性漂移。
   • 与分支数的标度关系： 临界值满足 $\beta_c \propto b - 1$。当除以 $b-1$ 后，对于 $b=2, 3, 4, 5$ 的相变线在百分之几的误差范围内塌缩到同一条曲线上。
   • 对 $a$ 的依赖性： 该线随 $a$ 增加而陡峭下降，在 $a \approx 0.8$ 附近消失（立即凝聚）。随着 $a$ 减小，$\beta_c$ 急剧上升。数据表明在 $b=2$ 时，在 $a \approx 0.42$ 附近存在发散，尽管作者并未证实存在严格的阈值 $a_c > 0$。

3. 可逆性与核心结构：

   • 在凝聚相中，高占据度的核心遵循冻结环境的可逆平稳测度：$\mu_v \propto w_v \sum_{u \sim v} w_u$。
   • 这导致了沿树骨架的占据轮廓的非线性三项递归。
   • 指数 $a=1/2$ 是核心轮廓形状的一个临界情况（允许几何尾部），但它不是相变阈值。相变由前沿动力学决定，这使得树的阈值与 $\mathbb{Z}$ 上 VRRW 已知的 Chen-Kozma 值（$a=1/2$）相解耦。

4. 相变的性质：

   • 在 $\beta_c$ 附近，系统表现出非自平均性和 $m_1$ 分布的双峰性，并具有负的 Binder 累积量。
   • 作者将其描述为一种“共存型现象学”（暗示是一阶相变），但保持了谨慎，指出有限的观测时间也可能反映了宽泛的捕获/逃逸时间分布。他们并未声称这是一个已证实的的一阶相变。

意义与主张
本文声称解决了正则树上增强驱动凝聚的机制，提供了一个可计算的模型，其中增强强度不会像在平均场完全图中那样被标度掉。

• 机制： 相变是由熵逃逸率（由分支数 $b$ 决定）与记忆增长率（由 $a$ 和 $\beta$ 决定）之间的竞争驱动的。
• 几何 vs. 记忆： 研究表明相图在 $(a, \beta)$ 空间中是二维的，且几何结构（分支）是必不可少的；平均场近似无法捕捉到这一相变。
• 凝聚 vs. 局部化： 本文的主要理论贡献是区分了占据凝聚（支持集保持无界但质量集中）与有限范围局部化。作者认为，在树上，游走发生凝聚，但在范围有界的意义上并未实现局部化。
• 标度律： 本文确立了标度律 $\beta_c \propto b - 1$，从而组织了不同树拓扑下的相图，并表明临界阈值受分支数控制，而非与 1D 格子相关的临界指数 $a=1/2$ 相关。

作者总结道，虽然核心轮廓受可逆性支配，但相变位置由前沿动力学决定，且关于小 $a$ 边界的性质以及相变的精确热力学阶数仍是未来工作的开放问题。