1. Introduction
With the globalization of economies and intensifying competition in financial markets, more and more investors and scholars have realized the importance of hedging against volatility risk. Against this background, volatility derivatives have emerged as instruments that hedge against this volatility risk. Their most prominent feature is that they provide direct exposure to asset volatility without continuous delta-hedging [
1]. Among the various volatility derivatives, variance or volatility swaps are considered the main and essential ones. Essentially, they are forward contracts for future realized volatility or variance. The long position of a volatility or variance swap pays a fixed delivery price at expiry and receives the floating amounts of annualized realized volatility or variance, whereas the short position is just the opposite [
2]. In other words, the payoff of a variance swap at maturity
T is
, where
is the realized annual variance of the return,
K is the strike price such that the initial value of the variance swap equal to zero for both long and short positions, and
I is the notional amount of the swap in dollars per annualized volatility point squared.
In the literature, there are roughly two types of approaches for pricing volatility swaps or variance swaps: static replication methods and model-based stochastic methods. Static replication means that a variance swap could be replicated theoretically by using a portfolio of standard options and without specifying the volatility function. For example, Demeterfi [
1] produced some of the most influential pioneering work that demonstrated how to use a standard option portfolio to replicate variance swaps. Obviously, the outstanding advantage of static replication is that it does not require a specific model to solve pricing problems [
3], and this method has been widely used for the VIX index since 2003 [
4]. However, according to Sahalia et al. [
5], this method has drawbacks that make it unsuitable for the real market. These include the inability to promptly reflect jumping variance swap rates, as well as the need for constant adjustment in options. Therefore, we consider model-based stochastic methods in this paper.
The study of variance swaps pricing in stochastic models can be broadly classified into two approaches. The first method entails pricing under stochastic volatility with continuously sampled realized variance. It provides an effective way to reduce the computational difficulty, and recent work can be found in [
6,
7,
8,
9,
10]. However, this approach is an idealization that assumes continuous sampling time. According to the research conducted by Zhu and Lian [
11], using this method to estimate realized variance may lead to significant inaccuracies in variance swaps with longer tenors or less frequent sampling rates when compared to discretely sampled variance swaps traded in the OTC market. Fortunately, the second method was developed to overcome this limitation, which entails pricing the realized variance obtained directly from discrete sampling. Broadie and Jain [
12] used discretely sampled data and various stochastic models to examine the fair strike prices of variance swaps and volatility swaps. They also compared their pricing formulae with a continuously sampled situation and found that as observation frequency increases towards infinity, the discrete price converges to the continuous price. Rujivan and Zhu [
13] proposed a simplified analytical approach for pricing discretely sampled variance swaps when the realized variance is defined as the squared log return of the underlying price. Zhang [
14] obtained an analytical solution for variance swaps under the MRG stochastic volatility model. Lian [
15] first explained the conditions for discrete volatility swaps to exhibit convergence behavior. Recent work in this field can be seen in [
16,
17,
18].
However, the regime-switching model is rarely considered in most of the literature cited above. It is worth noting that the long-term trends in the financial market, which are influenced by structural changes in macroeconomic conditions, can have a significant impact on prices. That is, incorporating market trends into regime considerations can help convert market volatility uncertainty into regime switching to some extent, and it can also better assist investors in capturing market risks and opportunities [
19]. As Elliott and Lian [
20] pointed out, different market states (good or bad economy) can be reflected by different regimes, and the parameters of the regime-switching model are modulated by an observable Markov chain. In other words, the regime-switching model can be described as a Markov-modulated model. Because the market is generally incomplete, the Esscher transform has been extensively studied as an effective way to determine the unique equivalent martingale measure (see [
21,
22,
23]). To the best of our knowledge, although Markov-modulated Heston models have been widely applied and studied for pricing options [
4,
24,
25,
26,
27], the pricing of discrete-sampling variance swaps under an Markov-modulated MRG model has not yet been investigated. Compared to the Heston model and its reformulated version, the MRG volatility model has the following highlighted strengths. First, it has the mean-reverting property that ensures the autocorrelation of volatility [
28] when volatility clustering occurs. Second, it has a relatively explicit solution for pricing. Therefore, in this paper, we attempt to study the pricing of discrete-sampling variance swaps under the Markov-modulated MRG model for the first time. This research has significant meaning and value. It can enrich the existing theory of regime-switching models and help market participants better understand market trends.
At the same time, there is a part of the regime-switching model that has not been well explored, which is the underlying asset process with jump diffusion. Sudden and unexpected events in real-world financial markets can result in significant fluctuations in asset prices, which cannot be disregarded, especially when pricing volatility derivatives. Many researchers have examined the jump effect with numerical tests and have proven that it can increase the degree of matching to the reality of financial markets [
29]. In particular, when the jump diffusion follows the asymmetric double exponential distribution proposed by Kou [
30], the model can help capture market-implied volatility smiles and market jump behavior [
31]. However, little work has been done on pricing variance swaps under a regime-switching model with asymmetric double exponential jump diffusion, especially the Markov-modulated MRG model.
Motivated by the reasons stated above, we introduce a novel model in this paper—a Markov-modulated MRG model with asymmetric double exponential jump diffusion—to price discretely sampled variance swaps. To our knowledge, this is the first attempt to apply this new framework to pricing discretely sampled variance swaps. Hopefully, this attempt to price variance swaps under this new model could provide some thinking for other researchers and fill the research gap. By utilizing the Feynman–Kac theorem and transition probability function, we can greatly simplify the pricing process and theoretically derive the general expression of the characteristic function. This in turn allows us to solve the pricing problem more effectively. The main contributions of this paper can be summarized as follows: First, one potential issue for the regime-switching model that incorporates various economic factors is that the state space of the Markov chain can become quite large. This can lead to significant computational complexity in pricing and may pose challenges in practice. For practical purposes, it is significant to reduce the additional burden caused by the Markov chain. In this paper, we creatively obtain analytic pricing formulae that can reduce the computational burden and give the analytical solution of the ordinary differential equations satisfied by the Feynman–Kac theorem. To achieve this, we propose a new method of dividing the sampling time into more detailed switching intervals and non-switching intervals. Our Proposition 2 presents a specific expression for this solution. Second, we have made a creative improvement to the semi-Monte Carlo simulation method that verifies the accuracy and efficiency of our pricing formulae. Specifically, we use independent and identically exponentially distributed random numbers to generate repeated simulation paths randomly, as opposed to the traditional method of generating repeated simulation paths. Hopefully, it can be viewed as a step forward for the semi-Monte Carlo simulation in the regime-switching field.
2. Model Setup
In this section, we will formally introduce the MRG model with regime-switching jump diffusion and the computation methods for calculating realized variance. Additionally, we will discuss the specific details of the regime-switching framework.
2.1. The Regime-Switching Frame and Markov Chain
Regime-switching models have generated considerable interest for their applications to a variety of financial problems. As the trend of the financial market significantly influences the movement of individual stocks, it is essential to enable the key parameters of the stock to react to general market fluctuations. The market regime is one possible way to reflect the state of the economy, the mood of investors, and other economically realistic factors. If the Markov chain is observable, then the state can be interpreted as the market regime [
20,
25,
32]. In other words, we can let parameters of the stochastic modes be modulated by the Markov chain, then the market factors can be considered in the model. It can help researchers capture the market movement by representing different market regimes.
Here, we follow the notation introduced in [
4]:
is a continuous-time finite-state
observable Markov chain on a complete probability
, where
presents the real-world probability measure,
represents
N different regimes of the financial market and the
ith state
. Moreover, the dynamic of
X is given by
where
is an
-valued martingale increment process with respect to the natural filtration generated by
X, and
is a transition matrix. The element
presents the transition rate from state
to
, satisfying
and
(it shows that the row sum of matrix
is always zero).
2.2. The MRG Model with Regime-Switching Jump Diffusion
Under the Markov-modulated MRG model with double exponential jump diffusion, the dynamic of the stock price process
with a short-term volatility process
is given by
All variables here are set on unless otherwise specified, where presents the real-world probability measure. Assume is the inner product and for every , the elements in are non-negative.
In the first equation of (
2),
is the appreciation rate;
is a standard Brownian motion;
is a Poisson process with rate
; and
is a sequence of independent identically distributed (i.i.d) nonnegative random variables such that
has an asymmetric double exponential distribution with the density
where
represent the probabilities of upward and downward jumps.
In the second equation of (
2),
is the long-term average value of volatility. To simplify our calculation and deduce an analytic solution for pricing variance swaps, we regard the volatility
and the mean recovery speed parameter
k as constants;
is another standard Brownian such that
.
Let the dynamic of the bank account process
follow
where
is the instantaneous market interest rate.
For simplicity, we also let , , and be the natural filtrations generated by , , and Markov chain , respectively.
Remark 1. This is the new model that we propose. We take the appreciation rate, long-term average value of volatility, and instantaneous market interest rate to be Markov-modulated in a creative way. In other words, the values of these parameters are not fixed compared to the traditional MRG model. Additionally, we add an asymmetric double exponential jump-diffusion to the asset process and also let the jump intensity be Markov-modulated parameters.
2.3. Variance Swaps
The realized variance in variance swaps is typically calculated through discrete sampling over the period
in financial markets, where
T represents maturity time. Two common methods for calculating realized variance are actual-return realized variance
and log-return realized variance
where
denotes the
kth observation of the stock price at time
, which implies that the total number of observations is
N. In addition,
is an annualized factor with a range from 4 to 252, which can be formed as
.
4. Pricing Variance Swaps
In this section, we will take the distribution characteristics of the new model as a breakthrough, and derive the pricing formulae of variance swaps under discrete sampling. Here, we employ the characteristic function method to accomplish this goal.
Now we are in a position to formally introduce the characteristic function method. Let
, and assume its characteristic function is given by
It should be noted that the imaginary number
is absorbed in the parameter
in (
14), and the result will not be changed without the explicit expression of the imaginary number [
20]. Similar to the method used in [
4], we first consider the conditional characteristic function under the given Markov chain information
. When the information
is given, the Markov-modulated parameters
,
, and
are deterministic functions, which means that the information available at time
is sufficient to evaluate these parameters. Namely, we can simplify (
14) as
4.1. Conditional Characteristic Function
In this section, we will focus on calculating the concrete expression of conditional characteristic function and give the following lemma and proposition.
Lemma 1. If the volatility process satisfies the dynamics (8), then it has an explicit expression as follows: Proof. Let
and
For the Ito–Doeblin formula, we shall need the following the partial derivatives of
:
We will also need the partial derivative of
, which is
Then, the Ito–Doeblin formula states that
which shows that the dynamics (
8) satisfies the expression (
16). □
Proposition 1. Given the information , if the stochastic process follows (8) and the parameters satisfy , then the characteristic function iswithwhere satisfies the following ordinary differential equations: Proof. The proof of this proposition is presented in
Appendix A. □
4.2. Characteristic Function
In this section, we will consider the effect of Markov chain
on the characteristic function. Using Formula (
17) in the previous section, we can derive that
which means the core calculation of characteristic function is
.
Note that
is defined by (
18) and
satisfies the ordinary differential equations in (
19), then we can rewrite
as
In addition, if we try to derive the partial derivatives of (
21), then
where
Therefore, the partial derivatives of
satisfy
In other words,
can be rewritten as
According to the Proposition 3.2 in [
20], if
denotes the transposition of
Q, then
where
and
Therefore, we can calculate (
24) in theory if
and
. So far, the theoretical derivation of the characteristic function has been completed in this paper. In the next section, we will first give the analytical solution of the ordinary differential equations in (
19) that ensure high calculation precision.
4.3. Derivation of ODEs
In this section, we will focus on deriving the analytical solutions of (
19) in Proposition 1. These solutions can greatly reduce errors caused by directly simulating numerical solutions of differential equations and improve computational efficiency.
Proposition 2. For any , if function follows the ordinary differential equationthen we can let denote the times of regime switching in , and obtain the following analytical solutions of and :where Proof. The proof of this proposition is presented in
Appendix B. □
Remark 2. The basic idea of this proposition is to simplify the numerical solution process of the differential equation. Considering that the Markov-modulated parameters remain constant before the next regime-switching moment, this proposition creatively divides the integration interval of the ordinary differential equation into the regime-switching interval and non-regime-switching one. Therefore, one can solve the simple function definite integral within these small intervals.
4.4. Pricing Formulae
In this section, we will use the relationship between the characteristic function and the expectation to derive the pricing formulae. Recall that the fair strike price
of variance swaps should satisfy
. In other words, we can obtain that
Similarly, we can compute
where
.
5. Numerical Tests and Discussion
In this section, we will use semi-Monte Carlo simulation to verify the accuracy and efficiency of our pricing formulae. It provides direct numerical evidence of the comparison between the Monte Carlo simulation of the model and our formulae. We will also examine the rationality of the promotion term in the new model, as it will help readers understand the effects of these realistic economic factors. In addition, we will analyze the sensitivity of (
28) and (
29) by studying a
change of model parameters. This analysis not only helps to assess the stability and reliability of the model, but also allows us to identify the parameters that exert the greatest influence on its results.
We assume a typical market with only two regimes, namely, the market state switches between regime 1 and regime 2. These market states can be interpreted as bull markets (a good economy) and bear markets (a bad economy), respectively. Thus, the state space in
is
, and the transition rate matrix is given by
Moreover, we will adopt the values of parameters in
Table 1 that are used in [
4,
14].
5.1. Semi-Monte Carlo Stimulation
In this section, we will test the accuracy and efficiency of our pricing Formulae (
28) and (
29) by comparing the results generated by semi-Monte Carlo simulation. This comparison more directly and effectively demonstrates the rationality of the new model and the advantages of the pricing formulae.
We utilize Euler–Maruyama discretization for our model (
8) to Monte Carlo simulation:
where
and
are two independent standard normal random variables. If we divide the sampling interval
into
n sub-intervals, then the time length of each part satisfies
.
Compared to traditional Monte Carlo simulation, semi-Monte Carlo simulation is an effective method for improving the simulation efficiency of regime-switching models [
4,
20,
38]. In traditional Monte Carlo simulation, randomness exists not only in the simulation of regime-switching paths but also in the simulation of asset price and volatility processes. This means that two identical regime paths can have different calculation results. In contrast, the randomness of semi-Monte Carlo simulation only exists in the regime-switching path. Therefore, if the paths are the same, the fair price of the simulation must be the same. Given that semi-Monte Carlo simulations require fewer computation paths than Monte Carlo simulations, this paper leverages independent and identically exponentially distributed random variables to effectively generate simulation paths and capture the necessary randomness. Specifically, we will adopt the following four steps to implement the semi-Monte Carlo simulation process:
Step 1: Generate N = 20,000 independent identically exponentially distributed (i.i.d) random variables with parameter . These variables represent the time until the first regime-switching event occurs for each of the N sample paths. Based on the initial market regime and the probability of remaining in this regime , we can estimate the number of sample paths that do not experience a regime switch before the maturity time as . Note that the actual number of such paths may differ due to randomness in the generation of the i.i.d random variables. Therefore, we only need to run the following computation for sample paths, where is the number of paths that do experience a regime switch before maturity, and represents the one path that does not experience any regime switches.
Step 2: For every sample path
, calculate the conditional characteristic function (
17) according to Propositions 1 and 2. When there is no switch in the regime switching path, it is worth noting that the difference between the analytical solutions of the differential equation system (Proposition 2) only depends on the interval length and the initial regime at the starting point (keeping in mind that there may have been a switch before the path interval, and the initial regime may be different). Additionally, the difference in the initial regime is proportional. Therefore, before starting step 2, we assume that there are no regime switches in the whole process. Then, we calculate and save the results under different sampling frequencies. After entering the loop simulation, for any path, we divide it into intervals with and without switching based on the sampling frequency. In the subsequent calculations, we call and adjust these results when there is no switching in the initial regime.
Step 3: Utilize the conditional characteristic function in step 2 to calculate the fair strike price , .
Step 4: Calculate the mean of the fair strike price
as the ultimate fair strike price
, that is
Moreover, we provide an algorithm (Algorithm 1) that calculates the fair strike price for a variance swap. Here, we only take as an example.
Figure 1 shows images of the fair strike price obtained from semi-Monte Carlo simulation of Equation (
28) and those calculated by Monte Carlo simulation under different observation frequencies ranging from 4 to 252. Actually, the observation frequencies represent the number of trading days in a year, with 4 and 252 representing quarterly and daily observations, respectively.
Algorithm 1: Calculation of the strike price |
|
It is evident that our solutions match the simulations well, which effectively verifies the accuracy of our pricing formulae. Furthermore, our pricing formula has higher computational efficiency, with a simulation time of 11.30401 s. In contrast, the Monte Carlo simulation time for the model requires 3.962049 min (using the simulation running time of as an example).
5.2. Regime-Switching Effect
In this section, we will compare the fair strike price
computed from the regime-switching hybrid model and the non-regime-switching hybrid model, as shown in
Figure 2. Here, the
without regime switching mainly refer to the situation where the model parameters are kept constant at the initial regime value throughout the entire contract period
T. For example, if the initial regime is bearish, then the model parameters during the maturity
T are all the values under regime 2.
Figure 2 shows that the price
can be significantly impacted by regime switching. More specifically, the regime tends to lower
to some level when it switches from a bull to a bear regime. On the other hand, when the regime switches from a bear to a bull regime, the price tends to rise somewhat. At the same time, we also observe that
are relatively higher when the initial regime is in a bull state (regime 1) compared to when the market regime is in a bear state (regime 2), which is consistent with the findings in [
4] and real financial markets. Particularly, when the observation frequencies rises from quarterly to daily, we keep track of the price under regime 1 (see
Table 2).
Furthermore, we can see from
Table 2 that
is considerably higher than
, and regime switching can result in large variations in
. This is mostly due to the fact that (
6) prioritizes describing the relative volatility of asset prices, whereas (
5) focuses more on the intensity of asset price swings.
5.3. Jump Diffusion Effect
In this section, we will investigate the impact of jump diffusion by comparing the fair strike prices
with different jump intensities
. Unless otherwise specified, we assume that there is no jump effect in the stock price process when
. Since the major concern in this section is not the regime-switching effect, we choose to explore the jump diffusion effect under regime 1.
Figure 3 shows
at different observation frequencies as the jump intensity
increases from 0 to 0.4.
From
Figure 3, we can observe that at each observation frequency, the jump in
leads to higher prices, and as the jump intensity
increases, the prices further increase. One possible explanation for this is that jump diffusion may lead to changes in the underlying asset prices, resulting in higher realized variances and higher prices.
5.4. Sensitivity Analysis
In this section, we will examine the sensitivity of our model parameters by computing the fair strike prices
on the
change in it. As regime switching and observation frequency are not the main concern in this section, we consider the sensitivity of parameters under daily sampling and the initial state of regime 1.
Table 3 shows the percentage change in
caused by changes in model parameters.
From
Table 3, we can see that the three parameters—long-term volatility mean
, initial volatility
, and mean reversion speed of volatility
k—have a significant impact on the volatility of fair execution prices
. Furthermore, changes in
k have the most significant impact on the volatility of prices. This is consistent with the real financial market, where asset price volatility is primarily related to the stochastic volatility process of the asset price. Additionally, we find that the initial stock price
has no impact on the volatility of fair execution price
.
5.5. Comparison of the New Model with the Traditional MRG Model
In this section, we will consider the fair strike price under different models and examine how the inclusion of more realistic factors influences price. In other words, we will compare the fair strike price under the new model illustrated in (
8) with prices derived from the traditional MRG model.
Both the Heston model and the MRG model are classic models in the stochastic field, and it would be comprehensive to compare our new model with both of them. However, as pointed out by Zhang [
14], the connection between the MRG model and the Heston model can be made clear under a specific assumption—that the long-term average value of volatility
is equal to zero. Therefore, in this section, we will focus on comparing our new model with the traditional MRG model only. Recall that the traditional MRG models assume a continuous price process with stochastic volatility, and all the parameters in the models are constant. To maintain this setting, we suppose that the parameters in the traditional models are all set in regime 1, and non-regime switching occurs in maturity.
It should be noted that MRG model may not adequately capture the long-term trends in the financial market that are influenced by structural changes in macroeconomic conditions. However, empirical finance literature extensively documents the correlation between the financial market and the macroeconomic environment. Moreover, sudden and unexpected events in real-world financial markets can result in significant fluctuations in asset prices, which cannot be disregarded, especially when pricing volatility derivatives. Therefore, adding the assumption in our new model can make it more comprehensive and meaningful.
Figure 4 shows the strike prices obtained from Zhang’s [
14] formula and ours. (Here, we only take
as an example to show the difference.)
From
Figure 4, we can see that the price under our new model is cheaper than the one in the traditional MRG model. This can be explained as follows. First, both regime switching and jump diffusion contribute to the variation of the strike price, as discussed in the previous subsection. Second, this can be interpreted by the assumption that the parameters in the traditional MRG model are all set in regime 1 (good economy), which is larger than the ones in regime 2 (bad economy). In other words, a bearish market can drag down the price to some extent. Moreover, the small value of the long-term volatility mean may lead to a lower jump probability, meaning that the regime-switching effect has a greater effect on the price than the jump-diffusion effect. The sensitivity analysis in
Table 3 can provide some evidence.
6. Conclusions
In this paper, we investigate the fair strike prices for discretely sampled variance swaps under the Markov-modulated MRG model with asymmetric double exponential jump diffusion. Compared to recent work in the variance swaps pricing field, the major contribution of this paper is to propose an effective approach to obtain a general pricing formula and improve the computational efficiency of semi-Monte Carlo simulation. The pricing problem can be reduced to an expectation calculation under the risk-neutral T-forward measure. We obtain analytical solutions for the pricing problem using the characteristic function and the conditional transition characteristic function. Additionally, we carry out a series of semi-Monte Carlo simulations to verify the accuracy and effectiveness of the pricing formulae. Our experiments aim to explore the influence of regime-switching and jump-diffusion terms, as well as analyze the sensitivity of pricing formulae parameters to the fair strike price. We find that when regime switching occurs, either the fair strike price for good economics will decrease or the fair strike price for bad economics will increase. When there is a jump in market asset prices, the diffusion of the jump will accelerate changes in stock prices, resulting in higher fair execution prices. In addition, we also find that the stochastic fluctuation parameters of the propagation model are more sensitive to price fluctuations.
Furthermore, we would like to emphasize the practical applications of our model. The pricing method proposed in our proposition is generalizable and provides a new way to accurately compute Markov-modulated variables. This approach can also be extended to other regime-switching stochastic models. For example, it could simplify the pricing formulae presented in Elliott [
25] and Yang [
4] without the need for overly complex expressions. Our attempt to price variance swaps under this new model can provide valuable insights and numerical evidence for financial institutions and investors to better manage and control market volatility risk, ultimately reducing trading losses.
However, it is a pity that we do not deeply consider the effect of a long contract life where the interest rate can change stochastically. In other words, if the interest rate is assumed to follow a specific stochastic process instead of a deterministic one, the ultimate risk-neutral measure and pricing formulae may differ. It would be worthwhile to explore how to adjust the pricing approach in such cases for future studies. Moreover, it would be worthwhile to explore reasonable methods for calibrating model parameter values, which we lack in this paper and hope that future research can address.