Next Article in Journal
Does Fintech-Driven Inclusive Finance Induce Bank Profitability? Empirical Evidence from Developing Countries
Next Article in Special Issue
Relations among Bitcoin Futures, Bitcoin Spot, Investor Attention, and Sentiment
Previous Article in Journal
Comparison of the Asymmetric Relationship between Bitcoin and Gold, Crude Oil, and the U.S. Dollar before and after the COVID-19 Outbreak
Previous Article in Special Issue
Interconnectedness of Cryptocurrency Uncertainty Indices with Returns and Volatility in Financial Assets during COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pricing Path-Dependent Options under Stochastic Volatility via Mellin Transform

Department of Mathematical Sciences, School of Engineering, Computer and Mathematical Sciences, Auckland University of Technology, Private Bag 92006, Auckland 1142, New Zealand
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2023, 16(10), 456; https://doi.org/10.3390/jrfm16100456
Submission received: 15 September 2023 / Revised: 18 October 2023 / Accepted: 19 October 2023 / Published: 20 October 2023
(This article belongs to the Special Issue Featured Papers in Mathematics and Finance)

Abstract

:
In this paper, we derive closed-form formulas of first-order approximation for down-and-out barrier and floating strike lookback put option prices under a stochastic volatility model using an asymptotic approach. To find the explicit closed-form formulas for the zero-order term and the first-order correction term, we use Mellin transform. We also conduct a sensitivity analysis on these formulas, and compare the option prices calculated by them with those generated by Monte-Carlo simulation.

1. Introduction

A standard option gives its owner the right to buy (or sell) some underlying asset in the future for a fixed price. Call options confer the right to buy the asset, while put options confer the right to sell the asset. Path-dependent options represent extensions of this concept. For example, a lookback call option confers the right to buy an asset at its minimum price over some time period. A barrier option resembles a standard option except that the payoff also depends on whether or not the asset price crosses a certain barrier level during the option’s life. Lookback and barrier options are two of the most popular types of path-dependent options
Following the lead set by Black and Scholes (1973) and assuming that the underlying asset price follows a geometric Brownian motion with constant volatility, Merton (1973) derived a closed-form pricing formula for down-and-out call options. Reiner and Rubinstein (1991) extended Merton’s results to other types of barrier options. Goldman et al. (1979) and Conze and Vishwanathan (1991) provided closed-form pricing formulas for lookback options. For a good summary for research on path-dependent options under the Black–Scholes framework, refer to Clewlow et al. (1994). As we know, the assumption that an asset price process follows a geometric Brownian motion with constant volatility does not capture the empirical observations, due to the volatility smile effect. So, it is desirable to overcome this drawback. There are different ways of extending the Black–Scholes model to incorporate the “smile” feature: one way is to consider “local volatility”, and the other is to consider “stochastic volatility”.
One popular local volatility model is the constant elasticity of variance (CEV) model introduced by Cox (1975, 1996), where a closed-form pricing formula for European call options was presented. Davydov and Linetsky (2001) derived solutions for barrier and lookback option prices under the CEV process in closed form and demonstrated that barrier and lookback option prices and hedge ratios under the CEV process can deviate dramatically from the lognormal values. In Boyle and Tian (1999), the pricing of certain path-dependent options was re-examined when the underlying asset follows the CEV diffusion process, by approximating the CEV process using a trinomial method.
Heston (1993) assumes that volatility reverts to a long-term mean at a specified rate. Bates (1996) builds upon the Heston model by introducing a jump component for asset prices, which is represented as a compound Poisson process with normally distributed jumps. In a further refinement of the Heston model, the jumps are characterized by infinite activity jumps generated by a tempered stable process, as demonstrated in Zaevski et al. (2014). Despite these advancements, the pricing challenges associated with path-dependent options in the context of stochastic volatility persist, as no analytical solutions are available for these models.
Chiarella et al. (2012) considered the problem of numerically evaluating barrier option prices when the underlying dynamics are driven by the Heston stochastic volatility model and developed a method of lines approach to evaluate the price as well as the delta and gamma of the option. Park and Kim (2013) investigated a semi-analytic pricing method for lookback options in a general stochastic volatility framework. The resultant formula is well connected to the Black–Scholes price that is the first term of a series expansion, which makes computing the option prices relatively efficient. Furthermore, a convergence condition for the expansion was provided with an error bound. Leung (2013) and Wirtu et al. (2017) derived an analytic pricing formula for floating strike lookback options under the Heston model by means of the homotopy analysis method. The price is given by an infinite series whose value can be determined once an initial term is given well.
In addition, Kato et al. (2013) derived a new semi-closed-form approximation formula for pricing an up-and-out barrier option under a certain type of stochastic volatility model, including an SABR model. In a more recent paper by Funahashi and Higuchi (2018), a unified approximation scheme was proposed for a single-barrier option under local volatility models, stochastic volatility models, and their combinations. The basic idea of their approximation is to mimic a target underlying an asset process using a polynomial of the Wiener process. They then translated the problem of solving the first hit probability of the asset price into the problem of solving that of a Wiener process whose distribution of the passage time is known. Finally, utilizing Girsanov’s theorem and the reflection principle, they showed that single-barrier option prices can be approximated in a closed form.
The main contribution of this paper is to derive new closed-form approximation formulas for pricing down-and-out put barrier options and floating strike lookback put options under a certain type of stochastic volatility model, which is similar to the one in Cao et al. (2023); Kato et al. (2013); Kim et al. (2023). To achieve our goal, we apply the asymptotic approach discussed in Fouque et al. (2011) and Mellin transform. Mellin transform techniques were used by Panini and Srivastav (2004) to derive integral equation representations for the price of European and American basket put options. Similarly, Yoon (2014) applied Mellin transform to derive a closed-form solution of the option price with respect to a European call option and a European put option with the Hull–White stochastic interest rate. Moreover, Kim and Yoon (2018) derived a closed-form formula of a second-order approximation for a European corrected option price under a stochastic elasticity of variance (SEV) model.
The rest of the paper is organized as follows. Section 2 discusses the model framework and the features of down-and-out and floating strike lookback put options. In Section 3, we provide detailed discussions on an asymptotic approach, which is used to derive approximations to the risk-netural values of these types of options. In Section 4, we apply Mellin transform to derive a closed-form formula of the first-order approximation for down-and-out barrier put options. In Section 5, we apply Mellin transform to derive a closed-form formula of the first-order approximation for floating strike lookback put options. Section 6 presents a sensitivity and comparison analysis and demonstrates that the results given by these closed-form formulas match well with those generated by Monte-Carlo simulation. Section 7 gives a brief summary. Details on Mellin transform and the derivation of the closed-form formulas in Section 4 and Section 5 are provided in Appendix A and Appendix B, respectively.

2. Basic Model Set-Up and Path-Dependent Options

2.1. Stochastic Volatility Model

Let { S t : t 0 } denote the price process of a risky asset on some filtered probability space Ω , F , ( F t ) t 0 , P , where P is the physical probability measure In this paper, we assume that { S t : t 0 } evolves according to the following system of stochastic differential equations:
d S t = μ S t d t + f Y t S t d W t s , d Y t = α m Y t d t + β ρ d W t s + 1 ρ 2 d W t y ,
where μ , α > 0 , β > 0 , and m are constants and f is a function having positive values and specifying the dependence on the hidden process { Y t : t 0 } . The processes W t s : t 0 and W t y : t 0 are independent standard Brownian motions. The constant correlation coefficient ρ with 1 < ρ < 1 captures the leverage effect. Here, μ is the drift rate. The mean-reversion process { Y t : t 0 } given in Equation (1) is characterized by its typical time to return back to the mean level m of its long-run distribution. The parameter α determines the speed of mean-reversion, and β controls the volatility of { Y t : t 0 } . In the sequel, we shall refer to the above system as the stochastic volatility (SV) model. In Section 2 and Section 3, we will not specify the concrete form of f, but assume that f is bounded and smooth enough, e.g., f C 0 2 ( R ) . Furthermore, f has to satisfy a sufficient growth condition in order to avoid bad behavior such as the non-existence of moments of { S t : t 0 } . For numerical results in Section 6, we choose f to take a special form, as used in Fouque et al. (2000, 2011) and Cao et al. (2021).
We apply the well-known Girsanov theorem to change the physical measure P to a risk-neutral martingale measure Q by letting
d W t s * = μ r f Y t d t + d W t s y and d W t y * = ξ Y t d t + d W t y ,
where ξ Y t represents the premium of volatility risk. Then, the model equations under the measure Q can be written as
d S t = r S t d t + f ( Y t ) S t d W t s * , d Y t = α m Y t β ρ μ r f ( Y t ) + ξ ( Y t ) 1 ρ 2 d t + β ρ d W t s * + 1 ρ 2 d W t y * .
Note that W t s * : t 0 and W t y * : t 0 are independent standard Brownian motions under Q . As an Ornstein–Uhlenbeck (OU) process, { Y t : t 0 } in Equation (1) has an invariant distribution, which is normal with mean m and variance β 2 / 2 α . Thus, we can expect that if mean reversion is very fast, i.e., α goes to infinity, the process { S t : t 0 } should be close to a geometric Brownian motion. This means that if mean reversion is extremely fast, then the model of Black and Scholes would become a good approximation. In reality, however, it may not be the case. For fast but not extremely fast mean-reversion, the Black–Scholes model needs to be corrected to account for the random characteristics of the volatility of a risky asset. For this purpose, we introduce another small parameter ϵ defined by ϵ = 1 / α , as performed by Fouque et al. (2000). For notational convenience, we put ν = β / 2 α . With the help of these notations, the model equations under Q are re-written as
d S t = r S t d t + f Y t S t d W t s * , d Y t = 1 ϵ m Y t 2 ν ϵ Λ Y t d t + 2 ν ϵ ρ d W t s * + 1 ρ 2 d W t y * ,
where Λ ( · ) , defined by
Λ ( y ) : = ρ μ r f ( y ) + ξ ( y ) 1 ρ 2 ,
is the combined market price of risk.

2.2. Path-Dependent Options

Let V ( T ) denote the payoff of a put option on the risky asset at its expiration T. Then, its risk-neutral price at time t [ 0 , T ] under our SV model is given by
P t , s , y = E Q e r T t V ( T ) | S t = s , Y t = y .
Note that V ( T ) depends on the type of options. In this paper, we consider two types of path-dependent options: down-and-out put options and floating strike lookback put options. For notational convenience, we put U t : = min 0 u t S u and Z t : = max 0 u t S u . The payoff of a down-and-out put option is given by
D O P ( T ) : = max { K S T , 0 } × 1 U T > B ,
where K is the strike price, B is the barrier level satisfying 0 < B < K , and 1 U T > B is the indicator function. For a floating strike lookback put option, its payoff has the form of L P f l o a t ( T ) : = Z T S T . Applying Itô’s lemma, we can obtain a partial differential equation (PDE) for P ( t , s , y ) as follows:
0 = P t + 1 2 s 2 f 2 ( y ) 2 P s 2 + r s P s P + 2 ρ ν s ϵ f ( y ) 2 P s y + ν 2 ϵ 2 P y 2 + 1 ϵ ( m y ) 2 ν ϵ Λ ( y ) P y .
The boundary conditions for Equation (3) vary depending on the type of options. For example, the boundary conditions for Equation (3) when V ( T ) = D O P ( T ) are
P ( T , s , y ) = max { K s , 0 } , s > B , P ( t , B , y ) = 0 , 0 t T .
When V ( T ) = L P f l o a t ( T ) , the boundary conditions become the following:
P z ( t , z , y , z ) = 0 , 0 t T , z > 0 , P ( T , s , y , z ) = z s , 0 s z .
Note that in this case, P is a function of t, s, y, and z (here, Z t = z ).
Remark: Since the Mellin transform of the payoff function of a call option is not defined, this paper primarily concentrates on evaluating put options. However, as outlined in Buchen (2001), the pricing of call options can be directly derived from put options through the put–call parity relationship.

3. Asymptotic Expansions

In this section, we apply an asymptotic expansion approach to establish partial differential equations, which will be used to derive an approximate solution to Equation (3) and thus find an approximated value of a put option.
We begin with re-organizing Equation (3) in terms of the orders of ϵ as follows:
1 ϵ L 0 P + 1 ϵ L 1 P + L 2 P = 0 ,
where the operators L 0 , L 1 and L 2 are defined by
L 0 : = m y y + ν 2 2 y 2 , L 1 : = 2 ρ ν s f y 2 s y 2 ν Λ y y ,   and L 2 : = t + 1 2 s 2 f 2 y 2 s 2 + r s s · .
In order to obtain an efficient approximate solution to P, as that in Fouque et al. (2011), we apply the following asymptotic expansion of P:
P = P 0 + ϵ P 1 + ϵ P 2 + ϵ ϵ P 3 + ,
where P 0 , P 1 , … are functions corresponding to varying orders of ϵ . Substituting P in Equation (5) into Equation (4) and re-organizing the terms, we obtain
0 = 1 ϵ L 0 P 0 + 1 ϵ L 1 P 0 + L 0 P 1 + L 0 P 2 + L 1 P 1 + L 2 P 0 + ϵ L 0 P 3 + L 1 P 2 + L 2 P 1 + .
Our aim is to find P 0 and P 1 .
Firstly, from the O ( 1 / ϵ ) -order term in Equation (6), we obtain L 0 P 0 = 0 . If we assume that P 0 does not grow as fast as e y 2 / 2 , as was assumed in Choi et al. (2013), we can show that P 0 is independent of y. Secondly, from the O ( 1 / ϵ ) -order term in Equation (6), we obtain L 1 P 0 + L 0 P 1 = 0 . Since P 0 is independent of y, then L 1 P 0 = 0 . It follows that L 0 P 1 = 0 . Again, if we assume that P 1 does not grow as fast as e y 2 / 2 , then we can deduce that P 1 is also independent of y.
Next, from the O ( 1 ) -order term in Equation (6), we obtain
L 0 P 2 + L 1 P 1 + L 2 P 0 = 0 .
Since P 1 is independent of y, we have L 1 P 1 = 0 , which implies that
L 0 P 2 + L 2 P 0 = 0 .
Seeing Equation (7) as a Poisson equation for P 2 in y, in order for it to have a solution, it is required to satisfy the centering condition
L 2 P 0 = L 2 P 0 = 0 ,
which is equivalent to
P 0 t + r s P 0 s + 1 2 s 2 f 2 2 P 0 s 2 r P 0 = 0 .
This is an equation for us to determine the P 0 term. Here, · denotes the expectation with respect to the invariant distribution of the process { Y t : t 0 } , i.e.,
h = + h ( y ) Φ ( y ) d y , where Φ ( y ) = 1 2 π ν 2 e ( y m ) 2 2 ν 2 .
Note that a small ϵ value corresponds to fast-mean reverting. In this case, Y t approaches a constant and f 2 can be regarded as constant variance, and then Equation (9) is the Black–Scholes PDE. Thus, for small ϵ , P 0 represents the put option price under the Black–Scholes model.
Following Equation (8), we have
L 2 P 0 = L 2 P 0 L 2 P 0 = 1 2 f 2 f 2 s 2 2 P 0 s 2 ,
which, together with Equation (7), implies
L 0 P 2 = 1 2 f 2 f 2 s 2 2 P 0 s 2 .
The solution to Equation (10) can be expressed as
P 2 = 1 2 ϕ + c s 2 2 P 0 s 2 ,
where ϕ is a function of y which only satisfies the equation L 0 ϕ = f 2 f 2 , and c is a function of other variables except y.
To derive an equation for P 1 , we consider the O ( ϵ ) -term in Equation (6) and obtain
L 0 P 3 + L 1 P 2 + L 2 P 1 = 0 .
This equation can be regarded as a Poisson equation for P 3 in y, and in order for it to have a solution, the following centering condition must be satisfied:
L 1 P 2 + L 2 P 1 = 0 .
After we substitute P 2 in Equation (11) into Equation (12) and make simplifications, we obtain
P 1 t + 1 2 f 2 s 2 2 P 1 s 2 + r s P 1 s r P 1 = c 1 s 3 3 P 0 s 3 + c 2 s 2 2 P 0 s 2 ,
where
c 1 : = 2 2 f ϕ ρ ν and c 2 : = 2 2 2 ρ f ϕ Λ ϕ ν .
This is an equation for us to determine the first correction term P 1 .
We summarize the previous formal analysis as the following theorem.
Theorem 1.
Under the SV model governed by Equation (1), an approximation of the risk-neutral value P of a path-dependent put option is given by
P = P 0 + ϵ P 1 + o ( ε ) ,
for small ϵ, where P 0 and P 1 are determined by Equations (9) and (13) with corresponding boundary conditions, respectively, such that P 0 is the put option price under the Black–Scholes model with constant effective volatility f 2 and P 1 is the first-order correction term.
Finally, as mentioned in Section 2, boundary conditions for Equations (8) and (13) depend on the types of options that we consider. We describe the corresponding boundary conditions and solve these equations in the next two sections.

4. Determining P 0 and P 1 for Down-and-Out Put Options

In this section, we use Mellin transform to derive analytical expressions of the P 0 and P 1 terms for down-and-out put options

4.1. P 0 Term for Down-and-Out Put Options

In order to use Mellin transform to calculate the P 0 term for down-and-out put options, noting that P 0 is independent of y under our assumption, we first follow the method in Buchen (2001) and use the boundary condition,
P ( T , s , y ) = max { K s , 0 } , for s > B ,
to set up the boundary condition of P 0 for s 0 as follows:
P 0 T , s : = K s 1 B < s < K B s k 1 1 K B 2 s 1 B 2 K < s < B ,
where k 1 = 2 r / f 2 . Now, we apply Mellin transform to Equation (9) to convert this PDE into the following ODE:
d P ^ 0 d t + 1 2 f 2 ( w 2 + w ) r w r P ^ 0 = 0 .
The solution to Equation (17) is given by
P ^ 0 t , w = θ ^ ( w ) e 1 2 f 2 w 2 + 1 k 1 w k 1 T t ,
where θ ^ is a function of w, determined by the boundary condition (16).
Next, we take inverse Mellin transform of Equation (18) and obtain
P 0 ( t , s ) = P 0 ( T , s ) M 1 e λ w + η 2 + δ ,
where
λ = 1 2 f 2 T t , η = 1 k 1 2 , δ = λ η 2 r T t
and the operation * means the convolution. Applying Table A1 in Appendix A and the boundary condition given in Equation (16), we have
P 0 t , s = P 0 T , s e δ s η 2 λ π e 1 4 λ ln s 2 = B K K u e δ s u η 1 2 λ π e 1 4 λ ln s u 2 d u u B 2 K B B u k 1 1 K B 2 u e δ s u η 1 2 λ π e 1 4 λ ln s u 2 d u u .
After some careful calculation, for down-and-out put options, we derive a closed-form expression of the P 0 term as follows:
P 0 ( t , s ) = K e r T t Φ Δ s K Φ Δ s B s Φ Δ + s K Φ Δ + s B K e r T t B s k 1 1 Φ Δ B s Φ Δ B 2 s K + B B s k 1 Φ Δ + B s Φ Δ + B 2 s K ,
where Φ ( · ) is the CDF of the standard normal distribution and
Δ ± ( x ) = 1 f 2 ( T t ) ln ( x ) + r ± 1 2 f 2 ( T t ) .
Note that P 0 given in Equation (20) is precisely the same as the price of a down-and-out put option given in the literature, e.g., Hull (2015, chp. 26, p. 606) or Haug (2006, chp. 4), if we let σ 2 = f 2 . For details of the derivation of formula (20), we refer the reader to Appendix B.

4.2. P 1 Term for Down-and-Out Put Options

For down-and-out put options, the boundary conditions for P 1 are
P 1 ( T , s ) = 0 , f o r s B , P 1 ( t , B ) = 0 , f o r 0 < t < T .
We again follow the method in Buchen (2001) and extend the boundary conditions P 1 ( T , s ) = 0 , for s B as P 1 T , s = 0 for all s 0 .
Next, we apply Mellin transform to Equation (13) to obtain
d P ^ 1 d t + 1 2 f 2 w 2 + w r w r P ^ 1 = c 1 w w + 1 w + 2 + c 2 w w + 1 P ^ 0 .
Solving this equation, we obtain
P ^ 1 t , w = c 1 T t w 3 c 2 3 c 1 T t w 2 c 2 2 c 1 T t w P ^ 0 t , w .
Finally, applying inverse Mellin transform, we obtain an explicit closed-form expression of P 1 as follows:
P 1 t , s = M 1 P ^ 1 t , w = c 1 T t s d d s P 0 t , s 3 s 2 d 2 d s 2 P 0 t , s s 3 d 3 d s 3 P 0 t , s c 2 3 c 1 T t s d d s P 0 t , s + s 2 d 2 d s 2 P 0 t , s c 2 2 c 1 T t s d d s P 0 t , s ,
where P 0 is given in the previous section and c 1 and c 2 are given in Equation (14).
We summarize the above analysis and calculation on down-and-out put options in the following theorem.
Theorem 2.
Under the SV model governed by Equation (1), an approximation of the risk-neutral value P of a down-and-out barrier put option is given by
P = P 0 + ϵ P 1 + o ( ε ) ,
where P 0 and P 1 are given by Equations (20) and (21), respectively.

5. Determining P 0 and P 1 for Lookback Put Options

In this section, we use Mellin transform to derive analytical expressions of the P 0 and P 1 terms for floating strike lookback put options

5.1. P 0 Term for Lookback Put Options

For lookback floating strike put options, the boundary conditions of P 0 are
P 0 z ( t , z , z ) = 0 , P 0 z T , s , z = 1 , f o r 0 < s < z .
Similar to the case of down-and-out put options, we extend the second boundary condition to 0 < s < as follows:
P 0 z T , s , z : = 1 s < z z s k 1 1 · 1 z < s , for 0 < s < .
Then, by integrating each side of the last equation, we can obtain
P 0 T , s , z = s z ξ s k 1 1 d ξ = 1 k 1 z s k 1 s + 1 k 1 s
for s > z . For convenience, we let u = s / z and Q 0 = P 0 / z . With these notations, Equation (9) becomes
Q 0 t + 1 2 u 2 f 2 2 Q 0 u 2 + r u Q 0 u r Q 0 = 0 ,
with boundary conditions
Q 0 T , u = 1 k 1 u 1 k 1 + 1 k 1 u , for u > 1 ,
and Q 0 ( T , u ) = 1 , for 0 < u < 1 .
Note that except the boundary conditions, Equation (24) is identical to Equation (9). Applying Mellin transform in the same way as that for the case of down-and-out put options, we can derive the solution to Equation (24) as follows:
Q 0 ( t , u ) = θ ^ ( w ) M 1 e λ w + η 2 + δ .
Again, applying Table A1 and P 0 given in Equation (16), we have
Q 0 ( t , u ) = Q 0 ( T , u ) e δ z η 1 2 π λ 1 2 e 1 4 λ ln z 2 = 0 1 1 ξ e δ u ξ η 1 2 π λ 1 2 e 1 4 λ ln u ξ 2 d ξ ξ + 1 1 k 1 ξ 1 k 1 + ξ k 1 e δ u ξ η 1 2 π λ 1 2 e 1 4 λ ln u ξ 2 d ξ ξ .
After calculating integrals, for floating strike lookback put options, we derive a closed-form expression of the P 0 term as follows:
P 0 ( t , s , z ) = z e r T t Φ Δ s z s Φ Δ + s z z k 1 s z 1 k 1 e r T t Φ Δ z s + s k 1 Φ Δ + s z ,
where Φ ( · ) is the CDF of the standard normal distribution. Note that P 0 given in Equation (27) is precisely the same as the price of a floating strike put option given in the literature, e.g., Hull (2015, chp. 26, p. 608) or Haug (2006, chp. 4), if we let σ 2 : = f 2 . Details of the derivation of this formula can be found in Appendix B.

5.2. P 1 Term for Lookback Put Options

For floating strike lookback put options, the boundary conditions for P 1 are
P 1 ( T , s , z ) = 0 , f o r 0 < s < z , P 1 z t , z , z = 0 , f o r 0 < t < T a n d z > 0 .
Just like that for the P 0 -term for floating strike lookback put options, we let u = s / z and Q 1 = P 1 / z . With these notation changes, Equation (13) is converted to the following:
Q 1 t + 1 2 f 2 u 2 2 Q 1 u 2 + r u Q 1 u r Q 1 = c 1 u 3 3 Q 0 u 3 + c 2 u 2 2 Q 0 u 2
with Q 1 ( T , u ) = 0 for 0 < u < 1 .
Note that Equation (28) is essentially the same as Equation (13), except the notational difference. So, we have
Q 1 t , u = c 1 T t u d d u Q 0 t , u 3 u 2 d 2 d u 2 Q 0 t , u u 3 d 3 d u 3 Q 0 t , u c 2 3 c 1 T t u d d u Q 0 t , u + u 2 d 2 d z 2 Q 0 t , u c 2 2 c 1 T t u d d u Q 0 t , u ,
where Q 0 is given previously. Consequently, we have
P 1 t , s , z = c 1 T t s d d s P 0 t , s , z 3 s 2 d 2 d s 2 P 0 t , s , z s 3 d 3 d s 3 P 0 t , s , z c 2 3 c 1 T t s d d s P 0 t , s , z + s 2 d 2 d s 2 P 0 t , s , z c 2 2 c 1 T t s d d s P 0 t , s , z ,
where c 1 and c 2 are the same as those defined previously.
We summarize the above analysis and calculation on floating strike lookback put options in the following theorem.
Theorem 3.
Under the SV model governed by Equation (1), an approximation of the risk-neutral value P of a floating strike lookback put option is given by
P = P 0 + ϵ P 1 + o ( ε ) ,
where P 0 and P 1 are given by Equations (27) and (30), respectively.

6. Numerical Results and Sensitivity Analysis

In this section, we conduct a numerical study to investigate the sensitivity of the first-order correction term P 1 and our approximation results P 0 + ϵ P 1 with respect to the initial value of underlying asset. This means that we set t = 0 throughout this section. We also compare the results given by our closed form formulas with those generated by the Monte-Carlo simulation.
First of all, as conducted by Fouque et al. (2000, 2011) and Cao et al. (2021), we choose f to take the following form:
f ( y ) = 0.35 tan 1 ( y ) + π 2 / π + 0.05 .
Secondly, the values of other parameters used in this section are given in Table 1 whenever they are required to be fixed.
Here, we do not choose precise values of β and ρ , and particular forms of ξ ( y ) (in Section 2) and ϕ ( y ) (in Section 3) to calculate the above values of c 1 and c 2 . Instead, c 1 and c 2 are calibrated from the term structure of the implied volatility surface as described in the book of Fouque et al. (2000). Specifically, the implied volatility I ϵ of a European vallina call option with fast mean-reverting stochastic process can be approximated by the following formula:
I ϵ = a ln ( K s ) T t + b + o ( ϵ )
with
a = c 1 f 2 3 / 2 and b = f 2 + c 1 f 2 3 / 2 r + 3 2 f 2 c 2 f 2 .
The parameters a and b are estimated as the slope and intercept of the regression fit of the observed implied volatilities as a linear function of logmoneyness-to-maturity-ratio ln ( K / s ) / ( T t ) . From the calibrated values a and b on the observed implied volatility surface, the parameters c 1 and c 2 are obtained as
c 1 = a σ f 2 3 / 2 and c 2 = f 2 ( ( f 2 b ) a ( r + 3 2 f 2 ) ) .
Thirdly, note that when t = 0 , s = z . Hence, in this case, the formula for P 0 given by Equation (27) is simplified.
Figure 1a shows how the ϵ P 1 -term for a down-and-out barrier put option changes with respect to a variation in ϵ values As we can see, for fixed ϵ , when s increases, P 1 decreases first, and then increases after it hits its trough. When ϵ becomes smaller (equivalently, the mean-reverting speed becomes larger), ϵ P 1 approaches to a zero. Figure 1b shows how the value of P 0 + ϵ P 1 for a down-and-out put option varies with respect to the change in ϵ values. As we can see, when the value of ϵ changes from 0.01 to 0.0001, the value of P 0 + ϵ P 1 does not vary much. In fact, the values of P 0 + ϵ P 1 match well with the result of Monte-Carlo simulation in all cases. Furthermore, in all cases, the value of P 0 + ϵ P 1 declines as s increases.
Figure 2a shows how the ϵ P 1 -term for a floating strike lookback put changes with respect to a variation in ϵ values. In a similar pattern, for a fixed ϵ -value, when s increases, P 1 decreases first and then increases after it hits its trough. Similar to the case of down-and-out put options, when ϵ becomes smaller (equivalently, the mean-reverting speed becomes larger), ϵ P 1 approaches to zero. Figure 2b shows how the value of P 0 + ϵ P 1 for a floating strike put varies with respect to the change in ϵ values. When the value of ϵ changes from 0.01 to 0.001, the value of P 0 + ϵ P 1 varies. But when the value of ϵ changes from 0.001 to 0.0001, the value of P 0 + ϵ P 1 does not vary much. The values of P 0 + ϵ P 1 match well with the result of Monte-Carlo simulation when ϵ = 0.001 or 0.0001 . Furthermore, in all cases, the value of P 0 + ϵ P 1 increases as s increases.
Figure 3a illustrates the variation in the value of P 0 + ϵ P 1 for a down-and-out put option in response to changes in ρ values. As ρ shifts from −0.6 to −0.4, there is a slight decrease in the value of P 0 + ϵ P 1 . Additionally, in all scenarios, the value of P 0 + ϵ P 1 shows an upward trend as s increases.
Figure 3b depicts the change in the value of P 0 + ϵ P 1 for a floating strike put concerning variations in ρ values. Similar to the previous case, a shift in ρ from −0.6 to −0.4 results in a minor decline in the value of P 0 + ϵ P 1 . Moreover, across all instances, an increase in s is associated with a rise in the value of P 0 + ϵ P 1 .

7. Concluding Remarks

This article establishes explicit closed-form solutions for first order approximations of down-and-out barrier and floating strike lookback put option prices under a stochastic volatility model by means of Mellin transform. The zero-order terms in the solutions for the prices of both types of put options coincide with those in Hull (2015) or Haug (2006) under the classical Back–Scholes model. Our numerical analysis shows that the results given by those explicit closed-form solutions match well with those generated by the Monte-Carlo simulation. This confirms the accuracy of the approximation. Furthermore, we also discussed the sensitivity of the first-order error terms and the approximation with respect to the underlying asset price and the mean-reverting speed of the OU-process which governs the volatility.
This model formula can be employed by financial professionals for the swift and precise pricing of barrier and lookback options. This is demonstrated by the efficiency of our formula in comparison to the conventional Monte-Carlo method. Our pricing formula offers an effective means of assessing barrier and lookback options. Looking ahead, we may extend our methodology to evaluate other path-dependent options in future works, including, but not limited to, Asian options, Russian options, and more.

Author Contributions

Conceptualization, J.C.; methodology, J.C. and W.Z.; software, X.L.; validation, J.C. and W.Z.; formal analysis, X.L.; investigation, X.L.; resources, X.L.; data curation, X.L.; writing—original draft preparation, X.L.; writing—review and editing, J.C. and W.Z.; visualization, X.L.; supervision, J.C. and W.Z.; project administration, W.Z.; funding acquisition, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank Jeong-Hoon Kim for valuable discussions and suggestions on Mellin transform method. The authors also thank the anonymous reviewers for their careful reading of our manuscript and their many insightful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Mellin Transform

The Mellin transform is an integral transform that may be regarded as the multiplicative version of the two-sided Laplace transform. It is often used in the theory of asymptotic expansions. For a locally Lebesgue integrable function h : R + R , the Mellin transform denoted by M h or h ^ is given by
h ^ ( w ) = M h ( w ) : = 0 + s w 1 h ( s ) d s , w C ,
and if a < Re ( w ) < b and c such that a < c < b exists, the inverse of the Mellin transform is expressed by
h ( s ) = M 1 h ^ ( s ) = 1 2 π i c i c + i s w h ^ ( w ) d w .
In this paper, we use the following properties of Mellin transform.
Table A1. List of properties of Mellin transform used in this paper.
Table A1. List of properties of Mellin transform used in this paper.
FunctionMellin Tansform
h h ^
s h w h ^
s 2 h w ( w + 1 ) h ^
s 3 h ( 3 ) w ( w + 1 ) ( w + 2 ) h ^
e δ s η 2 λ π e 1 4 λ ( ln s ) 2 e λ ( w + η ) 2 + δ
s h + s 2 h w 2 h ^
s h 3 s 2 h s 3 h ( 3 ) w 3 h ^
Here, λ , η , and δ are not related to w or s, and h , h , and h ( 3 ) are the first-order, second-order, and third-order derivatives of h, respectively.

Appendix B. Derivation of Formulas (20) and (27)

In this appendix, we give detailed derivation of the Formulas (20) and (27).

Appendix B.1. Derivation of Formula (20)

From Equation (19), we know that
P 0 t , s = B K K u e δ s u η 1 2 λ π e 1 4 λ ln s u 2 d u u B 2 K B B u k 1 1 K B 2 u e δ s u η 1 2 λ π e 1 4 λ ln s u 2 d u u .
By letting v = ln u , we convert the first integral to
ln B ln K ( K e v ) s η e δ e η v 1 2 λ π e 1 4 λ ( ln s v ) 2 d v = s η e δ 2 λ π ln B ln K K e 1 4 λ ( v 2 2 v ln s + ( ln s ) 2 + 4 λ η v ) d v ln B ln K e 1 4 λ ( v 2 2 v ln s + ( ln s ) 2 + 4 λ ( η 1 ) v ) d v = s η e δ 2 λ π ln B ln K K e 1 4 λ ( v ln s + 2 λ η ) 2 + λ η 2 η ln s d v ln B ln K e 1 4 λ [ v ln s + 2 λ ( η 1 ) ] 2 + λ ( η 1 ) 2 ( η 1 ) ln s d v ,
we further apply the following changes in variables:
x : = v ln s + 2 λ η 2 λ and x : = v ln s + 2 λ ( η 1 ) 2 λ
to obtain
ln B ln K ( K e v ) s η e δ e η v 1 2 λ π e 1 4 λ ( ln s v ) 2 d v = e δ 2 π K e λ η 2 ln ( B s ) + 2 λ η 2 λ ln ( K s ) + 2 λ η 2 λ e x 2 2 d x s e λ ( η 1 ) 2 ln ( B s ) + 2 λ ( η 1 ) 2 λ ln ( K s ) + 2 λ ( η 1 ) 2 λ e x 2 2 d x = K e δ + λ η 2 Φ ln ( K s ) + 2 λ η 2 λ Φ ln ( B s ) + 2 λ η 2 λ s e δ + λ ( η 1 ) 2 Φ ln ( K s ) + 2 λ ( η 1 ) 2 λ Φ ln ( B s ) + 2 λ ( η 1 ) 2 λ .
Now, if we plug into δ , η , and λ into the above formula, we derive
ln B ln K ( K e v ) s η e δ e η v 1 2 λ π e 1 4 λ ( ln s v ) 2 d v = K e r T t Φ Δ s K Φ Δ s B s Φ Δ + s K Φ Δ + s B .
Similarly, we can evaluate the second integral
B 2 K B B u k 1 1 K B 2 u e δ s u η 1 2 λ π e 1 4 λ ln s u 2 d u u
to obtain
K e r T t B s k 1 1 Φ Δ B s Φ Δ B 2 s K B B s k 1 Φ Δ + B s Φ Δ + B 2 s K .
Putting these two integrals together yields Formula (20).

Appendix B.2. Derivation of Formulas (27)

From Equation (26), we have
Q 0 ( t , u ) = 0 1 1 ξ e δ u ξ η 1 2 λ π e 1 4 λ ln u ξ 2 d ξ ξ + 1 1 k 1 ξ 1 k 1 + ξ k 1 e δ u ξ η 1 2 λ π e 1 4 λ ln u ξ 2 d ξ ξ .
We let v = ln ξ . For the first integral, we have
0 1 ( 1 ξ ) e δ u ξ η 1 2 λ π e 1 4 λ ln u ξ 2 d u u = 0 u η 1 e v e δ v η 1 2 λ π e 1 4 λ ( ln u v ) 2 d v = u η e δ 2 λ π 0 e 1 4 λ v 2 2 v ln u + ln u 2 + 4 λ η v d v 0 e 1 4 λ v 2 2 v ln u + ln u 2 + 4 λ η 1 v d v = u η e δ 2 λ π 0 e 1 4 λ v ln u + 2 λ η 2 + λ η 2 η ln u d v 0 e 1 4 λ v ln u + 2 λ η 1 2 + λ η 1 2 η 1 ln u d v .
Next, we let
v : = v ln u + 2 λ η 2 λ and v : = v ln u + 2 λ η 1 2 λ .
Then, we have
0 1 ( 1 ξ ) e δ u ξ η 1 2 λ π e 1 4 λ ln u ξ 2 d u u = e δ 2 π ln u + 2 λ η 2 λ e v 2 2 + λ η 2 d v u ln u + 2 λ ( η 1 ) 2 λ e v 2 2 + λ η 1 2 d v = e δ + λ η 2 Φ ln u + 2 λ η 2 λ u e δ + λ η 1 2 Φ ln u + 2 λ η 1 2 λ = e r T t Φ Δ s z s z Φ Δ + s z .
For the second integral, we have
1 1 k 1 ξ 1 k 1 + ξ k 1 e δ u ξ η 1 2 λ π e 1 4 λ ln u ξ 2 d ξ ξ = 0 1 k 1 e ( 1 k 1 ) v + 1 k 1 e v e δ u η e v η 1 2 λ π e 1 4 λ ln u v 2 d v = e δ u η 2 k 1 λ π 0 e η v 1 4 λ ln u v 2 + e v 1 η 1 4 λ ln u v 2 d v = e δ u η 2 k 1 λ π 0 e 1 4 λ v ln u 2 λ η 2 + λ η 2 + η ln u d v + 0 e 1 4 λ v ln u 2 λ ( 1 η ) 2 + λ ( 1 η ) 2 + ( 1 η ) ln u d v ,
where we use the fact that k 1 1 + η = η . Furthermore, we introduce a new variable
v : = v ln u 2 λ η 2 λ .
Then, we have
1 1 k 1 ξ 1 k 1 + ξ k 1 e δ u ξ η 1 2 λ π e 1 4 λ ln u ξ 2 d ξ ξ = e δ u η k 1 2 π ln u 2 λ η 2 λ e v 2 2 e λ η 2 + η ln u d v + ln u + 2 λ η 1 2 λ e v 2 2 e λ η 1 2 + 1 η ln u d v = 1 k 1 e δ + λ η 2 u 1 k 1 Φ ln u + 2 λ η 2 λ + 1 k 1 u e δ + λ η 1 2 Φ ln u + 2 λ 1 η 2 λ = 1 k 1 s z 1 k 1 e r T t Φ Δ z s + 1 k 1 s z Φ Δ + s z .
Putting these two integrals together and using the fact that P 0 = z Q 0 , we can obtain our Formula (27).

References

  1. Bates, David S. 1996. Jumps and stochastic volatility: Exchange rate processes implicit in Deutsche market options. The Review of Financial Studies 9: 69–107. [Google Scholar] [CrossRef]
  2. Black, Fischer, and Myron Scholes. 1973. The pricing of options and corporate liabilities. Journal of Political Economy 3: 637–54. [Google Scholar] [CrossRef]
  3. Boyle, Phelim P., and Yisong S. Tian. 1999. Pricing lookback and barrier options under the CEV process. Journal of Financial and Quantitative Analysis 34: 241–64. [Google Scholar] [CrossRef]
  4. Buchen, Peter. 2001. Image options and the road to barriers. Risk Magazine 14: 127–30. [Google Scholar]
  5. Cao, Jiling, Jeong-Hoon Kim, and Wenjun Zhang. 2021. Pricing variance swaps under hybrid CEV and stochastic volatility. Journal of Computational and Applied Mathematics 386: 113220. [Google Scholar] [CrossRef]
  6. Cao, Jiling, Jeong-Hoon Kim, and Wenjun Zhang. 2023. Valuation of barrier and lookback options under hybrid CEV and stochastic volatility. Mathematics and Computers in Simulation 208: 660–76. [Google Scholar] [CrossRef]
  7. Chiarella, Carl, Boda Kang, and Gunter H. Meyer. 2012. The evaluation of barrier option prices under stochastic volatility. Computers & Mathematics with Applications 64: 2034–48. [Google Scholar]
  8. Choi, Sun-Yong, Jean-Pierre Fouque, and Jeong-Hoon Kim. 2013. Option pricing under hybrid stochastic and local volatility. Quantitative Finance 13: 1157–65. [Google Scholar] [CrossRef]
  9. Clewlow, Les, Javier Llanos, and Chris Strickland. 1994. Pricing Exotic Options in a Black-Scholes World. Coventry: Financial Operations Research Centre, University of Warwick, vol. 54, 32p. [Google Scholar]
  10. Conze, Antoine, and R. Vishwanathan. 1991. Path-dependent options: The case of lookback options. The Journal of Finance 46: 1893–907. [Google Scholar] [CrossRef]
  11. Cox, John. 1975. Notes on Option Pricing I: Constant Elasticity of Variance Diffusions. Working Paper. Stanford: Stanford University. [Google Scholar]
  12. Cox, John. 1996. The constant elasticity of variance option pricing model. Journal of Portfolio Management 5: 15–17. [Google Scholar] [CrossRef]
  13. Davydov, Dmitry, and Vadim Linetsky. 2001. Pricing and hedging path-dependent options under the CEV process. Management Science 47: 881–1027. [Google Scholar] [CrossRef]
  14. Fouque, Jean-Pierre, George Papanicolaou, and Ronnie Sircar. 2000. Derivatives in Financial Markets with Stochastic Volatility. Cambridge: Cambridge University Press. [Google Scholar]
  15. Fouque, Jean-Pierre, George Papanicolaou, Ronnie Sircar, and Knut Sølna. 2011. Multiscale Stochastic Volatility for Equity, Interest Rate, and Credit Derivatives. Cambridge: Cambridge University Press. [Google Scholar]
  16. Funahashi, Hideharu, and Tomohide Higuchi. 2018. An analytical approximation for single barrier options under stochastic volatility models. Annals of Operations Research 266: 129–57. [Google Scholar] [CrossRef]
  17. Goldman, M. Barry, Howard B. Sosin, and Marry Ann Gatto. 1979. Path dependent options: “Buy at the low, sell at the high”. The Journal of Finance 34: 1111–27. [Google Scholar]
  18. Haug, Espen. 2006. The Complete Guide to Option Pricing Formulas, 2nd ed. New York: McGraw-Hill. [Google Scholar]
  19. Heston, Steven. 1993. A closed-form solution for options with stochastic volatility with applications to bond and currency options. The Review of Financial Studies 6: 327–43. [Google Scholar] [CrossRef]
  20. Hull, John C. 2015. Options, Futures and Other Derivatives, 9th ed. London: Pearson. [Google Scholar]
  21. Kato, Takashi, Akihiko Takahashi, and Toshihiro Yamada. 2013. An asymptotic expansion formula for up-and-out barrier option price under stochastic volatility model. Japan Society for Industrial and Applied Mathematics Letters 5: 17–20. [Google Scholar] [CrossRef]
  22. Kim, Hyun-Gyoon, Jiling Cao, Jeong-Hoon Kim, and Wenjun Zhang. 2023. A Mellin transform approach to pricing barrier options under stochastic elasticity of variance. Applied Stochastic Models in Business and Industry 39: 160–76. [Google Scholar] [CrossRef]
  23. Kim, So-Yeun, and Ji-Hun Yoon. 2018. An approximated European option price under stochastic elasticity of variance using Mellin transforms. East Asian Mathematical Journal 34: 239–48. [Google Scholar]
  24. Leung, Kwai Sun. 2013. An analytic pricing formula for lookback options under stochastic volatility. Applied Mathematics Letters 26: 145–49. [Google Scholar] [CrossRef]
  25. Merton, Robert C. 1973. The theory of rational option pricing. Bell Journal of Economics and Management Science 1: 141–83. [Google Scholar] [CrossRef]
  26. Panini, Radha, and Ram P. Srivastav. 2004. Option pricing with Mellin transforms. Mathematical and Computer Modelling 40: 43–56. [Google Scholar] [CrossRef]
  27. Park, Sang-Hyeon, and Jeong-Hoon Kim. 2013. A semi-analytic pricing formula for lookback options under a general stochastic volatility model. Statistics & Probability Letters 83: 2537–43. [Google Scholar]
  28. Reiner, Eric, and M. Rubinstein. 1991. Breaking down the barriers. Risk 4: 28–35. [Google Scholar]
  29. Wirtu, Teferi Dereje, Philip Ngare, and Ananda Kube. 2017. Pricing floating strike lookback put option under Heston stochastic volatility. Global Journal of Mathematical Sciences: Theory and Practical 9: 427–39. [Google Scholar]
  30. Yoon, Ji-Hun. 2014. Mellin transform method for European option pricing with Hull-White stochastic interest rate. Journal of Applied Mathematics 2014: 759562. [Google Scholar] [CrossRef]
  31. Zaevski, Tsvetelin S., Young Shin Kim, and Frank J. Fabozzi. 2014. Option pricing under stochastic volatility and tempered stable Lévy jumps. International Review of Financial Analysis 31: 101–8. [Google Scholar] [CrossRef]
Figure 1. Plots of ϵ P 1 and P 0 + ϵ P 1 with different values of ϵ against the initial value of the underlying asset, for the down-and-out put option.
Figure 1. Plots of ϵ P 1 and P 0 + ϵ P 1 with different values of ϵ against the initial value of the underlying asset, for the down-and-out put option.
Jrfm 16 00456 g001
Figure 2. Plots of ϵ P 1 and P 0 + ϵ P 1 with different values of ϵ , against the initial value of the underlying asset, for floating strike put options.
Figure 2. Plots of ϵ P 1 and P 0 + ϵ P 1 with different values of ϵ , against the initial value of the underlying asset, for floating strike put options.
Jrfm 16 00456 g002
Figure 3. Plots P 0 + ϵ P 1 with different values of ρ , against the initial value of the underlying asset, for down-and-out put option and floating strike put option.
Figure 3. Plots P 0 + ϵ P 1 with different values of ρ , against the initial value of the underlying asset, for down-and-out put option and floating strike put option.
Jrfm 16 00456 g003
Table 1. The role and numerical value of parameters.
Table 1. The role and numerical value of parameters.
ParameterRoleValue
rrisk-free interest rate0.035
Bbarrier level1500
Kput option strike price2700
Tmaturity time1
c 1 as defined in Section 3−0.004
c 2 as defined in Section 3−0.018
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cao, J.; Li, X.; Zhang, W. Pricing Path-Dependent Options under Stochastic Volatility via Mellin Transform. J. Risk Financial Manag. 2023, 16, 456. https://doi.org/10.3390/jrfm16100456

AMA Style

Cao J, Li X, Zhang W. Pricing Path-Dependent Options under Stochastic Volatility via Mellin Transform. Journal of Risk and Financial Management. 2023; 16(10):456. https://doi.org/10.3390/jrfm16100456

Chicago/Turabian Style

Cao, Jiling, Xi Li, and Wenjun Zhang. 2023. "Pricing Path-Dependent Options under Stochastic Volatility via Mellin Transform" Journal of Risk and Financial Management 16, no. 10: 456. https://doi.org/10.3390/jrfm16100456

APA Style

Cao, J., Li, X., & Zhang, W. (2023). Pricing Path-Dependent Options under Stochastic Volatility via Mellin Transform. Journal of Risk and Financial Management, 16(10), 456. https://doi.org/10.3390/jrfm16100456

Article Metrics

Back to TopTop