Search

We contribute to the recent studies of the so-called Bass martingale. Backhoff-Veraguas et al. (2020) showed it is the solution to the martingale Benamou-Brenier (mBB) problem, i.e., among all martingales with prescribed initial and terminal distributions it is the one closest to the Brownian motion. We link it with semimartingale optimal transport and deduce an alternative way to derive the dual formulation recently obtained in Backhoff-Veraguas et al. (2023). We then consider computational methods to compute the Bass martingale. The dual formulation of the transport problem leads to an iterative scheme that mirrors to the celebrated Sinkhorn algorithm for entropic optimal transport. We call it the measure preserving martingale Sinkhorn (MPMS) algorithm. We prove that in any dimension, each step of the algorithm improves the value of the dual problem, which implies its convergence. Our MPMS algorithm is equivalent to the fixed-point method of Conze and Henry-Labordere (2021), studied in Acciaio et al. (2023), and performs very well on a range of examples, including real market data.

We develop and implement a non-parametric method for joint exact calibration of a local volatility model and a correlated stochastic short rate model using semimartingale optimal transport. The method relies on the duality results established in Joseph, Loeper, and Obloj, 2023 and jointly calibrates the whole equity-rate dynamics. It uses an iterative approach which starts with a parametric model and tries to stay close to it, until a perfect calibration is obtained. We demonstrate the performance of our approach on market data using European SPX options and European cap interest rate options. Finally, we compare the joint calibration approach with the sequential calibration, in which the short rate model is calibrated first and frozen.

We develop a non-parametric, optimal transport driven, calibration methodology for local volatility models with stochastic interest rate. The method finds a fully calibrated model which is the closest to a given reference model. We establish a general duality result which allows to solve the problem via optimising over solutions to a non-linear HJB equation. We then apply the method to a sequential calibration setup: we assume that an interest rate model is given and is calibrated to the observed term structure in the market. We then seek to calibrate a stock price local volatility model with volatility coefficient depending on time, the underlying and the short rate process, and driven by a Brownian motion which can be correlated with the randomness driving the rates process. The local volatility model is calibrated to a finite number of European options prices via a convex optimisation problem derived from the PDE formulation of semimartingale optimal transport. Our methodology is analogous to Guo, Loeper, and Wang, 2022 and Guo, Loeper, Obloj, et al., 2022a but features a novel element of solving for discounted densities, or sub-probability measures. We present numerical experiments and test the effectiveness of the proposed methodology.