[论文解读] Central limit theorem for first-passage percolation time across thin cylinders

Sourav Chatterjee, Partha S. Dey|arXiv (Cornell University)|Nov 30, 2009

Stochastic processes and statistical mechanics参考文献 16被引用 9

一句话总结

该论文为在 $d$ 维格点中穿过细长圆柱体 $[0,n] \times [-h_n,h_n]^{d-1}$ 的首达通过时间建立了中心极限定理（CLT），证明当 $h_n = o(n^{1/(d+1)})$ 时存在高斯波动。在波动指数和横向指数存在且极限形状严格凸的附加假设下，该 CLT 在二维和三维中可扩展至无限制测地线高度，二维中存在数值证据支持该猜想。

ABSTRACT

We prove that first-passage percolation times across thin cylinders of the form $[0,n] imes [-h_n,h_n]^{d-1}$ obey Gaussian central limit theorems as long as $h_n$ grows slower than $n^{1/(d+1)}$. It is an open question as to what is the fastest that $h_n$ can grow so that a Gaussian CLT still holds. Under the natural but unproven assumption about existence of fluctuation and transversal exponents, and strict convexity of the limiting shape in the direction of $(1,0,...,0)$, we prove that in dimensions 2 and 3 the CLT holds all the way up to the height of the unrestricted geodesic. We also provide some numerical evidence in support of the conjecture in dimension 2.

研究动机与目标

• 在 $\mathbb{Z}^d$ 中为穿过细长圆柱体的首达通过时间建立中心极限定理（CLT）。
• 确定 CLT 仍然成立时圆柱体高度 $h_n$ 的最大增长率。
• 在波动指数和横向指数存在且极限形状严格凸的自然假设下，研究 CLT 是否可扩展至二维和三维中无限制测地线的高度。
• 提供数值证据，支持二维中方差按 $\sim n h_n^{-\gamma}$ 规律缩放的猜想，其中 $\gamma = 1/2$。

提出的方法

• 通过将通过时间沿圆柱体长轴分解为块，应用鞅和集中不等式。
• 应用李雅普诺夫条件及对中心矩的矩界，以证明收敛到正态分布。
• 利用递归估计和等周不等式，对方差进行上下界控制，以控制波动。
• 依赖形状定理及波动指数和横向指数的假设，将 CLT 扩展至更大的 $h_n$。
• 通过伯努利边权的数值模拟，估计方差缩放并利用 QQ 图检验渐近正态性。
• 对对数方差进行线性回归，以估计猜想缩放 $\operatorname{Var}(a_n(h_n)) \approx \beta n h_n^{-\gamma}$ 中的指数 $\gamma$。

实验结果

研究问题

• RQ1 当高度 $h_n$ 的增长速度慢于 $n^{1/(d+1)}$ 时，穿过细长圆柱体的首达通过时间是否满足中心极限定理？
• RQ2 在 $d$ 维格点中，Gaussian CLT 仍然成立的 $h_n$ 最快可能增长率是多少？
• RQ3 在二维和三维中，在波动指数和横向指数存在且极限形状严格凸的假设下，CLT 是否可扩展至无限制测地线的高度？
• RQ4 在二维中，首达通过时间的方差是否渐近正比于 $n h_n^{-\gamma}$，其中 $\gamma = 1/2$？
• RQ5 模拟结果是否支持 $h_n = n^{\alpha}$ 且 $\alpha \leq 2/3$ 时通过时间的猜想渐近正态性？

主要发现

• 当 $h_n = o(n^{1/(d+1)})$ 时，首达通过时间穿过细长圆柱体 $[0,n] \times [-h_n,h_n]^{d-1}$ 的中心极限定理成立。
• 在波动指数和横向指数存在且极限形状严格凸的假设下，二维和三维中 CLT 可扩展至无限制测地线的高度。
• 二维中的数值模拟支持猜想：当 $h_n = n^{\alpha}$ 且 $\alpha \leq 2/3$ 时，方差按 $\sim n h_n^{-1/2}$ 缩放，不同 $p$ 值下估计的 $\gamma \approx 1/2$。
• 在 $n = h_n^2 = 3025$ 且 $p \in \{0.6, 0.7, 0.8, 0.9\}$ 的模拟中，QQ 图显示与正态分布高度一致，支持渐近正态性。
• 该论文证明：当波动指数 $\chi > 0$ 时，若 $h_n = o(n^{\alpha})$ 且 $\alpha < \xi$，则 CLT 成立；当 $\chi = 0$ 时，若 $\alpha < 1/(d-1)$，则 CLT 成立，该结论基于矩界和方差控制。
• 作者证明了 $\left\| \frac{a_n(h_n) - \mathbb{E}[a_n(h_n)]}{\sqrt{\operatorname{Var}(a_n(h_n))}} - \frac{\sum_{i=1}^{2m+1}(X_i - \mathbb{E}[X_i])}{\sigma_n(h_n)} \right\|_2 \to 0$ 当 $n \to \infty$，确认在 $L^2$ 范数下收敛于正态分布。