The Pricing of Dual-Expiry Exotics with Mean Reversion and Jumps

This paper develops a new class of models for pricing dual-expiry options that are characterized by two expiry dates. The underlying asset price is modeled by a time changed exponential Ornstein Uhlenbeck (OU) process, where the time change process is a Lévy subordinator. The new models can capture both mean reversion and jumps often observed in various types of underlying assets of exotics. The pricing method exploits the observation that dual expiry options have payoffs that can be perfectly replicated by a particular set of first and second order binary options. The novelty of the paper is that we are able to derive the analytical solutions to the prices of these binaries through eigenfunction expansion method. Based on that, we can obtain the formulas for dual-expiry exotics through static replication. We also numerically investigate the sensitivities of prices of chooser, compound and extendable options with respect to the parameters of the models.


Introduction
Dual-expiry option is a class of exotic options that depends on a single underlying asset but whose payoff structure involves two fixed future dates.Usually, at the first expiry date, the option holder receives a contract that matures at the second expiry date.Chooser options, compound options, extendable options, shout options, American call options on an asset with a single known dividend payment and partial barrier options are good examples.The dual-expiry exotics Journal of Mathematical Finance can be easily generalized to multiple-expiry options, whose payoff structure involves several fixed future dates.Examples include ladder options, multiple shout options, Bermudan options, multiple extendable options and barrier options with discrete monitoring.[1] develops a new technique for pricing dual-expiry exotic options by assuming the underlying asset price dynamics follows continuous geometric Brownian motion.The method exploits the observation that dual expiry options have payoffs that can be perfectly replicated by a particular set of first and second order binary options.Hence, in order to avoid arbitrage, the exotic option prices are obtained by static replication with respect to this family of binaries.As pointed out by [1], the representation of prices in terms of binaries is quite general and does not depend on any particular underlying asset price dynamics.In the Black-Scholes (BS) framework, [2] extends the method of [1] for pricing multiple-expiry options using a concept of higher order binary options.
[3] [4] generalize [1] for dual-and multiple-expiry exotics to a Lévy environment.The binary option prices are derived using the mathematical methodology of [5] that employs the pseudo-differential operators whose symbol is expressed in terms of the characteristic exponent of the underly Lévy process.[6] further generalizes method of [1] to exotics with multiple-asset and multiple-expiry.
In this paper, we present a new class of models for pricing dual-expiry exotic options.We assume the underlying asset price is an exponential function of OU process subject to stochastic time change.The time change process is modeled by a Lévy subordinator.The resulting asset price then exhibits mean reversion and jumps.
We generalize the works of [1] [2] [3] [4] in two directions.In the previous works, the mean reversion is absent in the underlying asset price dynamics.As a result, these models are restricted to the dual-expiry exotics written on equities.
On the other hand, the empirical evidence on mean reversion in financial asset is abundant.There is a census that many underlying assets of exotic option contracts, such as currencies, commodities, energy, temperature and even some stocks, display mean reversion.Our models are built on the OU process and can therefore be applied to a wider range of asset classes.Time-changing a continuous time Markov process such as OU process can lead to a much wider class of models than classical jump-diffusion models.In the typical jump-diffusion models, such as [3] [4], the state-independent jumps are added to the diffusion process.In these models, upon arrival, the direction of the jump and the probability distribution of jump amplitude are independent of the current state of the process.Through Lévy subordination of OU process, our model, in contrast, can feature state-dependent mean reverting jumps with the jump direction and the jump amplitude dependent on the current state of the process.The subordinate OU process is found to be a better candidate for modeling phenomena where jumps depend on the state, such as mean-reverting jumps.See its successful applications for callable and putable bonds in [7], Journal of Mathematical Finance commodities in [8], electricity in [9] and variance swaps in [10].We also refer to [11] for a brief introduction of subordinate Markov processes and [12] [13] [14] [15] for their further applications.
To solve the dual-expiry exotics pricing problems, we extend the works of [1] to the case where the asset price follows the mean reverting process with jumps.
We are able to derive the analytical pricing formulas for first and second order binary options using eigenfunction expansion method.After that, we apply them to pricing of some dual-expiry exotics such as chooser options, compound options and extendable options by static replication.We need to emphasize that eigenfunction expansion method is particularly suitable for pricing contingent claims written on the subordinate processes.The subordinate process is as analytically tractable as the original process without time change.The subordinate process shares the same eigenfunctions with the original one and the only modification is the replacement of eigenvalues of the original process with the Laplace transform of the eigenvalues.To the best of our knowledge, this is the first piece of work that applies this method for pricing dual-expiry options and it can be extended with little effort to the case of multiple-expiry exotics.We refer to [11] [16] for the surveys on the eigenfunction expansion method and [7] [8] [9] [10] [12] [13] [14] [15] for its various applications.
The rest of the paper is organized as follows.In Section 2, we introduce the general framework for modeling asset price as a time changed exponential OU process, where the time change process is modeled by the Lévy subordinators.In Section 3, we introduce the eigenfunction expansion method for our new model and also discuss how to calculate some important integrals that are essential for the determination of the eigenfunction expansion coefficients.In Section 4, we apply the eigenfunction expansion method to the valuation of first and second order binary options.We provide the analytical formulas for these binaries.In Section 5, we express the payoffs of some dual-expiry options in terms of binary options and then demonstrate how to valuate these options by static replication.
We also implement the model and analyze the effect of parameters of the model on the selected option prices through specific numerical examples.

Let ( )
, ,Q Ω  be a probability space with an information filtration ( t  ).
Suppose under the risk neutral measure Q, the asset price process S is governed by a time-changed exponential OU process, that is, where where T is a time change process and X is an OU process where 0 κ > and ( ) B t is a standard Brownian motion.
To introduce jumps into the asset price dynamics, we will follow [7]- [15] and model the time change process T as a Lévy subordinator.The Lévy subordinator T is a nondecreasing process with positive jumps and non-negative drift with the Laplace transform: where φ is the Lévy exponent and given by the Lévy-Khintchine formula (see e.g., [17]) where 0 γ ≥ and the Lévy measure ν must satisfy To make sure that the expectation in ( 4) is finite for λ ∈  , we also impose the following restriction We can also reparameterize the exponent by setting ( ) ( ) ( ) According to [11], when the process X is time changed by a Lévy subordinator T, the resulting process Y will be a jump-diffusion process with mean-reverting diffusion drift and mean-reverting jumps if 0 γ > or a pure jump process with mean-reverting jumps if 0 γ = .

Eigenfunction Expansion Method
For the OU process X in (3), its infinitesimal generator  is defined by where f is transformation function.f ′ and f ′′ are first-and second-order derivatives of f, respectively. Let denote the space of functions on ( ) and endowed with the inner product ( ) and Then, for any , we have (see e.g., [11]) where ( ) For the OU process X defined in (3), its eigenvalues and eigenfunctions can be summarized in the following result (see e.g., [11]): Proposition 1 For the OU process X defined in (3), the eigenvalues n λ and eigenfunction n ψ , and where ( ) For the time changed OU process Y defined in (2), for any , we can also employ the eigenfunction expansion method to compute the following expectation (see e.g., [8]): where { } n λ and { } n ψ are the eigenvalues and eigenfunctions of OU process and can be obtained from ( 10) and (11), respectively.φ is the Lévy exponent for the time change process T.
It is clear that the eigenfunction expansion of Y has the same form as X, but with ( ) . Thus, the eigenfunction expansion method makes the time changed model as tractable as the original model.This important result explains why eigenfunction expansion method is the natural way of computing option prices for the time-changed processes.
For our new pricing model, the eigenvalues and eigenfunctions of OU process can be calculated easily from Proposition 1.To employ the eigenfunction expansion method to calculate the dual-expiry exotics prices, we still need to obtain the eigenfunction expansion coefficient n f .In this section, we provide the formulas for several integrals that will later be employed to calculate n f .
The following formulas can be found in [18]: The following integrals can be computed using the results of [7]: , n c u α can be computed recursively as follows: , exp , .
is the CDF of standard normal distribution.

( )
, m n u γ can be calculated recursively as follows:

Valuation of First and Second Order Binaries
Following [1], we define the first order up (or down) binary option as the option that delivers an agreed payoff on expiry date if the price of underlying asset is above (or below) a fixed exercise price and zero otherwise.Let Y be the time changed OU process defined in (2) conditioning on ( ) Clearly, the up and down binaries satisfy the parity relation We also denote three specific options as follows.
In the following theorem, we can derive the analytical formulas for the above three option prices.
Proposition 2 Assume the processes for the underlying asset price are given by ( 1)-( 3).Let T be a Lévy subordinator with Lévy exponent φ .Then, 1) where n λ and ( ) n y ψ are in (10) and (11), respectively.Furthermore, and where n α and n α can be calculated using Lemma 1 and Lemma 2, respectively. 2) where and where n β and n β can be calculated using Lemma 1 and Lemma 2, respectively. 3)

Q y t t k s A y t t kB y t t
Proof.To prove 1), using eigenfunction expansion, we have ( where the equality in the last line comes from Lemma 2. For ( ) The proof for 2) is similar to 1).The only differences are that we need to compute which can both be analytically solved using Lemma 1 and Lemma 2.
3) can be obtained using the definition of Q-options.
We also define the second order binary option as an expiry 0 t binary contract on an underlying expiry 1 t binary option, 0 1 t t < .Let ( ) ( ) a first order binary, 0 1 , s s be up-down indicators at times 0 t and 1 t , respectively and 0 1 , ξ ξ be their corresponding exercise prices.The time t value of the second order binaries, conditioning on ( ) It is clear that the following up-down parity relation exists:

+ =
Depending on if the first order binary is an asset binary, a bond binary or a Q-option, we can define the corresponding second order binaries.We derive the analytical formulas for the second order binaries in the following theorem.
Proposition 3. Assume the processes for the underlying asset price are given by ( 1)-( 3).Let T be a Lévy subordinator with Lévy exponent φ .Then, where functions s n f is given in Proposition 2 and n λ and ( ) n y ψ are in (10) and (11), respectively.Furthermore, 2) where .

Examples
Once we obtain the analytical formulas for the first and second order binary options, we can apply the theoretical framework for pricing dual-expiry options developed by [1].We can express the payoffs of dual-expiry exotics in terms of a portfolio of elementary binary options.According to the principle of static replication, if the payoff of a European style derivative can be expressed as a portfolio of elementary contacts, then the arbitrary free price of the derivative is the present value of this portfolio.Hereby, we utilize the results of the previous section to the pricing of some dual-expiry exotics.

Chooser Options
These exotics give the holder at time 0 t , the choice of either a European call option of strike c k and expiry c t , or a European put option of strike p k and expiry p t .The payoff at time 0 t , conditioning on ( ) ( ) where C and P are the call and put options, respectively.And ( ) We can calculate time t value of chooser options from the following result.Lemma 3 Assume the processes for the underlying asset price are given by ( 1)-( 3).Let T be a Lévy subordinator with Lévy exponent φ .Then, the time t value of chooser option, conditioning on ( ) , , ;  , , , ; , where d is the unique solution of ( ) , , ; , , ; .
Proof.Since the call option function , , ; is a monotonic decreasing function of ( ) 0 y t , we will have a unique solution d to (32).Then the value of chooser option at time t is given by ( ) { } ( ) ( )

Compound Options
For these exotics, the underlying are options.There are four basic types of compound options, generally referred to as call-on-call, call-on-put, put-on-all and put-on-put options.At time 0 t , the holder of a compound option has the right to buy (or sell) a standard European call (or put) option with strike price t .Then the time t value of a 0 1 s s -compound option can be obtained from the following result.
Lemma 4 Assume the processes for the underlying asset price are given by (1)-(3).Let T be a Lévy subordinator with Lévy exponent φ .Then, the time t value of compound option, conditioning on ( ) where ( ) s s s = and 2 k is the unique solution of ( ) Proof.Since the function ( ) , , , ; , , .

Extendable Options
These exotics allow expiry date to be extended to a future date for a fee with a different strike price at the extended expiry.The holder of an expendable call option has the right at time 0 t to exercise a standard European call option with strike price 0 k ; or for premium p, to extend the expiry date to time 1 0 t t > and change the strike from 0 k to 1 k .The time t value of an extendable option can be obtained using the following result.
Lemma 5. Assume the processes for the underlying asset price are given by ( 1)-( 3).Let T be a Lévy subordinator with Lévy exponent φ .Then the time t value of extendable option, conditioning on ( ) where ( ) and 1 I and 2 I are the solutions of ( ) ( ) , , ; .

Numerical Analysis
In this section, we numerically study chooser, compound and extendable options Journal of Mathematical Finance based on a specific time-changed process.We assume the time change process T is a Gamma subordinator with the Lévy exponent ( ) To calculate the prices for dual exotics, we need to truncate the eigenfunction expansion after a finite number of terms.Following [14], we truncate the infinite series when a given error tolerance level is reached.In practice, we find the convergence of the expansion is rather fast.
In Figures 1-3, we perform some sensitivity tests to demonstrate how sensitive the selected dual-expiry exotics are to the changes in the key parameters in the model.We can summarize the findings as follows: • The option prices decrease with the mean reverting parameter κ for chooser and compound options, but increase for extendable options.between conditional variance and mean reversion speed.The subordination complicates the issue further.Therefore, it is not surprising to find the impacts of model parameters differ across the option types.

Conclusion
This paper studies a new class of models for pricing dual-expiry exotic options.
The underlying asset price is modeled as exponential function of OU process time changed by a Lévy subordinator.The resulting asset price can display both mean reversion and jumps often observed in a large range of underlying assets of exotic option contracts.We employ the method of [1] to decompose the exotic option prices into a portfolio of first and second order binaries.We are able to employ the eigenfunction expansion technique to derive the analytical pricing formulas for the binaries.After that, we can compute dual-expiry exotics option prices by static replication.We also implement the model and analyze the effect

∫
An important sub-class of Lévy subordinators are the tempered stable subordinators.For such subordinators, the Lévy measure ( ) and 0 η > .For such subordinators, the Lévy exponent is given by the eigenvalues of − and { } n ψ are the corresponding eigenfunctions satisfying the following Sturm-Liouville equation . s

1 k and expiry 1 t
for + − for call or put at time 0 t and 1