An approximation of small-time probability density functions in a general jump diffusion model

We characterize the behavior of the Rough Heston model introduced by Jaisson and Rosenbaum (2016, Ann. Appl. Probab., 26, 2860–2882) in the small‐time, large‐time, and (i.e., ) limits. We show that the short‐maturity smile scales in qualitatively the same way as a general rough stochastic volatility model , and the rate function is equal to the Fenchel–Legendre transform of a simple transformation of the solution to the same Volterra integral equation (VIE) that appears in El Euch and Rosenbaum (2019, Math. Financ., 29, 3–38), but with the drift and mean reversion terms removed. The solution to this VIE satisfies a space–time scaling property which means we only need to solve this equation for the moment values of and so the rate function can be efficiently computed using an Adams scheme or a power series, and we compute a power series in the log‐moneyness variable for the asymptotic implied volatility which yields tractable expressions for the implied vol skew and convexity which is useful for calibration purposes. We later derive a formal saddle point approximation for call options in the Forde and Zhang (2017) large deviations regime which goes to higher order than previous works for rough models. Our higher‐order expansion captures the effect of both drift terms, and at leading order is of qualitatively the same form as the higher‐order expansion for a general model which appears in Friz et al. (2018, Math. Financ., 28, 962–988). The limiting asymptotic smile in the large‐maturity regime is obtained via a stability analysis of the fixed points of the VIE, and is the same as for the standard Heston model in Forde and Jacquier (2011, Finance Stoch., 15, 755–780). Finally, using Lévy's convergence theorem, we show that the log stock price tends weakly to a nonsymmetric random variable as (i.e., ) whose moment generating function (MGF) is also the solution to the Rough Heston VIE with , and we show that tends weakly to a nonsymmetric random variable as , which leads to a nonflat nonsymmetric asymptotic smile in the Edgeworth regime, where the log‐moneyness as , and we compute this asymptotic smile numerically. We also show that the third moment of the log stock price tends to a finite constant as (in contrast to the Rough Bergomi model discussed in Forde et al. (2020, Preprint) where the skew flattens or blows up) and the process converges on pathspace to a random tempered distribution which has the same law as the hyper‐rough Heston model, discussed in Jusselin and Rosenbaum (2020, Math. Finance, 30, 1309–1336) and Abi Jaber (2019, Ann. Appl. Probab., 29, 3155–3200).

We consider variance-optimal hedging in general continuous-time affine stochastic volatility models. The optimal hedge and the associated hedging error are determined semiexplicitly in the case that the stock price follows a martingale. The integral representation of the solution opens the door to efficient numerical computation. The setup includes models with jumps in the stock price and in the activity process. It also allows for correlation between volatility and stock price movements. Concrete parametric models will be illustrated in a forthcoming paper.

Two issues in using mixtures of polynomials for inference in hybrid Bayesian networks

A Hybrid Asymptotic Expansion Scheme: An Application to Long-Term Currency Options

A semi-analytic pricing formula for lookback options under a general stochastic volatility model

We develop an efficient Monte Carlo method for the valuation of financial contracts on discretely realized variance. We work with a general stochastic volatility model that makes realized variance dependent on the full path of the asset price. The variance contract price is a high-dimensional integral over the fundamental sources of randomness. We identify a two-dimensional manifold that drives most of the uncertainty in realized variance, and we compute the contract price by combining precise integration over this manifold, implemented as fine stratification or deterministic sampling with quasirandom numbers, with conditional Monte Carlo on the remaining dimensions. For a subclass of models and a class of nonlinear payoffs, we derive approximate theoretical results that quantify the variance reduction achieved by our method. Numerical tests for the discretized versions of the widely used Hull–White and Heston models show that the algorithm performs significantly better than a standard Monte Carlo, even for fixed computational budgets.

In this paper we consider the range of prices consistent with no arbitrage for European options in a general stochastic volatility model. We give conditions under which the infimum and the supremum of the possible option prices are equal to the intrinsic value of the option and to the current price of the stock, respectively, and show that these conditions are satisfied in most of the stochastic volatility models from the financial literature. We also discuss properties of Black–Scholes hedging strategies in stochastic volatility models where the volatility is bounded.

The Smile in Stochastic Volatility Models

A new method for the approximation of multivariate scalar probability density functions (PDFs) in turbulent reacting flow by means of a joint presumed discrete distribution (jPDD) is presented. The jPDDs can be generated with specified mean values and variances as well as covariances. Correlations between variables – e.g. fluctuating mixture fractions and/or reaction progress – can thereby be taken into account. In this way the new approach overcomes an important limitation of ordinary presumed PDF methods, where statistical independence between the variables is often assumed. Different methods are presented to generate discrete distributions, based either on biased random number generators or on mixing models familiar from PDF transport models. The new approach is extensively validated on a turbulent flow configuration with simultaneous mixing and reaction. Large eddy simulation data as well as results from a transported PDF model are used for the validation of the jPDD approach. The comparison shows that in particular distributions generated with mixing models are able to predict mean reaction rates accurately. For the configuration considered, the neglect of correlations results in significant underestimation of reaction rates. Moreover it is found that higher statistical moments (e.g. the skewness) can influence reaction rates. The consequences for the generation of jPDDs are discussed. In summary, the new jPDD model has the potential to be significantly more accurate than established presumed PDF methods, because correlations between fluctuating variables can be taken into account. At the same time, the new approach is nearly as efficient as standard presumed PDF formulations, since mean rates are computed in a pre-processing step and stored in look-up tables as a function of the first and second moments of the relevant variables.

Complete Analytical Solution of the American Style Option Pricing with Constant and Stochastic Volatilities: A Probability Density Function Approach

A probability density function discretization and approximation method for the dynamic load identification of stochastic structures

In this paper, we propose a fractional stochastic volatility jump-diffusion model which extends the Bates(1996) model. Where we model the volatility as a fractional process. Extensive empirical studies show that the distributions of the logarithmic returns of financial asset usually exhibit properties of self-similarity and long-range dependence and since the fractional Brownian motion has these two important properties, it has the ability to capture the behavior of underlying asset price. Further incorporating jumps into the stochastic volatility framework gives further freedom to financial mathematicians to fit both the short and long end of the implied volatility surface. We propose a stochastic model which contains both fractional and jump process. Then we price options using Monte Carlo simulations along with a variance reduction technique(antithetic variates). We use market data from the S&P 500 index and we compare our results with the Heston and Bates model using error measures. The results show our model greatly outperforms previous models in terms of estimation accuracy.

Pricing of interest rate derivatives, such as CMS spread or mid-curve options, depends on the modelling of the underlying single rates. For flexibility and realism, these rates are often described in the framework of stochastic volatility models. In this paper we allow rates to be modelled within a class of general stochastic volatility models, which includes the SABR, ZABR, free SABR and Heston models. We provide a versatile technique called Effective Markovian Projection, which allows a tractable model to be found that mimics the distribution of the more complex models used to price multi-rate derivatives. Three different numerical approaches are outlined and applied to relevant examples from practice. Finally, a new method that involves moment-matching of Johnson distributions is applied to facilitate closed-form pricing formulas.

In this paper, we present some results on Geometric Asian option valuation for affine stochastic volatility models with jumps. We shall provide a general framework into which several different valuation problems based on some average process can be cast, and we shall obtain closed form solutions for some relevant affine model classes.

Risk premium and fair option prices under stochastic volatility: the HARA solution