[論文レビュー] Binary branching processes with Moran type interactions

Alexander M. G. Cox, E Horton+1 more|arXiv (Cornell University)|Jul 7, 2022

Stochastic processes and statistical mechanics被引用数 4

ひとこと要約

本稿では、粒子の死が再生を引き起こし、分岐が他の粒子の死を引き起こす、モーラン型の相互作用を持つ新規な二値分岐粒子系を提案する。これにより、サイズ制約付きの集団動態が可能となる。主な貢献は、粒子系の占有測度と基礎となるマコフ過程との間で、Feynman-Kac型の双対性を明示的に確立したことである。集団サイズが増加するにつれ、真の過程への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.

研究の動機と目的

• 時間非定常な出生・死過程をよりよく近似できる、二値分岐過程とモーラン型リサンプリングを統合した統一的な粒子系モデルの構築を目的とする。
• 粒子系の経験的測度と基礎となるマコフ過程のFeynman-Kac半群との間の明確な数学的関係を確立すること。
• 粒子系の正規化された占有測度と正規化されたFeynman-Kac半群との間のL²距離を定量的に評価し、N → ∞ の際に収束することを示すこと。
• 固定サイズのモーランモデルを一般化し、動的変動する集団サイズと、状態依存の選択・リサンプリングルールを許容すること。
• 特に集団サイズの動態と中性子輸送の文脈において、古典的固定サイズのモーラン系に比べて優れた性能を示すことを実証すること。

提案手法

• モデルは、各粒子がマコフ過程 X に従って独立に発展する二値分岐過程を導入し、分岐率 b と死滅率 κ でイベントが発生する。
• 分岐または死滅イベントが発生すると、状態依存のリサンプリングまたは選択ルールが適用される：死が発生すると再生が引き起こされ、分岐が発生すると他の粒子の死が引き起こされる。
• リサンプリングおよび選択イベントの影響を明示的に特定するため、時間依存の重み Π_T^A および Π_T^B を定義し、それらが系の経験的測度に与える影響を特徴づける。
• 重要な双対性が確立される：時間 T における重み付き経験的測度の期待値は、累積的な分岐率および死滅率の指数関数的関数で重み付けられた f(X_T) の期待値に等しい。
• 粒子系の正規化された占有測度と正規化されたFeynman-Kac半群との間のL²距離が有界であり、初期粒子数 N → ∞ の際に0に収束することを示した。
• 集団サイズが N_min から N_max の範囲に制限される制約付きの状態において、境界で選択／リサンプリング確率を1に設定することで、サイズ制御を実現した。

実験結果

リサーチクエスチョン

• RQ1 状態依存のリサンプリングを組み込んだ二値分岐粒子系は、集団サイズを維持しつつ、目標過程への収束を保つようにどのように拡張可能か？
• RQ2 このような粒子系の経験的測度と基礎となるマコフ過程のFeynman-Kac半群との間の明確な数学的関係は何か？
• RQ3 粒子系の正規化された占有測度と正規化されたFeynman-Kac半群との間のL²距離を明示的に有界化でき、N → ∞ の際に0に収束することが示せるか？
• RQ4 このモデルは、特に収束速度とバイアスの観点から、古典的固定サイズのモーラン系に比べて出生・死過程の近似において優れているか？
• RQ5 粒子系は、基礎過程 X が持つ強いマルコフ性をどのような条件下で継承するか？

主な発見

• 時間 T における粒子系の重み付き経験的測度の期待値は、b および κ を含む時間依存の指数関数的重みを伴う f(X_T) のFeynman-Kac期待値に正確に等しい。
• 粒子系の正規化された占有測度と正規化されたFeynman-Kac半群との間のL²距離は、初期粒子数 N → ∞ の際に0に収束する。
• 適切なパrameter選択により、このモデルは固定サイズのモーラン過程における古典的な「多数から1つへの公式」と、[57]における不偏推定量の両方を特別な場合として回復する。
• N_min–N_max の制約付きモデルでは、境界で選択／リサンプリング確率を1に設定することで、集団サイズが境界内に保たれ、サイズ制約付き過程の安定なシミュレーションが可能となる。
• 特に非定常な集団サイズや時間非定常なダイナミクスを示す状況において、出生・死過程の近似において、古典的固定サイズのモーラン系に比べて優れた性能を示す。
• 粒子系の法則は、基礎過程が生成するフィルトレーションに関してマルコフ性を満たし、正則性条件の下では、強いマルコフ性を継承する。