Short-rate models with stochastic discontinuities: a PDE approach
Alessandro Calvia, Marzia De Donno, Chiara Guardasoni, Simona Sanfelici
TL;DR
The paper develops a PDE-based framework for pricing interest rate derivatives under short-rate dynamics that feature stochastic discontinuities at fixed times and a discontinuous numéraire at roll-over dates. By exploiting a Markovian representation and no-arbitrage, it derives a system of backward PDEs with jump boundary conditions and proves existence, uniqueness, and a Feynman-Kač representation for the solution, including exponential-growth controls. In the affine setting, it obtains closed-form ZCB prices and, for extended Vasicek with Gaussian jumps, closed-form call prices on ZCB, while for general cases it proposes two robust numerical schemes: a semi-analytic Green’s-function method and a finite-difference approach. Numerical results validate the methods against known solutions and demonstrate the framework’s flexibility to handle jumps in both the short-rate and the numéraire, with practical implications for pricing and risk management of RFR-based derivatives. The work thus provides a tractable, adaptable tool for post-LIBOR term structure modeling and derivative pricing in environments with regulatory-driven discontinuities.
Abstract
With the reform of interest rate benchmarks, interbank offered rates (IBORs) like LIBOR have been replaced by risk-free rates (RFRs), such as the Secured Overnight Financing Rate (SOFR) in the U.S. and the Euro Short-Term Rate (\euro STR) in Europe. These rates exhibit characteristics like jumps and spikes that correspond to specific market events, driven by regulatory and liquidity constraints. To capture these characteristics, this paper considers a general short-rate model that incorporates discontinuities at fixed times with random sizes. Within this framework, we introduce a PDE-based approach for pricing interest rate derivatives and establish, under suitable assumptions, a Feynman-Kač representation for the solution. For affine models, we derive (quasi) closed-form solutions, while for the general case, we develop numerical methods to solve the resulting PDEs.
