[论文解读] Stein's method for dependent random variables occurring in Statistical Mechanics

Peter Eichelsbacher, Matthias Löwe|ArXiv.org|Aug 13, 2009

Random Matrices and Applications参考文献 11被引用 10

一句话总结

本文针对统计力学中Curie-Weiss模型的相依随机变量部分和，发展了交换对的Stein方法，建立了最优Berry-Esseen界。研究证明了收敛速率趋近于非高斯极限分布，其密度与$\exp(-\mu|x|^{2k}/(2k)!)$成正比，包括在临界温度$\beta_c = 1$时中心极限定理的最优$O(n^{-1/2})$速率，即使标准中心极限定理在此失效。

ABSTRACT

We obtain rates of convergence in limit theorems of partial sums $S_n$ for certain sequences of dependent, identically distributed random variables, which arise naturally in statistical mechanics, in particular, in the context of the Curie-Weiss models. Under appropriate assumptions there exists a real number $α$, a positive real number $μ$, and a positive integer $k$ such that $(S_n- n α)/n^{1 - 1/2k}$ converges weakly to a random variable with density proportional to $\exp(-μ|x|^{2k} /(2k)!)$. We develop Stein's method for exchangeable pairs for a rich class of distributional approximations including the Gaussian distributions as well as the non-Gaussian limit distributions with density proportional to $\exp(-μ|x|^{2k} /(2k)!)$. Our results include the optimal Berry-Esseen rate in the Central Limit Theorem for the total magnetization in the classical Curie-Weiss model, for high temperatures as well as at the critical temperature $β_c=1$, where the Central Limit Theorem fails. Moreover, we analyze Berry-Esseen bounds as the temperature $1/ β_n$ converges to one and obtain a threshold for the speed of this convergence. Single spin distributions satisfying the Griffiths-Hurst-Sherman (GHS) inequality like models of liquid helium or continuous Curie-Weiss models are considered.

研究动机与目标

• 将Stein方法推广至交换对，以处理统计力学中出现的非高斯极限分布。
• 为Curie-Weiss模型中相依同分布随机变量的部分和建立精确的收敛速率。
• 分析当逆温度$\beta_n \to 1$（即临界温度）时收敛速率的行为。
• 为满足Griffiths-Hurst-Sherman（GHS）不等式的模型提供Berry-Esseen界的严格框架。
• 在总磁化量的中心极限定理中，即使在标准CLT失效的临界点，也推导出最优收敛速率。

提出的方法

• 将Stein方法用于交换对，适用于密度与$\exp(-\mu|x|^{2k}/(2k)!)$成正比的一类非高斯目标分布。
• 利用具有密度$p(x) \propto \exp(-a_k V(x))$的分布的Stein恒等式及Stein方程的解，其中$V(x) = |x|^{2k}$。
• 通过分部积分和矩估计建立Stein方程解$f_h$的界，确保方法所需的正则性条件。
• 将该方法应用于自旋服从$P_{n,\beta,h}$分布的Curie-Weiss模型，其中总磁化量$S_n$是关键的和。
• 为通过Gibbs抽样或局部移动动力学构造的交换对$(X,Y)$，推导$\mathbb{E}|f_h(X) - f_h(Y)|$的界。
• 利用GHS不等式确保势函数$V(x)$所需的凸性和对称性，从而支持矩估计和尾部估计。

实验结果

研究问题

• RQ1 在临界温度$\beta_c = 1$时，Curie-Weiss模型中总磁化量的中心极限定理的最优收敛速率是多少？
• RQ2 当逆温度$\beta_n$从上方或下方趋近于1时，收敛速率如何变化？
• RQ3 能否将Stein方法扩展至具有重尾或高阶多项式尾部（如$\exp(-\mu|x|^{2k}/(2k)!)$）的非高斯极限分布？
• RQ4 在平均场模型中，相依随机变量的古典中心极限定理在何种条件下会失效？
• RQ5 在Curie-Weiss模型中，如何构造交换对以支持对相依序列应用Stein方法？

主要发现

• 本文在临界温度$\beta_c = 1$时，为Curie-Weiss模型中的总磁化量建立了最优Berry-Esseen速率$O(n^{-1/2})$，此时经典CLT失效。
• 当$k \geq 2$时，归一化和$(S_n - n\alpha)/n^{1-1/2k}$弱收敛于密度与$\exp(-\mu|x|^{2k}/(2k)!)$成正比的非高斯极限，收敛速率为$O(n^{-1/2k})$。
• 该收敛速率被证明是最优的，即在给定的依赖结构下不可能存在更快的收敛速率。
• 当$\beta_n \to 1$时，收敛至非高斯极限的速率由一个依赖于$\beta_n$趋近1速度的阈值控制，且当$\beta_n$从上方趋近1时，速率逐渐恶化。
• 该方法适用于满足GHS不等式的模型，如连续Curie-Weiss模型和液氦模型，确保势函数所需的凸性和对称性。
• 针对非高斯目标分布的Stein方程解$f_h$满足$f_h$、$f_h'$和$f_h''$的统一有界性，这对控制Stein方法中的误差至关重要。