Abstract
As is well known, multi-factor stochastic volatility models are necessary to capture the market accurately in pricing financial derivatives. However, the multi-factor models usually require too many parameters to be calibrated efficiently and they do not lead to an analytic pricing formula. The double Heston model is one of them. The approach of this paper for this difficulty is to rescale the double Heston model to reduce the number of the model parameters and obtain a closed form analytic solution formula for variance swaps explicitly. We show that the rescaled double Heston model is as effective as the original double Heston model in terms of fitting to the VIX market data in a stable condition and yet the computing time is much less than that under the double Heston model. However, in a turbulent situation after the start of the COVID-19 pandemic in 2020, we acknowledge that even the double Heston model fails to capture the market accurately.
Similar content being viewed by others
1 Introduction
Volatility is a degree of uncertainty of an underlying risky asset and the volatility management and trading are considered to be a crucial task for private investors, institutions, and hedge funds. Volatility derivatives are a class of financial derivatives where the payoff is a function of some measure of the volatility. Variance swap is a prominent example of this type of derivatives. It is one of forward contracts for realized variance whose payoff is given by the spread between realized and implied variances. Variance swap users trade volatility levels directly and thus they may have an advantage over vanilla option traders for hedging purposes.
In order to derive a fair strike (delivery) price of variance swap, the Heston volatility model (Heston 1993) has been generally used as the classical one. Its volatility process is driven by the Cox–Ingersoll–Ross (CIR) process (Cox, Ingersoll & Ross 1985) which has a mean reversion property. For example, Swishchuk (2004) studied the pricing of volatility and variance swaps using a probabilistic approach. Broadie and Jain (2008) found the closed form prices of variance swap and volatility swap under the Heston model together with the impact of discrete sampling. Zhu and Lian (2011) derived a closed form solution for a discretely sampled variance swap by manipulating an additional variable introduced by Little and Pant (2001). Zheng and Kwok (2014b) used saddlepoint approximation methods for pricing variance derivatives. Also, Zheng and Kwok (2014a) derived the prices of generalized and exotic variance swaps under the Heston model with jumps. Kim and Kim (2020) obtained affine approximations for the fair strike prices of the generalized variance swaps using the projection techniques of Grzelak and Oosterlee (2011). However, the downside of the Heston model is that it often results in an unsuitable value compared with market prices, particularly when it is applied to the derivatives whose time-to-maturity is short. Also, it is difficult to capture the movement of volatility when the fluctuations of the volatility are relatively large.
So, there have been studies showing that multi-factor stochastic volatility models can explain the market more accurately than the Heston model. Gatheral (2008) proposed the double mean reverting model for consistent modeling of SPX and VIX options. Fouque, Papanicolaou, Sircar and Solna (2003, 2011) developed an asymptotic two-scale factor stochastic volatility model for pricing financial derivatives. Christoffersen, Heston and Jacobs (2009) proposed a two-factor Heston model, named the double Heston model, with fast and slow factors which are all mean reverting and showed that the model can explain the dynamics of market more flexibly than the classical Heston model. Using this model, Gauthier and Possamai (2011) improved the European option pricing formula given by Christoffersen et al. Moreover, Zhang and Feng (2019) and Fallah and Mehrdoust (2019) studied the pricing of American options under the double Heston model. However, it costs too much time for calibration under the double Heston model because the number of the model parameters is twice that of the Heston model.
The contribution of this article is as follows. First, we use a rescaling approach for the original double Heston model to reduce the number of the model parameters. Second, we obtain a closed form analytic solution formula for variance swaps explicitly under the rescaled double Heston model. Third, we show that the pricing error between the rescaled model and the original model is minimal while the computation time with the rescaled model is much reduced compared to the original model. So, the rescaling approach proposed in this article is verified as a possible efficient method for computing the prices of derivatives such as variance swaps.
This paper is organized as follows. We first introduce a reduced version of the double Heston stochastic volatility model and formulate a pricing problem for variance swaps in terms of partial differential equations (PDEs) in Section 2. We use the Fourier transform and Green function methods to solve the PDEs and obtain a closed form solution formula for variance swaps in Section 3. In Section 4, we check the validity of the formula using Monte Carlo simulation technique. Also, we investigate the impacts of the parameters of fast and slow volatility factors on the fair strike prices of variance swaps. We estimate the model parameters from VIX term structures provided by Chicago Board of Options Exchange (CBOE) and compare the Heston model, the double Heston model and the rescaled double Heston model and show the merits of the rescaled model in terms of computation time and accuracy. Finally, Sect 5 concludes.
2 Model and Problem Formulation
Given an underlying asset price process \(X_t\), the (original) double Heston model devised by Christoffersen et al. (2009) is given by
where \(c_i\)’s are mean-reversion rate, \(\theta _i\)’s are long-run mean variances and \(\sigma _y,\sigma _z\) are vol-of-vol, which are all positive constants. It is assumed that there is a zero correlation between \(Y_t\) and \(Z_t\). The market prices of volatility risk have been suppressed here. One volatility process moves as a fast mean reverting factor while another one is slowly mean reverting. The double Heston model may capture the movement of real market well but the drawback is that it has too many parameters for calibration. To calibrate parameters more efficiently, we propose in this paper the following rescaled version of the double Heston model.
where \(\alpha \), \(\beta \), \(\theta _1\) and \(\theta _2\) are all positive. Intuitively, we have in mind that \(\alpha ^{-1}\) plays a role as a fast-scale factor whereas \(\beta \) represents a slow-scale but mathematically, we do not assume that \(\alpha \) and \(\beta \) are small here. The above formulation satisfies the Feller condition, which is a sufficient condition for the volatility to be strictly positive. Note that the number of model parameters is reduced from 8 in (1) to 6 in (2):
As the discretely sampled realized variance of a variance swap over [0, T] is defined by
where \(\left[ 0,T\right] \) is divided into N periods of \([t_{j-1},t_{j}]\), \(t_N=T\) and \(j=1,\ldots ,N\), and
the fair strike, \(K_{\text{var}}\), of the variance swap is given by the strike K such that the expected value of the payoff \(\sigma ^2_{\text{RV}}-K\) is equal to 0 at time \(t=t_0\). Therefore, we have
Here, the tower property of the conditional expectation has been used. Note that the above realized variance is defined in terms of logarithmic return instead of simple return. It has an advantage against the simple return in that the multi-period log return is a sum of the one-period log returns, while the multi-period simple return is a product of the one-period simple returns, leading to computational problems for values close to zero. That is why the log return is preferred in financial industry.
The payoff function in (4) depends on two unknown random variables \(X_{t_{j-1}}\) and \(X_{t_j}\). To deal with this difficult problem, we use the idea proposed by Little and Pant (2001). Namely, considering the time interval \([0,t_j]\), we use a new independent variable \(I_t\) defined by
where the \(\delta (\cdot )\) denotes the generalized Dirac delta function. Note that \(I_t=0\) on the interval \([0,t_{j-1})\) and \(I_t=X_{t_{j-1}}\) on \([t_{j-1},t_{j}]\) and so \(I_t\) experiences a jump in value across time \(t_{j-1}\).
If \(U^j(t,x,y,z;I)\) is the solution of the PDE
with terminal condition
then by the well-known Feynman–Kac theorem (cf. Oksendal 2000), we have
From the property of Dirac delta function, one can get rid of the independent variable \(I_t\) in the PDE (5) except at the time \(t_{j-1}\). However, the terminal condition still depends on \(I_t\) which has a jump at \(t=t_{j-1}\). So, we need a continuity condition at \(t=t_{j-1}\) from the no-arbitrage pricing theory (cf. Wilmott 2013). The condition (a jump condition) is expressed as
Therefore, we divide the time domain \([0,t_{j}]\) into two subintervals \([0,t_{j-1}]\) and \([t_{j-1},t_{j}]\) as the variable \(I_t\) could be regarded as constant on each interval and solve the PDE system by two steps, the first step in \([t_{j-1},t_{j}]\) and the second step in \([0,t_{j-1}]\). The solution obtained in the first step will provide the terminal condition for the PDE problem of the second step through the jump condition. So, on the time interval \([t_{j-1},t_{j}]\), \(U_j(t,x,y,z;I)\) satisfies
with the terminal condition as stated above. On the time interval \([t_0=0,t_{j-1}]\), we use notation \(W^j(t,x,y,z;I)\) instead of \(U^j(t,x,y,z;I)\) for convenience. It satisfies
with terminal condition \(W^{j}(t,x,y,z)|_{t=t_{j-1}}=U^{j}(t_{j-1},x,y,z;x)\).
From the above two step process of solving the PDE problems, we have
and thus putting (4), (6) and (9) together leads to
3 Pricing Variance Swaps
In this section, we analytically solve the PDE system (7) for the inner expectation \(U^j\) and the PDE system (8) for the outer expectation \(W^j\) using the Fourier transform and Green function methods and obtain a closed form formula for the fair strike value of the variance swap. In fact, it will be expressed as elementary functions without any involvement of integral.
Proposition 1
Under the dynamics (2) of the underlying asset price, the solution of (7), \(j=1,2,\ldots ,N\), is given by
where \(\tau =t_{j}-t\) and \(q=r\tau +\log x\) and \(A_{1}\), \(A_{2}\), \(B_{11}\), \(B_{12}\), \(B_{21}\) and \(B_{22}\) are given by
respectively. Here, the functions \(\alpha _{1}\), \(\alpha _{2}\), \(\beta _{1}\) and \(\beta _{2}\) are explicitly given by
respectively.
Proof
Using the Fourier transform \(\hat{\phi }\) defined by
Equation (7) is transformed into
where the operator \(\hat{\mathcal{L}}\) and the function \(\hat{v}\) are given by
respectively. Here, \(v(e^q;I)\) denotes the payoff function.
Let \(\hat{G}(\tau ,w,y,z;I)\) be the solution of the PDE problem
Then it is the Fourier transform of the Green function of Eq. (7). Since the coefficients of \(\hat{\mathcal{L}}\) are linear in y and z, following Heston’s procedure in Heston (1993) for finding a (affine) solution, we suppose that
Putting (13) into (12), we can find the following ordinary differential equation (ODE) problems for A, \(B_1\) and \(B_2\).
Solutions of these Riccati type nonlinear differential equations can be explicitly obtained as
Note that the solutions are expressed as a combination of the single-factor Heston type solutions.
By the convolution theorem, which states that the Fourier transform of a convolution of two functions is the product of their Fourier transforms, one can deduce
Here, the payoff v(x; I) is given by
and thus the Fourier transform of it can be found to be
where \({\delta }_{0}'\) and \({\delta }_{0}''\) are the first and second derivatives of the generalized Dirac delta function, respectively, defined in the distributional sense and satisfy the property
Then from (13) and (15) we obtain
where \(k_1(\tau ,w;q)\) is given by
in terms of A, \(B_1\) and \(B_2\) in (14).
Now, if we define \(A_1\), \(A_2\), \(B_{k1}\) and \(B_{k2}\), \(k=1,2\), as
respectively, then one can find that \(\frac{\partial A}{\partial w}|_{w=0}\) and \(\frac{\partial B}{\partial w}|_{w=0}\) are pure imaginary numbers while \(\frac{\partial ^{2} A}{\partial w^2}|_{w=0}\) and \(\frac{\partial ^{2} B}{\partial w^2}|_{w=0}\) are real numbers and (16) leads to
where \(A_1\), \(A_2\), \(B_{k1}\) and \(B_{k2}\), \(k=1,2\) are given as in the statement of Proposition 1. Further direct calculation of this yields the desired form in Proposition 1. \(\hfill\square \)
Next, we solve the PDE system (8) for the outer expectation \(W^j\).
Proposition 2
Under the dynamics (2) of the underlying asset price, the solution of (8), \(j=1,2,\ldots ,N\), is given by
where \(\varDelta t =t_{j}-t_{j-1}\) and \(A_{1}\), \(A_{2}\), \(B_{11}\), \(B_{12}\), \(B_{21}\) and \(B_{22}\) are given as in Proposition 1and \(f_1\), \(f_2\), \(g_1\) and \(g_2\) are given by
respectively.
Proof
From the jump condition (2) and Proposition 1 with \(t=t_{j-1}\) (so, \(\tau =t_j-t=\varDelta t\)), the terminal condition of \(W^j\) is
which is a polynomial in y and z and independent of x. By the Feynman–Kac theorem, \(W^j\) is also independent of x on the interval \([0,t_{j-1})\) and (8) becomes
Then by the Feynman–Kac theorem, \(W^j\) becomes
Since \(Y_t\) and \(Z_t\) are both CIR processes, we have
Substituting (18) into (19) and exploiting (20), the desired result in Proposition 2 follows by linear property of conditional expectation and the independence of two processes \(Y_t\) and \(Z_t\). \(\hfill\square \)
Using the result in Proposition 2, we finally obtain the fair strike price of a variance swap as follows.
Theorem 1
Under the rescaled double Heston model (2), the fair strike of a variance swap at \(t=0\) defined in terms of the realized variance (3) is given by
where \(W^j(t,y,z)\) is calculated as in Proposition 2.
Proof
From the result (10) and Proposition 2, we obtain the theorem. \(\hfill\square \)
We note that the fair strike price of the variance swap does not depend on the spot price x. It depends only on the variance y of the underlying asset return. That is because the variance swap considered in this work is a vanilla variance swap. However, if it is an exotic variance swap such as gamma swap which has a different payoff structure, then it would depend on the spot price as studied by Kim and Kim (2020).
4 Numerical Results
In this section, we first check the validity of our theoretical formula in Theorem 1 using Monte Carlo (MC) simulation and then compare the rescaled double Heston model (2) (in brief, the rDH model) to the Heston model and the original double Heston model (1) (in brief, the DH model) based on the result of calibration to VIX market data. The theoretical formulas for the fair strike price under the Heston and DH models are given together with calibration results in “Appendix”.
4.1 MC Simulation
The analytic formula in Theorem is already given in a closed form but we confirm the validity of it by comparing it with a MC simulation result. For MC simulation, we use the Milstein method Mil’shtejn (1975) to reduce discretization errors since path dependent derivatives like variance swap may be vulnerable to errors created by discretizing the Brownian motions. Sample paths generated by the Milstein method under the rDH model are given by
where \(\varDelta t=T/N\), \(t_{j}=j\varDelta t\), \(\varDelta W_{i}(t_j)=W_{i}(t_{j})-W_{i}(t_{j-1})\), \(j=1,2,\ldots ,N\), \(i=1,2,3,4\) and the Brownian motions \(W_{i}\) are mutually independent. 2,000,000 number of sample paths are utilized for pricing the variance swap via MC simulation. The parameters are given by \(r=0.02\), \(\alpha =0.4\), \(\beta =0.1\), \(\theta _1=0.03\), \(\theta _2=0.5\), \(\rho _{xy}=-0.8\), \(\rho _{xz}=-0.7\), \(X_0=100\), \(Y_0=0.02\) and \(Z_0=0.04\).
The MC simulation result is given in Table 1. One can see that the formula in Theorem 1 is guaranteed to be derived correctly and the relative error is more reduced when sampling frequency is increased.
4.2 Sensitivity Analysis
In this section, we investigate the sensitivity of the fair strike value to the parameters of the fast and slow variance factors.
Figure 1a and b represent the impact of the long-run mean levels \(\theta _1\) and \(\theta _2\) of the fast and slow variances on the fair strike price under the rDH model, respectively. We set \(\theta _2=0.2\) in Fig. 1a and \(\theta _1=0.2\) in Fig. 1b. The other parameter values in Fig. 1 are fixed as \(r=0.02\), \(\delta =0.3\), \(\epsilon =0.1\), \(\rho _1=-0.8\), \(\rho _2=-0.9\), \(y(0)=0.15\) and \(z(0)=0.12\), where y(0) and z(0) are the initial values of the fast and slow variances, respectively. The figures indicate that the fair strike value increases as \(\theta _1\) or \(\theta _2\) increases as it should be. In Fig. 1a, the fair strike price decreases as tenor becomes shorter when \(\theta _1\) is bigger than or equals to \(y(0)=0.15\) while it decreases and then increases as tenor gets shorter when \(\theta _1\) is smaller than y(0). Meanwhile, in Fig. 1b, the fair strike price decreases or increases and then decreases as tenor becomes shorter depending on whether \(\theta _2 > z(0)\) or not, respectively. Comparing Fig. 1a with b, it can be viewed that the fair strike of the variance swap is more affected by the fast variance factor than the slow variance factor and the slow factor effect increases gradually as time-to-maturity gets longer.
Figure 2a and b show the influence of the mean reversion rates \(1/\alpha \) and \(\beta \) on the fair strike price under the rDH model. We set the parameter values as \(\theta _1=0.1\), \(\theta _2=0.2\), \(y(0)=0.15\) and \(z(0)=0.12\) in Fig. 2. In this case, we note that \(\theta _1 - y(0) < 0\) and \(\theta _2 - z(0) > 0\) hold. One can notice from Fig. 2a that the lower \(\alpha \) (higher mean reversion rate) is, the lower the fair strike price is. It could be explained that the higher mean reversion rate makes the fast variance process \(Y_t\) move to the lower mean level \(\theta _1\) than the initial level y(0) rapidly. On the contrary, as seen in Fig. 2b, the higher mean reversion rate (higher \(\beta \)) results in the higher strike price because the slow variance process \(Z_t\) moves to the higher mean level \(\theta _2\) than the initial level z(0) fast.
4.3 Calibration
There is no market data of variance swap available as it is traded in the over-the-counter market. However, the CBOE VIX term structure could be an alternative to calibrate the parameters. The definition of VIX is given by
Under the rDH model, it is reduced by the Ito formula to
Thus we could use the VIX term structure data as an instrument for calibration, regarding the VIX as a square root of variance. One can find a detailed study of VIX term structure, for example, in Luo and Zhang (2012).
In this paper, we quote weekly VIX term structure data from 2015 to 2019 and from 03/01/2020 to 30/10/2020 in the CBOE website (https://www.cboe.com), making a distinction between before and after the COVID-19 pandemic. For calibration purpose, the data from 2015 to 2019 is divided into five periods of 1 year. We perform a nonlinear regression analysis to estimate the parameters by taking the following steps. First, we estimate y(0) and z(0) for each day with other parameters fixed (Step1). Second, we estimate \(\alpha \), \(\beta \), \(\theta _1\), \(\theta _2\), \(\rho _{xy}\) and \(\rho _{xz}\) with the fixed y(0) and z(0) obtained in Step 1 (Step2). We repeat Step 1 and Step 2 with updated parameters until the mean-square-error (MSE) is minimized, i.e., the procedure stops when
and
where \(\text{MSE}_{k}\) is defined by
Here, K(t) is the theoretical price of the variance swap, \(\text{VIX}(t_i,T_j)\) is the market data of VIX term structure at \(t_i\) with time-to-maturity \(T_j\), M is the total number of data, k is the number of times that the calibration steps are repeated. The calibrations for the three models are implemented based upon MATLAB 2018c and executed on Intel(R) Core(TM) i5-6300HQ.
Table 2 shows the estimated parameters of the rDH model for the periods of 2015, 2016, 2017, 2018, 2019 and 2020. The estimated parameters of the Heston and DH models are quoted in Table 4 in “Appendix”. To show more credibility of the excellence of the rDH model, we also perform a calibration with a bigger sample size of VIX term structure covering the whole period from 01/01/2015 to 30/10/2020. As one can see from Table 2, it is reasonable to say that \(1/\alpha \) could work as a fast mean reversion rate while \(\beta \) as a slow mean reversion rate.
The MSEs and the elapsed time of calibration are quoted in Table 3. In this table, for example, ‘3 m 20 s’ means 3 min and 20 s. One can notice that the MSEs of calibration of the rDH and DH models are much smaller than that of the Heston model. Also, as shown in Fig. 3, in which ‘Difference’ refers to the absolute error (%) between the computed square root of the variance swap value and the VIX data, the rDH model captures the market volatility much better than the Heston model. It is interesting to see in Table 3 that there is not much difference of MSE between the rDH model and the DH model but the calibration time cost of the rDH model is much smaller than that of the DH model. It suggests that the role of vol-of-vol is very limited for model calibration or pricing of variance swaps under the double Heston framework. One noticeable fact is that the MSEs have far increased in the pandemic situation of COVID-19 around the world in 2020 regardless of the model of choice. Even the DH model has a difficulty of catching the market movement if the market is in extremely volatile situation.
Figure 4 presents the fitting results of the Heston, rDH and DH models to the VIX term structures in 2019 and 2020. As illustrated in Fig. 4a and b, in 2019 when the market was relatively stable, the VIX term structure was usually concave and upward and both the rDH and DH models captured the market quite well while the Heston model still followed the market trend but there was a significant gap. However, after the World Health Organization (WHO) declared COVID-19 pandemic on March 11, the situation has been changed drastically. As demonstrated in Fig. 4c–d, the VIX market itself was so irregular and unstable that even the multi-factor stochastic volatility models such as the rDH and DH models could not capture the market behavior.
In order to investigate how much the slow variance factor has effects on the variance swap value with the calibrated parameters in Table 2, we first decompose the fair strike price \(K_{var}\) given by (21) into three terms as follows.
where \(K_y\) is the term that depends only on the fast variance factor and \(K_z\) is the term that depends only on the slow variance factor, and \(R_{yz}\) is the remainder term that depends on both the fast and slow factors. Figure 5 illustrates the graphs of \(K_{var}\), \(K_y\), \(K_z\) and \(R_{yz}\) under the rDH model with the parameters in Table 2 with \(y(0)=0.0083\) and \(z(0)=0.062\) for the year 2020, where \(\theta _2\) is particularly very small as it is given by \(1.08\times 10^{-8}\). The figure shows that the fast and slow terms are both influential but the remainder term is negligible (average at \(1.95\times 10^{-2}\)). Even if the long-run mean level of the slowly mean reverting variance is very small, the slow term still can have a sufficient impact on the price formation of variance swaps.
5 Conclusion
This paper has proposed a computationally more efficient rescaled version of the double Heston model and applied to the pricing problem of variance swaps. We obtain a closed form analytic formula of the fair strike price and investigate the impacts of the fast and slowly mean reverting variance factors on the price. The pricing formula contains neither integral nor terms required to be numerically calculated, making pricing or calibration much easier. We exploit VIX term structure data for calibration and divide the time period into before and after the COVID-19 pandemic started. We find that volatility after the start of the pandemic in 2020 is so irregular that any model cannot cope with the market VIX term structure. However, we verify that the proposed rescaled double Heston model and the original double Heston model are equally superior to the Heston model in terms of fitting to the market data in a stable market condition, while the proposed model is much more efficient than the original double Heston model in view of calibration time cost. Generally speaking, our study in this article supports why the market practitioners had better view the volatility factors from the two different time scale point of view (fast and slow factors).
References
Broadie, M., & Jain, A. (2008). The effect of jumps and discrete sampling on volatility and variance swaps. International Journal of Theoretical and Applied Finance, 11(8), 761–797.
Christoffersen, P., Heston, S. L., & Jacobs, K. (2009). The shape and term structure of the index option smirk: Why multifactor stochastic volatility models work so well. Management Science, 55(12), 1914–1932.
Cox, J. C., Ingersoll, J. E., Jr., & Ross, S. A. (1985). A theory of the term structure of interest rates. Econometrica, 53(2), 385–407.
Fallah, S., & Mehrdoust, F. (2019). American option pricing under double Heston stochastic volatility model: Simulation and strong convergence analysis. Journal of Statistical Computation and Simulation, 89(7), 1322–1339.
Fouque, J.-P., Papanicolaou, G., Sircar, R., & Solna, K. (2003). Multiscale stochastic volatility asymptotics. Multiscale Modeling and Simulation, 2(1), 22–42.
Fouque, J.-P., Papanicolaou, G., Sircar, R., & Solna, K. (2011). Multiscale stochastic volatility for equity, interest rate, and credit derivatives. Cambridge: Cambridge University Press.
Gatheral, J. (2008). Consistent modeling of SPX and VIX options. Bachelier Congress, 37, 39–51.
Gauthier, P., & Possamai, D. (2011). Efficient simulation of the double Heston model. The IUP Journal of Computational Mathematics, 4(3), 23–73.
Grzelak, L. A., & Oosterlee, C. W. (2011). On the Heston model with stochastic interest rates. SIAM Journal on Financial Mathematics, 2(1), 255–286.
Heston, S. L. (1993). A closed form solution for options with stochastic volatility with applications to bond and currency options. The Review of Financial Studies, 6(2), 327–343.
Kim, S.-W., & Kim, J.-H. (2020). Pricing generalized variance swaps under the Heston model with stochastic interest rates. Mathematics and Computers in Simulation, 168, 1–27.
Little, T., & Pant, V. (2001). A finite difference method for the valuation of variance swaps. Journal of Computational Finance, 5(1), 81–103.
Luo, X., & Zhang, J. E. (2012). The term structure of VIX. Journal of Futures Markets, 32(12), 1092–1123.
Mil’shtejn, G. N. (1975). Approximate integration of stochastic differential equations. Theory of Probability and its Applications, 19(3), 557–562.
Oksendal, B. (2000). Stochastic differential equations: An introduction with applications. Berlin: Springer.
Swishchuk, A. (2004). Modeling of variance and volatility swaps for financial markets with stochastic volatilities. Wilmott Magazine, 2, 64–72.
Wilmott, P. (2013). Paul Wilmott on quantitative finance. Hoboken: Wiley.
Zhang, S. M., & Feng, Y. (2019). American option pricing under the double Heston model based on asymptotic expansion. Quantitative Finance, 19(2), 211–226.
Zheng, W., & Kwok, Y. K. (2014). Closed form pricing formulas for discretely sampled generalized variance swaps. Mathematical Finance, 24(4), 855–881.
Zheng, W., & Kwok, Y. K. (2014). Saddlepoint approximation methods for pricing derivatives on discrete realized variance. Applied Mathematical Finance, 21(1), 1–31.
Zhu, S.-P., & Lian, G.-H. (2011). A closed-form exact solution for pricing variance swaps with stochastic volatility. Mathematical Finance, 21(2), 233–256.
Acknowledgements
We thank the anonymous referees for their careful reading of our earlier version of the manuscript and their insightful comments and suggestions. The research of J.-H. Kim was supported by the National Research Foundation of Korea NRF2021R1A2C1004080.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
We declare that we have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix A
Appendix A
1.1 A.1 Fair Strike Price Formulas under the Heston and DH Models
The Heston model is known as
Using the same method as in the derivation of Theorem 1, the fair strike, \(K_{\text{H}}\), of a variance swap under the Heston model can be obtained analogously as follows.
where \(C_0\), \(C_1\) and \(C_2\) are given by
respectively, Here, \(\alpha _1\), \(\alpha _2\), \(\beta _1\) and \(\beta _2\) are the same as in Proposition 1. Note that this formula is based on the realized variance defined by log return, which is practically used in financial industry, while the formula of Zhu and Lian (2011) underlies simple (arithmetic) return.
Following the similar procedure to the derivation of Theorem 1, the fair strike, \( K_{\text{DH}}\), of a variance swap under the original double Heston model (1) can be obtained as
where \(\varDelta t =t_j - t_{j-1}\) and \(\tilde{f}_1\), \(\tilde{f}_2\), \(\tilde{g}_1\), \(\tilde{g}_2\), \(\tilde{A}_{1}\), \(\tilde{A}_{2}\), \(\tilde{B}_{11}\), \(\tilde{B}_{12}\), \(\tilde{B}_{21}\) and \(\tilde{B}_{22}\) are given by
respectively. Here, \(\alpha _1\), \(\alpha _2\), \(\beta _1\) and \(\beta _2\) are the same functions as in Proposition 1.
1.2 A.2 Calibration Parameter Results of the Heston and DH Models
See the Table 4.
Rights and permissions
About this article
Cite this article
Yoon, Y., Kim, JH. A Closed Form Solution for Pricing Variance Swaps Under the Rescaled Double Heston Model. Comput Econ 61, 429–450 (2023). https://doi.org/10.1007/s10614-021-10214-6
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10614-021-10214-6