[论文解读] Perfect simulation and finitary coding for multicolor systems with interactions of infinite range

Antonio Galves, Nancy L. Garcia+1 more|ArXiv.org|Sep 20, 2008

Markov Chains and Monte Carlo Methods参考文献 11被引用 7

一句话总结

本文提出了一种针对 ℤᵈ 上具有无限范围相互作用的多色相互作用粒子系统的完美模拟算法，通过将无限范围的转移率表示为可数个有限范围分量的混合，实现了从 i.i.d. 过程到系统不变测度的有限编码，建立了指数遍历性，并将先前结果推广至非马尔可夫性、非齐次性及长程模型，且无需单调性或对偶性假设。

ABSTRACT

We consider a particle system on $Z^d$ with finite state space and interactions of infinite range. Assuming that the rate of change is continuous and decays sufficiently fast, we introduce a perfect simulation algorithm for the stationary process. The algorithm follows from a representation of the multicolor system as a finitary coding from a sequence of independent uniform random variables. This implies that the process is exponentially ergodic. The basic tool we use is a representation of the infinite range change rates as a mixture of finite range change rates.

研究动机与目标

• 为 ℤᵈ 上具有无限范围、连续且快速衰减的转移率的多色系统开发一种完美模拟算法。
• 建立从有限值 i.i.d. 过程到此类系统不变测度的有限编码的存在性。
• 将已知的有限编码与完美模拟结果推广至非马尔可夫性、非有限范围或非单调模型之外。
• 在转移率的温和正则性条件下，证明该过程的指数遍历性。
• 通过向后与向前模拟步骤，提供一个自包含的平稳过程构造，避免依赖对偶性或吸引性假设。

提出的方法

• 通过定理 1 将无限范围的转移率 $ c_i(a,\eta) $ 表示为可数个有限范围率的混合，实现局部近似。
• 构建一个向后黑白草图算法，以向后时间模拟过程历史，由于几何尾部界的存在而几乎必然停止。
• 使用向前算法，基于混合表示与条件率，为黑白草图分配实际颜色。
• 通过耦合论证与次鞅不等式，建立几乎必然的有限停止时间，依赖于长程影响的衰减。
• 通过有界停止时间 $ V_{\text{STOP}} $ 的矩，证明指数遍历性，得到 $ \mathbb{E}[V_{\text{STOP}}] \leq \frac{1}{1 - \bar{\lambda}} $。
• 通过定义一个从 i.i.d. 输入的平移不变、可测映射，构造有限编码，使得原点处的输出仅依赖于输入的有限随机子集，利用完美模拟输出。

实验结果

研究问题

• RQ1 在转移率的连续性与衰减性条件较弱时，能否为具有无限范围相互作用的多色系统构造完美模拟算法？
• RQ2 此类系统是否即使在无马尔可夫性或单调性假设下，也能从有限 i.i.d. 过程实现到其不变测度的有限编码？
• RQ3 能否在不依赖对偶性或吸引性假设的前提下，构造性地证明平稳过程的存在性？
• RQ4 在无限范围相互作用粒子系统中，何种条件可确保指数遍历性？
• RQ5 如何将无限范围率结构分解为可数个有限范围分量的混合，以实现模拟？

主要发现

• 完美模拟算法几乎必然在有限时间内停止，且 $ \mathbb{E}[V_{\text{STOP}}] \leq \frac{1}{1 - \bar{\lambda}} $，确保计算成本有限。
• 通过停止时间的矩有界性，建立了指数遍历性，意味着快速收敛至不变测度。
• 存在从有限 i.i.d. 过程到系统不变测度的有限编码，如定理 4 所示，将结果推广至非马尔可夫与无限范围模型。
• 该构造推广了 Van den Berg 与 Steif (1999) 及 Häggström 与 Steif (2000) 的结果，现可适用于非马尔可夫性、长程吉布斯测度。
• 对于热浴动力学，条件 $ \alpha(-1) > \frac{2d}{2d+1} $ 确保停止时间具有有限期望，其中 $ \alpha(-1) = \sum_{a \in A} \min_{\zeta(\partial 0) \in A^{\partial 0}} \mathcal{Q}(X(0)=a \mid X(\partial 0)=\zeta(\partial 0)) $。
• 将无限范围率表示为有限范围率混合（定理 1）是实现完美模拟与有限编码的关键技术工具。