Branching-selection particle systems and inverse first passage problems

Technische Zusammenfassung: Verzweigungs-Selektions-Partikelsysteme und inverse First-Passage-Probleme

Problemstellung
Diese Arbeit befasst sich mit dem verallgemeinerten inversen First-Passage-Problem (IFPP). Das klassische IFPP, das aus einer Frage von Shiryaev hervorging, fragt, ob es für eine Standard-Brownsche Bewegung $(W_t)_{t \ge 0}$ und eine Zielwahrscheinlichkeitsverteilung $p$ auf $[0, \infty]$ eine Grenze $b(t)$ gibt, sodass die First-Passage-Zeit $\tau := \inf\{t : W_t \le b(t)\}$ der Verteilung $p$ folgt.

Die vorliegende Arbeit konzentriert sich auf eine Verallgemeinerung, die als Soft-Killed Inverse First Passage Problem (SKIFPP) bezeichnet wird. Die Haltezeit ist definiert als:
$$\tau := \inf \left\{ t : \Lambda \int_0^t \omega(W_s - b(s)) ds \ge U \right\}$$
wobei $U \sim \text{Exp}(1)$, $\Lambda \in (0, \infty)$ und $\omega: \mathbb{R} \to [0, 1]$ eine monoton fallende Funktion ist. Im Gegensatz zum klassischen Fall, in dem ein Ausfall unmittelbar beim Überschreiten einer Grenze eintritt, modelliert diese Formulierung den Ausfall als kumulative Funktion der Historie des Pfades relativ zur Grenze. Das Ziel besteht darin, eine Grenze $b(t)$ zu konstruieren, sodass $\tau \sim p$, gegeben eine spezifische Anfangsverteilung für $W_0$ und eine Ziel-Überlebenswahrscheinlichkeit $p([t, \infty))$.

Methodik
Der Autor konstruiert ein Verzweigungs-Selektions-Partikelsystem (speziell eine Variante der $N$-Branching-Brownian-Motion oder $N$-BBM), um die Lösung des SKIFPP zu approximieren. Die Methodik verläuft wie folgt:

1. Konstruktion des Partikelsystems:

   • Betrachte $N$ Teilchen, die sich als unabhängige Brownsche Bewegungen bewegen.
   • Teilchen verzweigen mit einer zeitabhängigen Rate $f(t)$.
   • Gleichzeitig mit der Verzweigung wird ein Teilchen gelöscht. Die Löschung erfolgt „soft-killed“ basierend auf einer Selektionsgewichtungsfunktion $\omega$. Konkret wird bei Verzweigungszeiten ein Teilchen an der Position $x$ mit einer Wahrscheinlichkeit gelöscht, die proportional zu $\omega(x - b_N(t))$ ist, wobei $b_N(t)$ eine bewegliche Grenze ist, die durch die empirische Verteilung der Teilchen bestimmt wird.
   • Die Grenze $b_N(t)$ ist implizit definiert, um sicherzustellen, dass die gesamte Verzweigungsrate der vorgegebenen Rate $Nf(t)$ entspricht:
     $$b_N(t) := \inf \left\{ x : \Lambda \langle \mu_N^t, \omega(\cdot - x) \rangle \ge f(t) \right\}$$
     wobei $\mu_N^t$ das empirische Maß der $N$ Teilchen ist.

2. Analyse des hydrodynamischen Limits:

   • Der Autor analysiert das Limit des empirischen Maßes $\mu_N^t$ für $N \to \infty$.
   • Unter Verwendung von Martingalmethoden (speziell der Itō-Formel und der Burkholder-Davis-Gundy-Ungleichung) leitet der Autor die Evolutionsgleichung für die Testfunktionen des empirischen Maßes her.
   • Die Tightness der Folge der Gesetze wird im Raum der maßwertigen Prozesse mit der Skorokhod-Topologie etabliert.
   • Es wird gezeigt, dass das Limit ein Free-Boundary-Problem (FBP) erfüllt, das durch eine partielle Differentialgleichung (PDE) gesteuert wird.

Wichtigste Ergebnisse

1. Konvergenz gegen das Free-Boundary-Problem (Theorem 1):
   Das empirische Maß $\mu_N^t$ des konstruierten $(\omega, f, \Lambda, N)$-BBM konvergiert schwach gegen ein deterministisches Maß $\mu_t$ mit Dichte $u(x, t)$. Diese Dichte ist die eindeutige schwache Lösung des folgenden Free-Boundary-Problems:
   $$\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}$$
   wobei die freie Grenze $b(t)$ durch die Integralbedingung $\int_{\mathbb{R}} \omega(x - b(t)) u(x, t) dx = f(t)/\Lambda$ bestimmt wird.

2. Verbindung zum SKIFPP (Theorem 2):
   Die Arbeit zeigt, dass, falls $(u, b)$ das obige Free-Boundary-Problem löst, die Grenze $b(t)$ die eindeutige Lösung des SKIFPP ist. Speziell gilt, wenn $\tau$ durch die Haltezeit unter Einbeziehung von $\omega$ und $b$ definiert ist, dass die Überlebenswahrscheinlichkeit die folgende Form hat:
   $$P(\tau > t) = \exp\left( -\int_0^t f(s) ds \right)$$
   Dieses Ergebnis stützt sich auf bestehende Unizitätsresultate für das SKIFPP (zitiert aus [8] und [9]) und ein Skalierungsargument, um die PDE-Lösung auf die probabilistische Haltezeit abzubilden.

3. Eindeutigkeit und Regularität:
   Es wird bewiesen, dass die Arbeit die Eindeutigkeit der schwachen Lösung des FBP für kontinuierliche $\omega$ (Proposition 9) besitzt und stützt sich auf die bestehende Literatur für den Fall $\omega(x) = \mathbb{1}_{x \le 0}$. Zudem wird etabliert, dass das Grenzmaß eine Dichte besitzt und fast sicher keine Masse im Unendlichen aufweist (Proposition 5).

Bedeutung und Ansprüche
Die Arbeit beansprucht, eine konstruktive probabilistische Approximation für das verallgemeinerte inverse First-Passage-Problem bereitzustellen. Durch die Verknüpfung des SKIFPP mit dem hydrodynamischen Limit eines spezifischen Verzweigungs-Selektions-Partikelsystems bietet die Arbeit einen numerischen und theoretischen Rahmen zur Lösung dieser Probleme.

Der Autor stellt fest, dass dies die bisherige Arbeit von Klump [13] und De Masi et al. [4] verallgemeinert, welche das klassische IFPP (mit exponentiellen Verteilungen) mit der Standard-$N$-BBM verknüpften. Diese Arbeit erweitert diese Verbindung auf das „Soft-Killed“-Setting mit allgemeinen monoton fallenden Gewichtungsfunktionen $\omega$ und zeitabhängigen Verzweigungsraten $f(t)$, wodurch eine breitere Klasse von Kreditausfallmodellen abgedeckt wird, bei denen der Ausfall nicht instantan erfolgt, sondern über die Zeit akkumuliert. Die Bedeutung liegt in der Identifizierung der spezifischen Dynamik des Partikelsystems (Verzweigungsrate $f(t)$ und Selektionsgewicht $\omega$), die erforderlich ist, um die Lösung $b(t)$ für eine gegebene Zielverteilung zu approximieren.