[논문 리뷰] Stein's method for dependent random variables occurring in Statistical Mechanics

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

Random Matrices and Applications참고 문헌 11인용 수 10

한 줄 요약

이 논문은 통계역학의 큐리-바이스 모델에서 종속된 랜덤 변수의 부분합에 대해 최적의 베르누이-에세른 경계를 확립하기 위해 교환가능한 쌍에 대한 스텐의 방법을 개발한다. 비정규 근사 분포에 비례하는 밀도를 가진 비정규 근사 분포로의 수렴 속도를 증명하며, 특히 임계 온도 $\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.

연구 동기 및 목표

• 통계역학에서 발생하는 비정규 근사 분포에 대해 스텐의 방법을 교환가능한 쌍으로 확장하기 위해.
• 큐리-바이스 모델에서 종속적이고 동일하게 분포된 랜덤 변수의 부분합에 대해 수렴 속도의 정밀한 경계를 설정하기 위해.
• 역온도 $\beta_n \to 1$이 되는 경우, 즉 임계 온도로 수렴할 때 수렴 속도의 행동을 분석하기 위해.
• 그리피스-허스트-셔먼(GHS) 부등식을 만족하는 모델에 대해 베르누이-에세른 경계의 엄밀한 프레임워크를 제공하기 위해.
• 표준 중심극한정리가 실패하는 임계점에서조차도 총 자화율에 대한 중심극한정리의 최적 수렴 속도를 유도하기 위해.

제안 방법

• 밀도가 $\exp(-\mu|x|^{2k}/(2k)!)$ 비례하는 비정규 목표 분포의 클래스에 대해 교환가능한 쌍에 대한 스텐의 방법을 적응하여 적용한다.
• 밀도 $p(x) \propto \exp(-a_k V(x))$를 가진 분포에 대해 스텐의 항등식과 스텐 방정식의 해를 사용한다. 여기서 $V(x) = |x|^{2k}$이다.
• 부분적분과 모멘트 추정을 통해 스텐 방정식의 해 $f_h$에 대한 경계를 확립하여, 방법에 필요한 정규성 조건을 보장한다.
• 스핀이 $P_{n,\beta,h}$에 따라 분포하는 큐리-바이스 모델에 이 방법을 적용하며, 총 자화율 $S_n$이 핵심 합이다.
• Gibbs 샘플링 또는 국소 이동 동역학을 통해 구성된 교환가능한 쌍 $(X,Y)$에 대해 $\mathbb{E}|f_h(X) - f_h(Y)|$에 대한 경계를 도출한다.
• GHS 부등식을 사용하여 잠재 에너지 함수 $V(x)$의 필요한 볼록성 및 대칭성 성질을 확보함으로써, 모멘트 및 尾 확률 추정을 가능하게 한다.

실험 결과

연구 질문

• RQ1 큐리-바이스 모델에서 임계 온도 $\beta_c = 1$에서 총 자화율에 대한 중심극한정리의 최적 수렴 속도는 무엇인가?
• RQ2 역온도 $\beta_n \to 1$이 위 또는 아래에서 수렴할 때 수렴 속도는 어떻게 행동하는가?
• RQ3 스텐의 방법은 첨도가 높거나 고차 다항 꼬리가 있는 비정규 근사 분포, 예를 들어 $\exp(-\mu|x|^{2k}/(2k)!)$에 대해 확장될 수 있는가?
• RQ4 평균장 모델에서 종속된 랜덤 변수에 대해 고전적 중심극한정리가 실패하는 조건은 무엇인가?
• RQ5 큐리-바이스 모델에서 종속된 수열에 대해 스텐의 방법을 적용할 수 있도록 교환가능한 쌍을 어떻게 구성할 수 있는가?

주요 결과

• 논문은 고전적 중심극한정리가 실패하는 임계 온도 $\beta_c = 1$에서 큐리-바이스 모델의 총 자화율에 대해 최적의 베르누이-에세른 수렴 속도 $O(n^{-1/2})$를 확립한다.
• 모든 $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 부등식을 만족하는 모델에 적용 가능하며, 잠재 에너지 함수에 필요한 볼록성 및 대칭성을 보장한다.
• 비정규 목표 분포에 대한 스텐 방정식의 해 $f_h$는 $f_h$, $f_h'$, $f_h''$에 대해 균일한 경계를 만족하며, 이는 스텐 방법에서 오차를 유한하게 유지하는 데 필수적이다.