Functional limit theorems for elephant random walks on general periodic structures

基本信息

• 论文ID: 2511.10347
• 标题: Functional limit theorems for elephant random walks on general periodic structures
• 作者: Shuhei Shibata (Kyushu University)
• 分类: math.PR (Probability Theory)
• 发表时间: 2025年11月13日 (arXiv预印本)
• 论文链接: https://arxiv.org/abs/2511.10347

摘要

本文研究了一般周期结构上大象随机游走(Elephant Random Walk, ERW)的泛函极限定理，扩展了Bertenghi在标准整数格$\mathbb{Z}^d$上的结果。研究揭示了在经典$\mathbb{Z}^d$设定中不出现的新的结构依赖量，突出了底层结构如何影响游走的渐近行为。

研究背景与动机

问题背景

1. 核心问题：研究具有长程记忆的随机游走的渐近行为。大象随机游走(ERW)由Schütz和Trimper于2004年引入，用于研究一维离散时间随机游走中的记忆效应，该模型展现了从扩散到超扩散行为的相变现象。
2. 问题重要性：
   • ERW是研究记忆效应的重要模型，具有完整的历史记忆
   • 该模型在不同参数区域表现出本质不同的行为（扩散、临界、超扩散）
   • 理解结构对随机游走行为的影响具有理论和应用价值
3. 现有方法局限性：
   • 大多数现有工作集中在标准格$\mathbb{Z}^d$上
   • Bercu和Laulin、Bertenghi等人建立了$\mathbb{Z}^d$上的渐近理论
   • 缺乏对更一般周期结构（如三角格、六边形格、砖墙格等）的系统研究
4. 研究动机：
   • 将ERW理论推广到更一般的周期结构
   • 发现结构依赖的新量，揭示底层几何结构对渐近行为的影响
   • 通过Pólya型瓮模型方法建立统一的分析框架

核心贡献

1. 理论扩展：将Bertenghi在$\mathbb{Z}^d$上的泛函极限定理推广到一般周期结构，包括三角格、六边形格和砖墙格等。
2. 新结构量的发现：识别并分析了在经典$\mathbb{Z}^d$设定中不出现的结构依赖量，如协方差矩阵$\Sigma(U)$和$\Sigma(W)$，它们编码了底层几何结构的信息。
3. 统一分析框架：通过Pólya型瓮模型建立了统一的分析框架，适用于：
   • Type-I ERW（单色顶点集，$U=W$）
   • Type-II ERW（双色顶点集，$U \neq W$）
4. 完整的相图分析：在所有参数区域（扩散、临界、超扩散）建立了强大数定律和泛函极限定理，给出了精确的临界值$p_c^m = \frac{m+1}{2m}$。
5. 具体实例的计算：为多个典型结构（标准格、三角格、六边形格、砖墙格等）提供了关键量的显式计算。

方法详解

任务定义

研究对象是定义在一般周期结构$\Gamma$上的大象随机游走$\{S_n\}_{n=0}^{\infty}$：

输入：

• 步向量集合$U = \{u_1, \ldots, u_m\} \subset \mathbb{R}^d$和$W = \{w_1, \ldots, w_{m'}\} \subset \mathbb{R}^d$
• 记忆参数$p \in (0,1)$
• 初始步向量$u_{i_0}$和$w_{j_0}$

输出：

• 强大数定律：$\frac{S_n}{n} \to \frac{1}{2}(\bar{u} + \bar{w})$ a.s.
• 泛函极限定理：在不同参数区域的适当归一化下的极限分布

约束条件：

• $\Gamma$必须是$\mathbb{R}^d$中的格
• 对于Type-II ERW，需满足交替规则（bipartite结构）

模型架构

1. 状态空间构造

Case U = W（Type-I ERW）： $\Gamma = \left\{\sum_{i=1}^m k_i u_i : k_i \in \mathbb{N} \cup \{0\}\right\}$ 单色顶点集，所有顶点结构等价。

Case U ≠ W（Type-II ERW）： 首先定义： \Gamma_0 = \left\{\sum_{i=1}^m k_i u_i + \sum_{j=1}^{m'} l_j w_j : \sum_{i=1}^m k_i = \sum_{j=1}^{m'} l_j, k_i, l_j \in \mathbb{N} \cup \{0}\right\}

然后： $\Gamma = \Gamma_0 \sqcup (\Gamma_0 + U)$

定义顶点类：$Z_U := \Gamma_0$和$Z_W := \Gamma_0 + U$，形成二部图结构。

2. ERW的定义

Type-II ERW（$U \neq W$）：

• 位置定义： $S_{2n} = \sum_{i=1}^n (\sigma_i + \tau_i), \quad S_{2n-1} = S_{2(n-1)} + \sigma_n$
• 步选择机制：在时刻$n+1$，给定历史$\sigma_1, \ldots, \sigma_n$和$\tau_1, \ldots, \tau_n$： $P(\sigma_{n+1} = \sigma_{U_n}) = p, \quad P(\sigma_{n+1} = \sigma) = \frac{1-p}{m-1} \text{ for } \sigma \in U \setminus \{\sigma_{U_n}\}$ 其中$U_n$在$\{1,2,\ldots,n\}$上均匀分布。

Type-I ERW（$U = W$）： $S_n = \sum_{i=1}^n \sigma_i$ 步选择机制类似，但从整个历史中采样。

技术创新点

1. Pólya瓮模型的连接

建立ERW与瓮过程的分布等价：

Type-II ERW： $\{S_{2n}\}_{n=0}^{\infty} \stackrel{d}{=} \left\{\sum_{i=1}^m X_n^i u_i + \sum_{j=1}^{m'} Y_n^j w_j\right\}_{n=0}^{\infty}$

Type-I ERW： $\{S_{2n}\}_{n=0}^{\infty} \stackrel{d}{=} \left\{\sum_{i=1}^m X_{2n}^i u_i\right\}_{n=0}^{\infty}$

其中$X_n$和$Y_n$是独立的Pólya瓮过程。

2. 替换矩阵分析

关键的$m \times m$替换矩阵： $A = \frac{1-p}{m-1}J_m + \frac{mp-1}{m-1}I_m$

特征值：

• $\lambda_1 = 1$（重数1）
• $\lambda_2 = \frac{mp-1}{m-1}$（重数$m-1$）

临界值由$\lambda_2/\lambda_1 = 1/2$确定： $p_c^m = \frac{m+1}{2m}$

3. 结构依赖协方差矩阵

定义关键的$d \times d$协方差矩阵： $\Sigma(U) = \frac{1}{m}\sum_{i=1}^m (u_i - \bar{u})(u_i - \bar{u})^{\top}$$\Sigma(W) = \frac{1}{m'}\sum_{j=1}^{m'} (w_j - \bar{w})(w_j - \bar{w})^{\top}$

这些矩阵编码了步向量的几何配置，在$\mathbb{Z}^d$情况下简化为$I_d/d$，但在一般结构中可能是非对角的，反映了不同坐标轴上扩散的相关性。

4. 与Janson一般理论的联系

本文巧妙地应用了Janson 2004关于多型分支过程和广义Pólya瓮的泛函极限定理（特别是Theorem 3.21, 3.24, 3.31），并通过连续映射定理将瓮过程的渐近行为转化为ERW的渐近行为。

实验设置

理论验证

本文是纯理论研究，不涉及数值实验或数据集。验证通过严格的数学证明完成。

典型实例计算

论文在第6节提供了6个具体实例的关键量计算：

1. 标准格$\mathbb{Z}^d$：$m=2d$, $p_c^m = \frac{2d+1}{4d}$, $\Sigma(U) = I_d/d$
2. 三角格：$m=6$, $p_c^m = 7/12$, $\Sigma(U) = I_2/2$
3. 六边形格：$m=m'=3$, $p_c^m = 2/3$, $\Sigma(U) = I_2/2$
4. 砖墙格：$m=m'=3$, $p_c^m = 2/3$, $\Sigma(U) = \frac{2}{9}\begin{pmatrix}3 & 0\\0 & 1\end{pmatrix}$
5. 混合结构1：$U = \{\pm u_1, \pm u_2\}$, $W = \{\pm e_1, \pm e_2\}$，$m=m'=4$
6. 混合结构2：$U = \{\pm e_1, \pm e_2, (1,2)^{\top}\}$, $W = \{\pm e_1, \pm e_2\}$，$m=5, m'=4$

主要结果

1. 强大数定律（Theorem 5.1）

对所有$p \in (0,1)$： $\frac{S_n}{n} \to \frac{1}{2}(\bar{u} + \bar{w}) \quad \text{a.s. as } n \to \infty$

关键特点：

• 极限不依赖于边数$m, m'$
• 仅依赖于步向量的平均值
• 对Type-I和Type-II ERW都成立

2. 扩散区域（Theorem 5.3）

条件：$0 < p < p_c^m \leq p_c^{m'}$

结果： $\left\{\frac{S_{\lfloor 2nt \rfloor} - nt(\bar{u} + \bar{w})}{\sqrt{n}}\right\}_{t \geq 0} \Rightarrow \{W_t\}_{t \geq 0}$

其中$\{W_t\}_{t \geq 0}$是中心化的$\mathbb{R}^d$值连续高斯过程，协方差结构为： $\mathbb{E}[W_s W_t^{\top}] = C_a s\left(\frac{t}{s}\right)^a \Sigma(U) + C_{a'} s\left(\frac{t}{s}\right)^{a'} \Sigma(W)$

其中$a = \frac{mp-1}{m-1}$, $C_a = \frac{1}{1-2a}$。

创新点：

• 协方差矩阵$\Sigma(U), \Sigma(W)$体现结构依赖性
• 在$\mathbb{Z}^d$情况下退化为Bertenghi的结果
• 当$\Sigma(U), \Sigma(W)$非对角时，不同坐标轴的扩散相关

3. 临界区域（Theorem 5.5）

条件：$0 < p_c^m = p = p_c^{m'}$

结果： $\left\{\frac{S_{\lfloor 2nt \rfloor} - nt(\bar{u} + \bar{w})}{n^{t/2}\sqrt{\log n}}\right\}_{t \geq 0} \Rightarrow \{W_t\}_{t \geq 0}$

协方差结构： $\mathbb{E}[W_s W_t^{\top}] = s\Sigma(U, W)$

其中$\Sigma(U, W) = \Sigma(U) + \Sigma(W)$。

特点：

• 需要额外的$\sqrt{\log n}$归一化因子
• 极限过程可表示为$W_t = \Sigma(U,W)^{1/2}B_t$，其中$B_t$是标准布朗运动
• 在$\mathbb{Z}^d$情况下，$\{W_t/\sqrt{d}\}$是标准布朗运动

4. 超扩散区域（Theorem 5.8）

条件：$0 < p_c^m = p_c^{m'} < p$

Type-II ERW结果： $\left\{\frac{S_{\lfloor 2nt \rfloor} - nt(\bar{u} + \bar{w})}{n^a}\right\}_{t \geq 0} \Rightarrow \{t^a L\}_{t \geq 0}$

Type-I ERW结果： $\left\{\frac{S_{\lfloor nt \rfloor} - nt\bar{u}}{n^a}\right\}_{t \geq 0} \Rightarrow \{t^a \tilde{L}\}_{t \geq 0}$

其中$L, \tilde{L}$是非零的$\mathbb{R}^d$值随机向量。

重要观察：

• Type-I和Type-II ERW在超扩散区域有不同的极限分布
• 极限分布依赖于初始步的选择
• 在假设$\bar{u} = \bar{w} = 0$和初始步均匀选择下： $\mathbb{E}[LL^{\top}] = \frac{1}{(2a-1)\Gamma(2a)}\Sigma(U,W)$

5. 混合参数区域（Remarks 5.7, 5.9）

论文还讨论了$p_c^m \neq p_c^{m'}$的情况：

• 当$0 < p_c^m = p < p_c^{m'}$时，$X$过程主导，极限仅依赖于$\Sigma(U)$
• 当$0 < p_c^m < p_c^{m'} < p$时，需根据$a$和$a'$的大小关系确定主导项

相关工作

一维ERW的研究

• Schütz和Trimper 2004：引入ERW模型，发现$p=3/4$的相变
• Baur和Bertoin 2016：建立ERW与Pólya瓮的联系
• Bercu 2017：鞅方法分析ERW
• Coletti等 2017：中心极限定理
• Kubota和Takei 2019：超扩散区域的高斯波动

多维ERW的研究

• Bercu和Laulin 2019：通过鞅方法研究MERW的渐近行为
• Bertenghi 2022：在$\mathbb{Z}^d$上建立泛函极限定理（本文直接推广的工作）
• González-Navarrete 2020：具有随机倾向的多维游走
• Chen和Laulin 2023：平滑遗忘增强的多维ERW
• Curien和Laulin 2024：平面ERW的常返性
• Qin 2025：多维ERW的常返性和暂态性

碰撞问题

• Roy, Takei和Tanemura 2024：$\mathbb{Z}$上两个ERW的碰撞问题
• Shibata和Shirai 2025：不同记忆参数情况下的碰撞问题和距离渐近

Pólya瓮理论

• Janson 2004：多型分支过程和广义Pólya瓮的泛函极限定理（本文的核心理论工具）
• Athreya和Karlin 1968：瓮模型嵌入连续时间分支过程
• Chauvin等 2011：大型Pólya瓮的极限分布

本文的优势

1. 一般性：不限于$\mathbb{Z}^d$，覆盖广泛的周期结构
2. 结构依赖性：揭示了新的结构依赖量
3. 统一框架：Type-I和Type-II ERW的统一处理
4. 完整性：覆盖所有参数区域的完整理论

结论与讨论

主要结论

1. 理论推广成功：成功将Bertenghi在$\mathbb{Z}^d$上的结果推广到一般周期结构，证明了Pólya瓮方法的普适性。
2. 结构影响的量化：通过协方差矩阵$\Sigma(U)$和$\Sigma(W)$精确量化了底层几何结构对ERW渐近行为的影响。
3. 相变现象的保持：临界值$p_c^m = \frac{m+1}{2m}$由替换矩阵的特征值决定，相变现象在一般结构中依然存在。
4. Type-I与Type-II的差异：除超扩散区域外，两类ERW有相同的极限分布；在超扩散区域，初始条件的影响导致不同的极限。

局限性

1. 结构限制：
   • 仅考虑可表示为格的周期结构
   • 对于多部图（$l \geq 3$），如kagome格，需要不同的分析技术
   • 排除了$m=1$的平凡情况
2. 初始条件：
   • 假设从原点出发
   • 前两步是确定性的
   • 超扩散区域的结果依赖于初始步的分布
3. 参数范围：
   • 仅考虑$p \in (0,1)$
   • $p=1$对应平凡情况（确定性游走）
   • $p=0$未讨论
4. 分布细节：
   • 超扩散区域的极限随机变量$L$的分布未完全刻画
   • 仅给出了前两阶矩的表达式

未来方向

1. 扩展到更一般的结构：
   • 非周期结构
   • 多部图（$l \geq 3$）如kagome格
   • 随机环境中的ERW
2. 极限分布的精细分析：
   • 超扩散区域$L$的完整分布
   • 应用Guérin等2023, 2025的不动点方程方法
3. 碰撞问题：
   • 一般结构上两个ERW的碰撞问题
   • 不同记忆参数的情况
4. 其他性质：
   • 常返性和暂态性的完整刻画
   • 大偏差原理
   • 局部时的研究

深度评价

优点

1. 数学严谨性：
   • 证明完整且严格，充分利用了Janson的一般理论
   • 对不同参数区域的分类清晰完整
   • 技术细节处理得当（如协方差计算）
2. 理论贡献：
   • 首次系统研究一般周期结构上的ERW
   • 发现了结构依赖的新量$\Sigma(U), \Sigma(W)$
   • 建立了Type-I和Type-II ERW的统一框架
3. 清晰度：
   • 论文组织合理，从简单到复杂
   • 图示（三角格、六边形格、砖墙格）直观
   • 第6节的具体实例计算非常有用
4. 完整性：
   • 覆盖所有参数区域（扩散、临界、超扩散）
   • 讨论了多种混合参数情况
   • 与现有文献的联系清楚

不足

1. 实例的有限性：
   • 虽然理论一般，但仅给出6个具体实例
   • 缺少对某些有趣结构（如kagome格）的讨论
   • 可以增加更多非标准结构的例子
2. 物理直觉：
   • 缺少对结构依赖性的物理或几何直觉
   • $\Sigma(U), \Sigma(W)$的非对角元素的意义未充分解释
   • 不同结构导致不同行为的机制讨论不够
3. 超扩散区域：
   • 极限分布$L$的刻画不完整
   • 仅给出二阶矩，更高阶矩或完整分布未知
   • Type-I与Type-II差异的深层原因未充分阐明
4. 应用讨论：
   • 缺少对实际应用场景的讨论
   • 这些结果对理解真实系统有何启示未提及

影响力

1. 对领域的贡献：
   • 为ERW研究开辟了新方向（一般结构）
   • 提供了可用于其他记忆增强随机游走的方法论
   • 丰富了Pólya瓮理论的应用
2. 实用价值：
   • 理论性强，但可为研究材料中的扩散、网络上的随机过程等提供理论基础
   • 第6节的公式可直接用于具体计算
3. 可复现性：
   • 证明完整，可验证
   • 关键量的计算公式明确
   • 理论结果可应用于新结构

适用场景

1. 数学研究：
   • 概率论中的极限定理
   • 随机过程理论
   • 组合随机结构
2. 物理系统：
   • 晶格上的扩散过程
   • 具有记忆的输运现象
   • 相变现象的研究
3. 网络科学：
   • 图上的随机游走
   • 信息传播模型
   • 网络探索算法
4. 统计物理：
   • 非马尔可夫过程
   • 长程相关系统
   • 异常扩散

参考文献（精选）

• 2 Baur & Bertoin (2016): Elephant random walks and their connection to Pólya-type urns. Physical Review E.
• 5 Bercu & Laulin (2019): On the multi-dimensional elephant random walk. J. Stat. Phys.
• 7 Bertenghi (2022): Functional limit theorems for the multi-dimensional elephant random walk. Stoch. Models.
• 17 Janson (2004): Functional limit theorems for multitype branching processes and generalized Pólya urns. Stochastic Process. Appl. （核心理论工具）
• 23 Schütz & Trimper (2004): Elephants can always remember: Exact long-range memory effects in a non-Markovian random walk. Physical Review E. （ERW起源）

总体评价：这是一篇高质量的概率论理论论文，成功地将ERW理论推广到一般周期结构，揭示了结构对渐近行为的深刻影响。数学上严谨完整，理论贡献显著，为该领域的进一步研究奠定了基础。主要不足在于对结果的物理直觉和应用场景的讨论不够充分，以及超扩散区域极限分布的刻画有待深化。