[논문 리뷰] Central limit theorem over non-linear functionals of empirical measures with applications to the mean-field fluctuation of interacting particle systems

Benjamin Jourdain, Alvin Tse|arXiv (Cornell University)|2020. 02. 04.

Stochastic processes and statistical mechanics참고 문헌 31인용 수 3

한 줄 요약

이 논문은 선형 기능도수의 정규성 조건 하에서 비선형 기능도수의 경험 측도에 대한 중심극한정리(CLT)를 수립하여, 상호작용 확산 시스템에서 평균장 변동을 분석할 수 있게 한다. 변동을 마팅게일 타입의 합과 하드러르 연속성에 의해 통제되는 나머지 항으로 분해함으로써, 저자들은 가우시안 극한으로의 약한 수렴을 증명한다. 이는 입자 시스템과 McKean-Vlasov SDE에 응용된다.

ABSTRACT

In this work, a generalised version of the central limit theorem is proposed for nonlinear functionals of the empirical measure of i.i.d. random variables, provided that the functional satisfies some regularity assumptions for the associated linear functional derivative. This generalisation can be applied to Monte-Carlo methods, even when there is a nonlinear dependence on the measure component. We use this result to deal with the contribution of the initialisation in the convergence of the fluctuations between the empirical measure of interacting diffusion and their mean-field limiting measure (as the number of particles goes to infinity), when the dependence on measure is nonlinear. A complementary contribution related to the time evolution is treated using the master equation, a parabolic PDE involving L-derivatives with respect to the measure component, which is a stronger notion of derivative that is nonetheless related to the linear functional derivative.

연구 동기 및 목표

• 비독립 identical 분포(이.I.D.) 설정을 초월하여 경험 측도의 비선형 기능도수에 대한 중심극한정리를 확장하는 것, 특히 기능도수가 측도에 비선형적으로 의존할 경우에 대한 것.
• 상호작용 입자 시스템의 평균장 변동에서 초기화 효과를 해결하기 위해, 극한 분산이 초기 측도 구성에 어떻게 의존하는지 분석함으로써 이 효과를 정량화하는 것.
• 시간에 따른 변동의 진화를 기술하기 위해 선형 기능도수와 L-도수(주로 마스터 방정식을 통해)를 사용하는 프레임워크를 개발하는 것.
• 변동 과정이 가우시안 과정으로 약한 수렴함을 확립하여, 초기화와 시간 진동의 기여를 분리하는 것.
• 특히 고차원 또는 복잡한 확률 시스템에서 경험 측도에 비선형적으로 의존하는 몬테카를로 방법에 대한 엄밀한 확률적 기반을 제공하는 것.

제안 방법

• 경험 측도 $ \frac{1}{N}\sum_{i=1}^N \delta_{\zeta_i} $ 의 비선형 기능도수 $ U $ 에 대해 일반화된 CLT를 제안하며, 선형 기능도수 $ \frac{\delta U}{\delta m}(m, y) $ 가 존재하고 정규성을 만족한다고 가정한다.
• 입자 $ i $ 에 따라 인덱싱된 마팅게일 증분의 합으로 $ \sqrt{N}(U(\mu_N) - U(m_0)) $ 를 새로운 방식으로 분해하며, 과거 입자들에 의존하는 랜덤 측도의 매개변수 $ m_{N,i}^s $ 를 사용한다.
• 총변동에 대한 나머지 항을 총변동에 대한 측도 매개변수에 대한 하드러르 연속성으로 통제하여, 확률적으로 0으로 수렴함을 보장한다.
• 마팅게일 CLT를 검증하기 위해 린데베르그 조건과 브라켓 수렴을 적용하여 渐近 정규성을 확보한다.
• 시간 진동의 과정을 모델링하기 위해 마스터 방정식(선형 기능도수보다 강력한 도수인 L-도수를 포함하는 포물형 PDE)을 사용한다.
• 모멘트 추정, 커플링 방법, 레비의 연속성 정리의 조합을 통해 유한차원 분포의 유계성과 약한 수렴을 확립한다.

실험 결과

연구 질문

• RQ1 경험 측도의 비선형 기능도수 $ U $ 에 대해 $ \sqrt{N}(U(\mu_N) - U(m_0)) $ 가 가우시안 분포로 수렴하는 조건은 무엇인가?
• RQ2 초기 설정이 평균장 변동에 미치는 기여를 어떻게 엄밀하게 정량화하고 시간 진동 성분과 분리할 수 있는가?
• RQ3 선형 기능도수 $ \frac{\delta U}{\delta m} $ 의 하드러르 연속성이 CLT에서 나머지 항이 사라지도록 보장하는 데 어떤 역할을 하는가?
• RQ4 L-도수를 포함하는 마스터 방정식은 상호작용 확산 시스템에서 변동의 시간 진동을 어떻게 기술하는가?
• RQ5 한계 변동 과정이 초기화에서 비롯하는 독립적인 성분과 시간에 의존하는 동역학에서 비롯하는 성분으로 분해될 수 있는가?

주요 결과

• 논문은 선형 기능도수의 미약한 정규성 조건 하에서 경험 측도의 비선형 기능도수에 대한 중심극한정리를 수립하여 프레셰 도수의 필요성을 피한다.
• 선형 기능도수 $ \frac{\delta U}{\delta m} $ 가 총변동에 대해 지수 $ \alpha > 1/2 $ 를 가진 하드러르 연속성을 만족할 경우, $ \sqrt{N}(U(\mu_N) - U(m_0)) $ 전개의 나머지 항이 $ N \to \infty $ 일 때 확률적으로 0으로 수렴한다.
• 주요 변동 항은 초기 측도에서 평가된 기능도수의 공분산에 의해 결정되는 중심가우시안 랜덤 변수로 약한 수렴을 보인다.
• 시간 진동 성분은 마스터 방정식의 해와 분산 계수 $ \sigma $ 를 포함한 공분산構조를 가진 가우시안 과정으로 수렴한다.
• 한계 변동 과정은 두 개의 독립적인 가우시안 성분의 합으로 특징지어지며, 하나는 시간에 무관한 초기화에서 기인하고, 다른 하나는 McKean-Vlasov 과정의 경로에 의존하는 시간 진동에서 기인한다.
• 유한차원 분포의 변동 과정은 특성함수 수렴과 레비의 연속성 정리를 통해 초기 기여와 동적 기여를 조합한 명시적 공분산 행렬을 가진 다변량 정규분포로 약한 수렴을 보인다.