[논문 리뷰] Binary branching processes with Moran type interactions

Alexander M. G. Cox, E Horton외 1명|arXiv (Cornell University)|2022. 07. 07.

Stochastic processes and statistical mechanics인용 수 4

한 줄 요약

이 논문은 모리안형 상호작용을 갖는 새로운 이진 분열 입자 시스템을 제안하며, 입자 사망 시 복제가 발생하고 분열 시 다른 입자의 사망이 유도되어 크기 제약이 있는 인구구성 동역학을 가능하게 한다. 주요 기여는 입자 시스템의 점유 측도와 기저 마코프 과정 사이의 명시적 페인만-카크 유형의 이중성으로, 입자 수가 증가함에 따라 실제 과정으로의 L² 수렴 속도를 확보하며, 고정된 크기의 모리안 모델보다 타깃 생식-사망 과정을 더 잘 근사함을 보여준다.

ABSTRACT

The aim of this paper is to study the large population limit of a binary branching particle system with Moran type interactions: we introduce a new model where particles evolve, reproduce and die independently and, with a probability that may depend on the configuration of the whole system, the death of a particle may trigger the reproduction of another particle, while a branching event may trigger the death of an other one. We study the occupation measure of the new model, explicitly relating it to the Feynman-Kac semigroup of the underlying Markov evolution and quantifying the L2 distance between their normalisations. This model extends the fixed size Moran type interacting particle system discussed in [18, 19, 6, 7, 57] and we will indeed show that our model outperforms the latter when used to approximate a birth and death process. We discuss several other applications of our model including the neutron transport equation [36, 15] and population size dynamics.

연구 동기 및 목표

• 시간에 따라 변화하는 생식-사망 과정을 더 잘 근사하기 위해 이진 분열 과정과 모리안형 재표본 추출을 통합한 통합 입자 시스템 모델을 개발한다.
• 입자 시스템의 경험 측도와 기저 마코프 과정의 피네만-카크 준거군 사이의 정밀한 수학적 관계를 확립한다.
• 입자 시스템의 정규화된 점유 측도와 정규화된 피네만-카크 준거군 사이의 L² 거리를 정량화하고, N → ∞일 때 수렴함을 보인다.
• 고정된 크기의 모리안 모델을 일반화하여 동적 인구구성 크기와 구성에 의존하는 선택/재표본 추출 규칙을 허용한다.
• 특히 인구구성 크기 동역학과 중성자 이동 문제에서, 기존의 고정된 크기의 모리안 시스템보다 모델의 우수성을 입증한다.

제안 방법

• 모델은 각 입자가 마코프 과정 X에 따라 독립적으로 진화하는 이진 분열 과정을 포함하며, 분열 및 사망 사건은 각각 비율 b와 κ에서 발생한다.
• 분열 또는 사망 사건이 발생할 경우, 구성에 따라 재표본 추출 또는 선택 규칙을 적용한다: 사망 시 복제가 발생하고, 분열 시 다른 입자의 사망이 유도된다.
• 시간에 따라 변하는 가중치 Π_T^A와 Π_T^B를 정의하여 재표본 추출 및 선택 사건의 영향을 보정하며, 이를 통해 시스템의 경험 측도에 미치는 영향을 명시적으로 기술한다.
• 핵심적인 이중성은 다음과 같다: 시간 T에서 가중치가 부여된 경험 측도의 기대값은 누적 분열 및 사망 비율의 지수 함수로 가중된 f(X_T)의 기대값과 정확히 일치한다.
• 입자 시스템의 정규화된 점유 측도와 정규화된 피네만-카크 준거군 사이의 L² 거리는 유계이며, 초기 입자 수 N → ∞일 때 0으로 수렴함을 보였다.
• 모델은 N_min에서 N_max 사이의 제약된 인구구성 크기 영역에 적용되며, 경계에서 선택/재표본 추출 확률을 1로 설정하여 크기 제어를 달성한다.

실험 결과

연구 질문

• RQ1 구성에 따라 재표본 추출이 이루어지는 이진 분열 입자 시스템은 어떻게 확장되어야 하며, 이를 통해 인구구성 크기를 유지하면서도 목표 과정으로의 수렴을 보장할 수 있는가?
• RQ2 이러한 입자 시스템의 경험 측도와 기저 마코프 과정의 피네만-카크 준거군 사이의 정밀한 수학적 관계는 무엇인가?
• RQ3 입자 시스템의 정규화된 점유 측도와 정규화된 피네만-카크 준거군 사이의 L² 거리는 명시적으로 유계이며, N → ∞일 때 0으로 수렴하는가?
• RQ4 이 모델은 고정된 크기의 기존 모리안 모델보다 생식-사망 과정을 근사하는 데서 우월한가, 특히 수렴 속도와 편향 측면에서?
• RQ5 입자 시스템은 기저 과정 X로부터 강한 마코프 성질을 어떤 조건에서 상속하는가?

주요 결과

• 시간 T에서 입자 시스템의 가중 경험 측도의 기대값은 b와 κ를 포함하는 시간에 따라 변하는 지수 함수로 가중된 f(X_T)의 피네만-카크 기대값과 정확히 일치한다.
• 입자 수 N → ∞일 때, 입자 시스템의 정규화된 점유 측도와 정규화된 피네만-카크 준거군 사이의 L² 거리는 0으로 수렴한다.
• 적절한 매개변수를 선택함으로써, 이 모델은 고정된 크기의 모리안 과정에 대해 기존의 다수의 한 명식과 비편향 추정량 [57]을 특수한 경우로 포함한다.
• 제약된 N_min–N_max 모델은 경계에서 선택/재표본 추출 확률를 1로 설정함으로써 인구구성 크기가 항상 범위 내에 유지되도록 보장하며, 크기 제약이 있는 과정의 안정적인 시뮬레이션을 가능하게 한다.
• 특히 비정상적인 인구구성 크기 또는 시간에 따라 변화하는 동역학을 갖는 시나리오에서, 이 모델은 기존의 고정된 크기의 모리안 시스템보다 생식-사망 과정을 더 잘 근사한다.
• 입자 시스템의 법칙은 기저 과정에 의해 생성된 필터레이션에 대해 마코프 성질을 만족하며, 정규성 조건이 만족될 경우 시스템은 강한 마코프 성질을 상속한다.