We derive a stochastic version of the Magnus expansion for the solution of linear systems of Ito stochastic differential equations (SDEs). The goal of this paper is twofold. First, we prove existence and a representation formula for the logarithm associated to the solution of the matrix-valued SDEs. Second,we propose a new method for the numerical solution of stochastic partial differential equations (SPDEs) based on spatial discretization and application of the stochastic Magnus expansion. A notable feature oft he method is that it is fully parallelizable. We also present numerical tests in order to asses the accuracy of the numerical schemes.
ESR6 Kevin Kamm
We derive the stochastic version of the Magnus expansion for linear systems of stochastic differential equations (SDEs). The main novelty with respect to the related literature is that we consider SDEs in the Itô sense, with progressively measurable coefficients, for which an explicit Itô-Stratonovich conversion is not available. We prove convergence of the Magnus expansion up to a stopping time τ and provide a novel asymptotic estimate of the cumulative distribution function of τ . As an application, we propose a new method for the numerical solution of stochastic partial differential equations (SPDEs) based on spatial discretization and application of the stochastic Magnus expansion. A notable feature of the method is that it is fully parallelizable.We also present numerical tests in order to asses the accuracy of the numerical schemes.
In this paper, we propose a new model to address the problem of negative interest rates that preserves the analytical tractability of the original Cox–Ingersoll–Ross (CIR) model without introducing a shift to the market interest rates, because it is defined as the difference of two independent CIR processes. The strength of our model lies within the fact that it is very simple and can be calibrated to the market zero yield curve using an analytical formula.We run several numerical experiments at two different dates, once with a partially sub-zero interest rate and once with a fully negative interest rate. In both cases, we obtain good results in the sense that the model reproduces the market term structures very well. We then simulate the model using the Euler–Maruyama scheme and examine the mean, variance and distribution of the model. The latter agrees with the skewness and fat tail seen in the original CIR model. In addition, we compare the model’s zero coupon prices with market prices at different future points in time. Finally, we test the market consistency of the model by evaluating swaptions with different tenors and maturities.