[논문 리뷰] Simulation of particle systems interacting through hitting times

Vadim Kaushansky, Christoph Reisinger|arXiv (Cornell University)|2018. 05. 29.

Markov Chains and Monte Carlo Methods참고 문헌 10인용 수 12

한 줄 요약

이 논문은 입자가 경계에 도달할 때의 도달 시간을 통해 평균장 상호작용을 갖는 McKean–Vlasov 방정식을 시뮬레이션하기 위한 새로운 오일러형 입자 방법을 제안한다. 입자가 경계에서 흡수될 때 피드백이 발생한다. 시간 간격에 대해 수렴 속도 1/2를 확립하였으며, 브라운 운동 다리와 국소 메시 조절을 통해 수렴 속도를 1로 향상시켰다. 이는 수치적으로 검증되었으며, 붕괴 행동을 보이는 시스템적 리스크 및 신경망 모델에 적용되었다.

ABSTRACT

We develop an Euler-type particle method for the simulation of a McKean--Vlasov equation arising from a mean-field model with positive feedback from hitting a boundary. Under assumptions on the parameters which ensure differentiable solutions, we establish convergence of order $1/2$ in the time step. Moreover, we give a modification of the scheme using Brownian bridges and local mesh refinement, which improves the order to $1$. We confirm our theoretical results with numerical tests and empirically investigate cases with blow-up.

연구 동기 및 목표

• 경계에 도달함으로써 발생하는 양의 피드백을 갖는 McKean–Vlasov 방정식에 대한 체계적인 수치적 방법을 개발하는 것. 이는 이전에 체계적인 시뮬레이션 방법이 없는 방정식의 클래스이다.
• 해당 방정식의 입자 기반 시간 스텝 스킴의 수렴 속도를, 해가 미분 가능함을 보장하는 가정 하에 분석하는 것.
• 브라운 운동 다리 보간법과 국소 메시 조절을 사용하여 수렴 차수를 1/2에서 1로 향상시키는 것.
• 시스템적 리스크 및 신경망 모델에서의 붕괴 현상에 대한 경험적 조사 수행.
• 대규모 입자 시스템을 사용한 수치 실험을 통해 이론적 수렴 속도를 검증하는 것.

제안 방법

• 경계에 도달하는 입자의 도달 시간 밀도에 따라 변동성이 변화하는 McKean–Vlasov SDE를 시뮬레이션하기 위해 시간 이산화된 오일러형 입자 스킴을 제안한다.
• 입자 경로로부터 도달 시간 분포 $ L_t $ 의 경험적 추정치를 사용하여 피드백 드리프트를 계산한다.
• 흡수 경계 근처의 경로 근사 향상을 위해 수정된 스킴이 브라운 운동 다리를 사용한다.
• 도달 시간 근처에 적응적으로 국소 메시 조절을 적용하여 중요 영역의 정확도를 향상시킨다.
• 초기 법칙과 매개변수에 대한 규칙성 가정 하에 수렴을 분석하며, 붕괴 시간까지 해가 미분 가능함을 보장한다.
• 분산 감소 및 교환 가능성 논증을 사용하여 경험적 도달 시간 측도의 수렴성을 확률론적 및 $ L^2 $-노름에서 입증한다.

실험 결과

연구 질문

• RQ1 도달 시간 피드백이 있는 McKean–Vlasov 방정식에 대한 표준 오일러형 입자 스킴의 수렴 속도는 무엇인가?
• RQ2 브라운 운동 다리 보간법과 국소 메시 조절을 통해 수렴 속도를 1/2를 초월해 향상시킬 수 있는가?
• RQ3 도달 시간 분포에 비연속성이 있거나 붕괴가 발생하는 경우 스킴은 어떻게 행동하는가?
• RQ4 입자 시스템 크기가 경험적 도달 시간 추정치의 정확도에 어떤 영향을 미치는가?
• RQ5 이 방법은 금융 네트워크의 시스템적 리스크나 신경집단 역학과 같은 실제 모델에 적용될 수 있는가?

주요 결과

• 정규성 가정 하에 표준 오일러형 입자 스킴은 시간 간격에 대해 수렴 속도 1/2를 확보한다.
• 브라운 운동 다리와 국소 메시 조절을 사용하는 수정된 스킴은 수렴 속도를 1로 향상시킨다.
• 수치 실험은 이론적 수렴 속도를 확인하며 붕괴 행동을 잘 포착하는 데에 뛰어난 안정성을 보였다.
• $ N \to \infty $ 일 때 경험적 도달 시간 측도 $ \tilde{L}_t $ 가 진짜 도달 시간 분포로 확률적으로 수렴한다.
• $ N \to \infty $ 일 때 경험적 도달 시간 추정자 $ \hat{L}_t^N $ 의 분산은 0으로 수렴하여 일관성을 보장한다.
• 이 방법은 재무 및 신경망 응용에서 관찰되는 바와 같이, $ L_t $ 의 비연속성 형성도 잘 포착한다.