Branching-selection particle systems and inverse first passage problems

技術要約：分岐選択粒子系と逆第一通過問題

問題設定
本論文は、**一般化された逆第一通過問題（generalized inverse first passage problem: IFPP）**を扱う。古典的なIFPPは、Shiryaevによる問いに端を発しており、標準的なブラウン運動 $(W_t)_{t \ge 0}$ と $[0, \infty]$ 上の目標確率分布 $p$ が与えられたとき、境界 $b(t)$ が存在して、第一通過時刻 $\tau := \inf\{t : W_t \le b(t)\}$ が分布 $p$ に従うかどうかを問うものである。

本論文は、**ソフトキル型逆第一通過問題（soft-killed inverse first passage problem: SKIFPP）**として知られる一般化された問題に焦点を当てている。ここでは、停止時刻は次のように定義される：
$$\tau := \inf \left\{ t : \Lambda \int_0^t \omega(W_s - b(s)) ds \ge U \right\}$$
ここで、$U \sim \text{Exp}(1)$ であり、$\Lambda \in (0, \infty)$、そして $\omega: \mathbb{R} \to [0, 1]$ は単調減少関数である。境界を越えた瞬間にデフォルトが発生する古典的なケースとは異なり、この定式化では、デフォルトは経路の境界に対する履歴の累積関数としてモデル化される。目的は、特定の初期分布 $W_0$ および目標生存確率 $p([t, \infty))$ が与えられたときに、境界 $b(t)$ を構成することである。

手法
著者は、SKIFPPの解を近似するために、分岐選択粒子系（具体的には、$N$-分岐ブラウン運動（$N$-BBM）の変種）を構築する。手法は以下の通りである：

1. 粒子系の構築:

   • $N$ 個の独立なブラウン運動として動く粒子を考える。
   • 粒子は、時間依存のレート $f(t)$ で分岐する。
   • 分岐と同時に、粒子が一つ削除される。この削除は、選択重み関数 $\omega$ に基づいて「ソフトキル」される。具体的には、分岐時において、位置 $x$ にある粒子は、粒子の経験分布によって決定される移動境界 $b_N(t)$ を用いた確率 $\omega(x - b_N(t))$ に比例して削除される。
   • 境界 $b_N(t)$ は、総分岐レートが規定されたレート $Nf(t)$ と一致するように、以下のように暗黙的に定義される：
     $$b_N(t) := \inf \left\{ x : \Lambda \langle \mu_N^t, \omega(\cdot - x) \rangle \ge f(t) \right\}$$
     ここで $\mu_N^t$ は $N$ 個の粒子の経験測度である。

2. 流体力学的極限解析:

   • $N \to \infty$ における経験測度 $\mu_N^t$ の極限を解析する。
   • マルチンゲール法（具体的には伊藤の公式とブルカ・デ・ギディ・ベンゲリ（BDG）不等式）を用いて、経験測度のテスト関数の進化方程式を導出する。
   • スキロコフ（Skorokhod）位相を備えた測度値過程の空間における、法則のタイトネス（緊密性）を確立する。
   • 極限が、偏微分方程式（PDE）によって支配される**自由境界問題（FBP）**を満たすことが示される。

主要な結果

1. 自由境界問題への収束（定理1）:
   構築された $(\omega, f, \Lambda, N)$-BBM の経験測度 $\mu_N^t$ は、決定論的な測度 $\mu_t$（密度 $u(x, t)$ を持つ）に弱収束する。この密度は、以下の自由境界問題を解く唯一の弱解である：
   $$\begin{cases} u_t = \frac{1}{2}u_{xx} + f(t)u - \Lambda \omega(x - b(t))u, & x \in \mathbb{R}, t > 0 \\ \int_{\mathbb{R}} u(y, t) dy = 1, & t > 0 \\ u(x, 0) = \rho(x), & x \in \mathbb{R} \end{cases}$$
   ここで、自由境界 $b(t)$ は積分制約 $\int_{\mathbb{R}} \omega(x - b(t)) u(x, t) dx = f(t)/\Lambda$ によって決定される。

2. SKIFPPとの関連性（定理2）:
   上記の自由境界問題 $(u, b)$ を解が、SKIFPP の一意な解となることを本論文は示している。具体的には、$\tau$ が $\omega$ と $b$ を含む停止時刻で定義される場合、生存確率は次を満たす：
   $$P(\tau > t) = \exp\left( -\int_0^t f(s) ds \right)$$
   この結果は、既存のSKIFPPの一意性の結果（[8]および[9]より引用）と、PDEの解を確率的な停止時刻へと写像するスケーリング引数に基づいている。

3. 一意性と正則性:
   本論文は、連続な $\omega$ に対してFBPの弱解の一意性を証明し（命題9）、$\omega(x) = \mathbb{1}_{x \le 0}$ のケースについては先行文献に依拠している。また、極限測度が密度を持ち、無限遠に質量を持たないことを確立している（命題5）。

意義と主張
本論文は、一般化された逆第一通過問題に対する構成的な確率的近似を提供することを主張している。SKIFPPを特定の分岐選択粒子系の流体力学的極限と結びつけることで、本研究はこれらの問題を解くための数値的および理論的な枠組みを提示している。

著者は、古典的なIFPP（指数分布を用いたもの）を標準的な $N$-BBM に結びつけた Klump [13] や De Masiら [4] の先行研究を一般化していると述べている。本論文は、一般的な単調重み関数 $\omega$ と時間依存の分岐レート $f(t)$ を用いた「ソフトキル」の設定へとその接続を拡張しており、これにより、デフォルトが即時的ではなく累積的に発生するという、より幅広い信用リスクモデルをカバーしている。その意義は、与えられた目標分布に対して解 $b(t)$ を近似するために必要な、具体的な粒子系のダイナミクス（分岐レート $f(t)$ および選択重み $\omega$）を特定した点にある。