We construct a new stochastic interest rate model with two stochastic factors, by introducing a stochastic long-run equilibrium level into the Vasicek interest rate model which follows another Ornstein–Uhlenbeck process. With the interest rate under the Black–Scholes model being assumed to follow the newly proposed model, a closed-form representation of European option prices is successfully presented, when the analytical characteristic function of the underlying log-price under a forward measure is derived. To assess the model performance, a preliminary empirical study is conducted using S&P 500 index and its options, with the Vasicek model and an alternative two-factor Vasicek model taken as benchmarks.

Zhu and He [(2018). A new closed-form formula for pricing European options under a skew Brownian motion. The European Journal of Finance 24(12): 1063–1074] provided an innovative closed-form solution by replacing the standard Brownian motion in the Black–Scholes framework using a particular skew Brownian motion. Their formula involves numerically integrating the product of the Guassian density and corresponding distribution function. Being different from their pricing formula, we derive a much simpler formula that only involves the Gaussian distribution function and Owen's $T$ function.

We combine the rough Heston model and the CIR (Cox–Ingersoll–Ross) interest rate together to form a rough Heston-CIR model, so that both the rough behaviour of the volatility and the stochastic nature of the interest rate can be captured. Despite the convoluted structure and non-Markovian property of this model, it still admits a semi-analytical pricing formula for European options, the implementation of which involves solving a fractional Riccati equation. The rough Heston-CIR model is more general, taking both the rough Heston model and the Heston-CIR model as special cases. The influence of rough volatility and stochastic interest rate is shown to be significant through numerical experiments.

We study the value of European security derivatives in the Black–Scholes model when the underlying asset $\xi $ is approximated by random walks ${\xi }^{(n)} $. We obtain an explicit error formula, up to a term of order $ \mathcal{O} ({n}^{- 3/ 2} )$, which is valid for general approximating schemes and general payoff functions. We show how this error formula can be used to find random walks ${\xi }^{(n)} $ for which option values converge at a speed of $ \mathcal{O} ({n}^{- 3/ 2} )$.

We propose a variational analysis for a Black and Scholes equation with stochastic volatility. This equation gives the price of a European option as a function of the time, of the price of the underlying asset and of the volatility when the volatility is a function of a mean reverting Orstein-Uhlenbeck process, possibly correlated with the underlying asset. The variational analysis involves weighted Sobolev spaces. It enables to prove qualitative properties of the solution, namely a maximum principle and additional regularity properties. Finally, we make numerical simulations of the solution, by finite element and finite difference methods.