Numerical analysis for a unified 2 factor model of structural and reduced form types for corporate bonds with fixed discrete coupon

Technische Samenvatting: Numerieke Analyse voor een Unificerend 2-Factor Model van Bedrijfsobligaties met Vast Discrete Coupons

Probleemstelling
Het artikel behandelt de prijsbepaling van verzuimgevoelige bedrijfsobligaties die vaste discrete coupons uitkeren. Hoewel er aanzienlijk theoretisch werk bestaat voor nulcouponobligaties (Merton, 1977) en structurele modellen voor couponobligaties (Geske, 1977; Agliardi, 2011), is er een gat in de analytische oplossingen voor een unificerend model dat structurele en reduced-form benaderingen combineert, specifiek voor obligaties met vaste discrete coupons.

De auteurs merken op dat hoewel een unificerend één-factor model voor vaste discrete coupons een analytische prijsformule heeft, en een unificerend twee-factor model voor gedisconteerde discrete coupons is opgelost, het unificerende twee-factor model voor vaste discrete coupons schijnbaar geen analytische oplossing heeft. Dit komt door de complexiteit van de onderliggende partiële differentiaalvergelijkingen (PDE's), die parabolische vergelijkingen zijn met drie variabelen (bedrijfswaarde, korte rente en tijd) en gemengde afgeleiden. Bovendien moet het model rekening houden met twee soorten verzuim:

1. Verwacht Verzuim: Treedt alleen op op discrete coupondata wanneer het eigen vermogen van het bedrijf onvoldoende is om schulden en coupons te dekken.
2. Onverwacht Verzuim: Kan op elk moment voorkomen, gestuurd door een constante verzuimintensiteit (hazard rate).

Methodologie
Om dit probleem op te lossen, stellen de auteurs een unificerend twee-factor model voor onder de risiconeutrale martingale-maatstaf, waarbij het volgende wordt geïntegreerd:

• Vasicek-model: Voor de stochastische korte rente ($r$).
• Geometrische Brownse Beweging: Voor de bedrijfswaarde ($V$), met een constante dividendrendement.
• Correlatie: Een correlatiecoëfficiënt ($\rho$) tussen de Wiener-processen die de korte rente en de bedrijfswaarde aansturen.

De prijsbepaling wordt geformuleerd als een reeks beginwaardeproblemen voor PDE's (Vergelijking 2.4). De oplossing bestaat uit een sequentie van PDE's gedefinieerd over sub-intervallen tussen de coupondata, waarbij de randvoorwaarden wisselen tussen de obligatieprijs plus de volgende coupon en de herstelwaarde, afhankelijk van de vraag of de bedrijfswaarde de schuld dekt.

Omdat een analytische formule voor het algemene geval ($i < N-1$) als onhandelbaar wordt beschouwd, maken de auteurs gebruik van een expliciet eindige differentie-schema voor numerieke analyse. Belangrijke methodologische stappen zijn:

1. Variabele Transformatie: Het transformeren van de bedrijfswaarde-variabele $V$ naar $x = \ln V$ om het domein $V > 0$ te behandelen als $x \in \mathbb{R}$.
2. Discretisatie: Het gebruik van voorwaartse differenties voor tijdsafgeleiden en centrale differenties voor ruimtelijke afgeleiden. Voor gemengde afgeleiden ($\partial^2 B / \partial x \partial r$) verkennen de auteurs centrale, voorwaartse en gemengde voorwaarts-achterwaartse differentie-schema's, afhankelijk van de teken van de correlatie $\rho$.
3. Stabiliteitsanalyse: Een rigoureuze afleiding van de stabiliteitsvoorwaarden voor het expliciete schema. De auteurs tonen aan dat het standaard expliciete schema over het algemeen instabiel is vanwege de gemengde afgeleiden. Zij leiden specifieke beperkingen af voor de tijdstap ($\Delta t$) en de ruimtelijke stappen ($\Delta x, \Delta r$) relatief aan de diffusiecoëfficiënten en de correlatie $\rho$ om stabiliteit te garanderen.

Belangrijkste Bijdragen

1. Unificerend Model Formule: Het artikel formaliseert een unificerend twee-factor model dat een structureel verzuim (barrière-overschrijding op discrete tijdstippen) en reduced-form verzuim (gestuurd door intensiteit op elk moment) integreert voor obligaties met vaste discrete coupons.
2. Stabiliteitsvoorwaarden: De primaire theoretische bijdrage is de afleiding van de stabiliteitsvoorwaarden vereist om het expliciete eindige differentie-schema toepasbaar te maken voor dit specifieke 3-variabele parabolische PDE-systeem met gemengde afgeleiden.
   • Voor $\rho = 0$ vereist stabiliteit $\mu_x + \mu_r < 1$ en specifieke grenzen op de drift-termen.
   • Voor $0 < \rho < 1/2$ laten de auteurs zien dat stabiliteit vereist dat $0 < \rho < 1/2$ en specifieke ratio's tussen de ruimtelijke stapgroottes ($\sqrt{\mu_r/\mu_x}$).
   • Voor $-1 < \rho < 0$ worden vergelijkbare beperkingen afgeleid met behulp van een gemengde voorwaarts-achterwaartse differentie voor de gemengde afgeleide.
3. Numerieke Implementatie: De auteurs bieden een computationeel kader om obligatieprijzen, kredietspreads en duraties te berekenen waar analytische formules niet beschikbaar zijn.

Resultaten
Gebruikmakend van het afgeleide expliciete schema met $\rho = 0$, presenteren de auteurs numerieke resultaten gebaseerd op een specifieke set parameters (bijv. $N=2$ coupons, Vasicek-parameters, specifieke verzuimintensiteiten $\lambda_0, \lambda_1$). De analyse levert het volgende op:

• Obligatieprijzen:
  • Prijzen dalen naarmate de korte rente ($r$) stijgt.
  • Prijzen stijgen naarmate de bedrijfswaarde ($V$) stijgt, waarbij convexiteit ten opzichte van $V$ wordt vertoond.
  • Prijzen vertonen sprongen (discontinuïteiten) op coupondata, wat de betaling van coupons en de reset van de verzuimrisicovoorwaarden reflecteert.
• Kredietspread (CS):
  • CS daalt naarmate de bedrijfswaarde stijgt.
  • CS stijgt snel naarmate de looptijd nadert als de bedrijfswaarde laag is.
  • CS daalt naarmate de korte rente stijgt (door het verdisconterings-effect op het risicovrije component).
  • CS stijgt met een hogere verzuimintensiteit ($\lambda$).
  • CS vertoont sprongen op coupondata; de omvang van de sprong neemt toe met de grootte van de coupon.
• Duration:
  • Duration neemt toe met de bedrijfswaarde en de korte rente.
  • Duration vertoont sprongen na couponbetalingen, consistent met de definitie van duration als een gewogen gemiddelde tijd tot kasstromen.

Betekenis en Claims
Het artikel claimt bescheiden een noodzakelijk numeriek instrument te bieden voor een model dat Agliardi's structurele model generaliseert naar een unificerend kader voor obligaties met vaste discrete coupons. De auteurs stellen dat het verkrijgen van een analytische formule voor dit specifieke unificerende twee-factor model "moeilijk" en waarschijnlijk onmogelijk is in het algemene geval. Daarom ligt de betekenis van het werk in:

1. Het vaststellen van de stabiliteitsvoorwaarden die nodig zijn om het expliciete eindige differentie-schema levensvatbaar te maken voor dit complexe PDE-systeem.
2. Het aantonen dat de numerieke resultaten compatibel zijn met reële financiële situaties, zoals blijkt uit het intuïtieve gedrag van obligatieprijzen, kredietspreads en duraties onder variërende marktomstandigheden (rente, bedrijfswaarde, verzuimintensiteit).

Het werk dient als een brug tussen de theoretische unificatie van structurele/reduced-form modellen en de praktische prijsbepaling voor realistische obligatiestructuren (vaste discrete coupons) waar gesloten oplossingen niet bestaan.