Stochastic Volatility Effects on Correlated Log-Normal Random Variables

The transition density function plays an important role in understanding and explaining the dynamics of the stochastic process. In this paper, we incorporate an ergodic process displaying fast moving fluctuation into constant volatility models to express volatility clustering over time. We obtain an analytic approximation of the transition density function under our stochastic process model. Using perturbation theory based on Lie–Trotter operator splitting method, we compute the leading-order term and the first-order correction term and then present the left and right skew scenarios through numerical study.


Introduction
The multivariate log-normal distribution is a widely used stochastic model in social sciences.What is the probability of the sum or difference of log-normal random variables?The solution to this question has wide applications in many fields such as finance [1,2], actuarial science [3,4], and physics [5].Especially, in physics, examples include wave transmittance in random media, dissipation rate of turbulence energy, and temporal fluctuations of some nonlinear systems.However, to the best of the author of this paper's knowledge, almost nothing is known about the question yet.So, finding analytical approximations to it is as important as an alternative approach (see [6] for details).The main advantage of the analytical approximation approaches compared to other numerical methods is that in general the first ones are much faster and precise at least under certain model parameter regime.In addition, analytic approximation formulas retain qualitative model information and preserve an explicit dependence of the results on the underlying parameters.Many approximate methods can be categorized by their scope into three classes: generally correlated [7][8][9][10], independent [11][12][13][14], and independent and identically distributed [15,16].The probability density function of the correlated log-normal especially is carried out to propose various approximations of the sum distribution.These approximations can be divided into two categories.The first method is approximated by another log-normal random variable.For example, the lognormal parameters are obtained by moment matching [7], the mean and the variance are computed recursively [17], and the shifted log-normal distribution [18].A second method of these methods is introduced to construct efficient numerical integrations techniques such as a steepest-descent integration [19] and a Smolyak's algorithm [20].
In option pricing theory, it is well-known that constant volatility for the stock price in the Black-Scholes model [21] can hardly capture the accumulated empirical evidence in financial markets since the parameters contained in the model actually change over time.A main drawback in the assumption of this model lies in the flat implied volatility surface, which is contradictory to empirical results that the implied volatilities of the equity options exhibit the smile or skew curve, in particular, from the 1987 crash onward.Among those several methods of overcoming the above drawback and relaxing the assumption of the Black-Scholes constant volatility model, a variety of stochastic volatility models have been suggested to incorporate the volatility skew.Fast mean reversion especially is one of the notable features of volatility.For example, Fouque et al. [22] estimated that its volatility reverts to the mean with characteristic time about 1.5 days in high-frequency S&P 500 index data.Also, Engle [23] and Bollerslev [24] introduced the family of autoregressive conditional heteroskedasticity and generalized autoregressive conditional heteroskedasticity, respectively, to describe the evolution of the volatility of the stock price in discrete time setting and showed that econometric tests of these models reject the assumption of constant volatility and find evidence of volatility clustering over time.
Based on these observations together with the renowned contribution of stochastic volatility formula to option pricing, we propose the following fast mean-reverting (FMR) model in which asset prices are conditionally log-normal and the volatility process is a positive increasing function of a FMR process: where  1 and  2 are log-normal random variables and  1 ,  2 , and   are standard Brownian motions correlated as follows: with −1 ≤  12 ,  1 ,  2 ≤ 1.Here, the parameter  measures the rate at which the ergodic Markov process  reaches its long-term mean value ,  is a function,  is a constant, and the correlations  1 and  2 control the slope of the skew of each underlying.We assume that  is large so that volatility is sufficiently fast mean-reversion when looked over the time scale of options Fouque et al. [22].We do not specify the concrete forms of  1 ,  2 , and  since they will not play an essential role in the perturbation theory performed in this paper but the functions must satisfy a sufficient growth condition to avoid some kind of bad cases such as nonexistence of the moments of .This generalization is a clear advantage of this model in that the FMR model encompasses the models which reflect stylized facts such as a feedback effect between volatility and volatility of volatility.

Asymptotic Analysis
2.1.Problem Formulation.Now we will denote  1 ±  2 and  10 ±  20 as  ± and  ± 0 , respectively, for simplicity.The probability distribution of the sum or difference of the two correlated log-normal distributions can be obtained by calculating the integral where  is the joint probability distribution of the two lognormal random variables and ( ± −  ± 0 ) is the Dirac delta function.Unfortunately, a closed-form representation for this probability distribution still does not exist.FMR volatility enables us to make good use of perturbation theory.To end this, we use the small positive parameter , which denotes the inverse of the rate of mean-reversion  (which is assumed to be large).We suppose  = ] √ 2/√, where the variance ] 2 of the invariant distribution of  is a constant with respect to .If rewriting (2) in terms of , we have This gives a singular perturbation problem with respect to the small parameter , which furnishes us with analytic tractability.We define the scaled FMR (SFMR) model as the FMR model used with an first-order asymptotic method in this paper.
If we compare the Heston model [25], one of representative stochastic volatility models, with the SFMR model, we can get some strengths.The Heston model has five parameters required to be estimated from market data while the SFMR model has two.In addition, the Heston model has certain restrictions on the volatility functions  1 ,  2 , and  whereas the SFMR model does not have any.

Perturbation Theory Based on Lie-Trotter
Operator Splitting Method.Perturbation theory as developed by Fouque et al. [22] is a methodology utilized to find an approximated solution when the original problem is difficult to solve by separating it into more easily solvable, simple parts.If we apply the Feynman-Kac formula, we find that   ± satisfies the following Kolmogorov backward equation: where Here, L 0 is the infinitesimal generator of the SFMR process .L 1 contains the mixed partial derivatives due to the correlations of the two Brownian motions  1 and   and  2 and   , respectively.L 2 is the operator of a generalized version of the two-dimensional standard Brownian motion at the volatility levels  1 ( 0 ) and  2 ( 0 ) instead of constant volatilities, respectively.Before we solve problem (6), we write a useful lemma about the centering (or solvability) condition on the Poisson equation related to the operator L 0 as follows.)) of .Then, solutions of ( 10) are given by the form Proof.See Fouque et al. [22].
Using the resultant partial differential equation (PDE) (6) and expanding   ± in powers of √, one can approximate   ± to the sum of the leading term  (0) ± and the first correction term √ (1)  ± as follows: Since we focus on the first-order correction term for   ± , we reset (12) with respect to √ (1)  ± and denote it by ũ(1) ± .Using perturbation theory based on Lie-Trotter operator splitting method,  (0) ± and ũ(1) ± must satisfy the following PDEs with boundary conditions, respectively: with (,  1 ,  2 ;  0 → ,  10 ,  20 with ũ(1) ± (,  1 ,  2 ;  0 → ,  10 ,  20 ) = 0 (16) and then one obtains the solutions  (0) ± and ũ(1) ± of the PDEs, respectively, where and the constant parameters  1 ,  2 ,  12 , and  21 are defined as follows: See the Appendix for the relevant technique.As a result, the distributions of the sum and difference of correlated log-normal random variables under SFMR model are shown to follow combining a shifted log-normal distribution with mixed partial derivatives of it.Also, all the original parameters are absorbed in the group parameters  1 ,  2 ,  12 , and  21 and the present level  of the hidden process  driving the fast time-scale volatility needs not be specified in the present approximation.

Numerical Experiment.
In this subsection, we illustrate the effectiveness of our result (12) by showing numerical results.
Several parameter sets are obtained from Lo except for group parameters.Figure 1(a) depicts log-normal distribution (leading-order term or Lo's result) and log-normal distributions under SFMR model with positive or negative group parameters (leading-order term plus first-order correction term), respectively.The solid line denotes the shifted lognormal distribution and the dash lines show combining the shifted log-normal distribution and the first correction term, respectively.Figure 1(b) depicts the errors calculated by subtracting the log-normal distributions with positive or negative group parameters from the log-normal distribution, respectively.We also apply the same result to the difference of the probability density function in Figure 2. Our numerical results show that the left and right skew scenarios are presented through first-order correction term and the major discrepancies appear around the peak of the probability density function.These pictures are sensitive to the choice of the involved parameters and give a lot of flexibility to the shape of the transition densities.

Final Remarks
Stochastic processes are popular in modeling various economics and financial variables.The transition density function especially plays a key role in the analysis of continuoustime diffusion models.In this paper, we obtained an analytic approximation of correlated log-normal random variables under SFMR model.
This paper offers various possible directions for further development.Our result can be applied to pricing and hedging spread options and is to incorporate a slowly varying volatility-driving process into the SFMR model.Also, this result can provide a very useful guide for credit risk management (see [26,27]).We leave these issues as future research topics.