We propose a novel time discretization for the log-normal SABR model which is a popular stochastic volatility model that is widely used in financial practice. Our time discretization is a variant of the Euler–Maruyama scheme. We study its asymptotic properties in the limit of a large number of time steps under a certain asymptotic regime which includes the case of finite maturity, small vol-of-vol and large initial volatility with fixed product of vol-of-vol and initial volatility. We derive an almost sure limit and a large deviations result for the log-asset price in the limit of a large number of time steps. We derive an exact representation of the implied volatility surface for arbitrary maturity and strike in this regime. Using this representation, we obtain analytical expansions of the implied volatility for small maturity and extreme strikes, which reproduce at leading order known asymptotic results for the continuous time model.

In quantitative finance, it is often necessary to analyze the distribution of the sum of specific functions of observed values at discrete points of an underlying process. Examples include the probability density function, the hedging error, the Asian option, and statistical hypothesis testing. We propose a method to calculate such a distribution, utilizing a recursive method, and examine it using various examples. The results of the numerical experiment show that our proposed method has high accuracy.

In this paper, we establish a fluid limit for a two-sided Markov order book model. The main result states that in a certain asymptotic regime, a pair of measure-valued processes representing the “sell-side shape” and “buy-side shape” of an order book converges to a pair of deterministic measure-valued processes in a certain sense. We also test the fluid approximation on data. The empirical results suggest that the approximation is reasonably good for liquidly traded stocks in certain time periods.

In this paper, we present a pricing model for vulnerable options in discrete time. A Generalized Autoregressive Conditional Heteroscedasticity process is used to describe the variance of the underlying asset, which is correlated with the returns of the asset. As for counterparty default risk, we study it in a reduced form model and the proposed model allows for the correlation between the intensity of default and the variance of the underlying asset. In this framework, we derive a closed-form solution for vulnerable options and investigate quantitative impacts of counterparty default risk on option prices.

Variance reduction has always been a central issue in Monte Carlo experiments. Population Monte Carlo can be used to this effect, in that a mixture of importance functions, called a D-kernel, can be iterativelyoptimized to achieve the minimum asymptotic variance for a function of interest among all possible mixtures.The implementation of this iterative scheme is illustrated for the computation of the price of a Europeanoption in the Cox-Ingersoll-Ross model. A Central Limit theorem as well as moderate deviations are established for the D-kernel Population Monte Carlo methodology.


A generalized correlated random walk is a process of partial sums such that (X, Y) forms a Markov chain. For a sequence (Xn) of such processes in which each takes only two values, we prove weak convergence to a diffusion process whose generator is explicitly described in terms of the limiting behaviour of the transition probabilities for the Yn. Applications include asymptotics of option replication under transaction costs and approximation of a given diffusion by regular recombining binomial trees.

Let us consider n stocks with dependent price processes each following a geometric Brownian motion. We want to investigate the American perpetual put on an index of those stocks. We will provide inner and outer boundaries for its early exercise region by using a decomposition technique for optimal stopping.