Next Article in Journal
The Relationship between Governance Quality and the Cost of Equity Capital in Italian Listed Firms: An Update
Next Article in Special Issue
A Computational Approach to Sequential Decision Optimization in Energy Storage and Trading
Previous Article in Journal
Crisis and the Role of Money in the Real and Financial Economies—An Innovative Approach to Monetary Stimulus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pricing Perpetual American Put Options with Asset-Dependent Discounting

by
Jonas Al-Hadad
and
Zbigniew Palmowski
*,†
Faculty of Pure and Applied Mathematics, Wrocław University of Science and Technology, ul. Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Risk Financial Manag. 2021, 14(3), 130; https://doi.org/10.3390/jrfm14030130
Submission received: 26 February 2021 / Revised: 10 March 2021 / Accepted: 16 March 2021 / Published: 20 March 2021
(This article belongs to the Special Issue Stochastic Control and Optimization with Financial Applications)

Abstract

:
The main objective of this paper is to present an algorithm of pricing perpetual American put options with asset-dependent discounting. The value function of such an instrument can be described as V A Put ω ( s ) = sup τ T E s [ e 0 τ ω ( S w ) d w ( K S τ ) + ] , where T is a family of stopping times, ω is a discount function and E is an expectation taken with respect to a martingale measure. Moreover, we assume that the asset price process S t is a geometric Lévy process with negative exponential jumps, i.e., S t = s e ζ t + σ B t i = 1 N t Y i . The asset-dependent discounting is reflected in the ω function, so this approach is a generalisation of the classic case when ω is constant. It turns out that under certain conditions on the ω function, the value function V A Put ω ( s ) is convex and can be represented in a closed form. We provide an option pricing algorithm in this scenario and we present exact calculations for the particular choices of ω such that V A Put ω ( s ) takes a simplified form.
MSC:
60G40; 60J60; 91B28
JEL Classification:
G13; C61

1. Introduction

In this paper, we consider a perpetual American put option with asset-dependent discounting. We consider a standard stochastic background for this problem, i.e., we define a complete filtered risk-neutral probability space ( Ω , F , { F t } t 0 , P ) , on which we define the asset price process S t . Then, F t is a natural filtration of S t satisfying the usual conditions and P is a risk-neutral measure under which the discounted (with respect to the risk-free interest rate r > 0 ) asset price process e r t S t is a local martingale. A family of F t -stopping times is denoted by T while E s denotes the expectation with respect to P when S 0 = s = e x . The value function of the perpetual American put option with asset-dependent discounting can be represented by
V A Put ω ( s ) : = sup τ T E s e 0 τ ω ( S w ) d w ( K S τ ) + .
The asset-dependent discounting is reflected in the ω function, which is our key concept considered in this article. We underline here that the discount function ω for various economical reasons can be different from the risk-free interest rate r > 0 ; see Al-Hadad and Palmowski (2021) for further explanation. The way we choose discounting is to model strong dependence of discount factor with the asset price. The goal is to understand various economical phenomena that might appear in this extreme case. Our approach differs from typical studies considered in the literature, where the interest rate is independent from the asset price or there is a weak dependence between these two factors. Therefore, this research is noteworthy not only in the context of option pricing, but also in other areas where optimisation problems appear.
Our main motivation comes from the work of Linetsky (1999) who considers up-and-out put options traded mainly over-the-counter. In these options, the discounting factor depends exponentially on the cumulative excursion time above or below a given barrier during the entire life of the option. In other words, ω ( s ) = ϱ 𝟙 { s A ( H ) } where A ( H ) = { s 0 : ± ( s H ) 0 } . Similar step-like options were considered in Rodosthenous and Zhang (2018) where the discounting is related to the occupation time of the asset price below a fixed level y exceeding an independent exponential random variable with mean 1 / ϱ . This discounting corresponds to the special choice of ω ( s ) = r + ϱ 𝟙 { s y } , where r is a risk-free interest rate. In general, one can think of ω ( S t ) as the general knock-out (or knock-in) rate depending on the asset price of some vanilla options. This additional feature might allow us to customise the financial product according to the demand of the risk management.
Moreover, we assume that the asset price process S t is a geometric Lévy process with negative exponential jumps, i.e.,
S t : = s e X t
with
X t : = ζ t + σ B t i = 1 N t Y i ,
where ζ and σ > 0 are constant, N t is the Poisson process with intensity λ 0 independent of Brownian motion B t , and { Y i } i N are i.i.d. random variables independent of B t and N t having exponential distribution with mean 1 / φ > 0 . Under the martingale measure P , the drift parameter is of the form
ζ : = r σ 2 2 + λ φ + 1 .
Note that when λ = 0 , then we end up with the classical Black–Scholes model. However, empirical studies show that stock prices have a heavier left tail than normal distribution. Note that the negative jumps underline the market crashes that appear once per while. Therefore, nowadays many books and articles concern, as we do in this work, pricing of derivative securities in market models based on Lévy processes; see Cont and Tankov (2004) for more details. Still, to get a complete picture, we analyse here the case of λ = 0 as well. Our model of asset price is a particular case of the Kou model where the logarithm of the price process is a jump diffusion; see Kou (2002).
The main objective of this paper is to present an algorithm of pricing perpetual American put options with asset-dependent discounting with the value function defined in (1) and the asset price process S t given in (2). Furthermore, we take into account some specific scenarios (e.g., when σ = 0 or λ = 0 ), and, for these cases, we are able to derive analytical forms of the value function, while for more complex examples, we show how to handle them numerically.
Detailed theoretical results of the analysed problem were already developed in Al-Hadad and Palmowski (2021), where the authors presented the approach of deriving a closed form of value function (1) for even a more general setting than it is considered here. Therefore, in this paper, we focus more on numerical side of this problem and analyse, in detail, a few particular cases where more explicit results can be derived.
Still, before we present the option pricing method in our set-up, we recall the most important theoretical issues on which our article is based on. A key step in deriving a closed form of (1) is identifying the form of the optimal stopping rule τ * for which the supremum in (1) is attained. It turns out that under certain conditions on the discount function ω , which are presented in the next section, the value function is convex. By combining this fact with the classical optimal stopping theory presented, e.g., in Peskir and Shiryaev (2006), it allows us to conclude that the optimal stopping region is an interval [ l * , u * ] and hence
τ * = inf { t 0 : S t [ l * , u * ] }
for some optimal thresholds l * u * . Observe that for the non-negative discount function ω , we have l * = 0 (since waiting is not beneficial). Therefore, in this case a single continuation region appears. In general, for the negative ω , we can observe a double continuation region; for more details, see De Donno et al. (2020).
The optimal boundary levels l * and u * can be found by application of standard methods of maximising the function
v A Put ω ( s , l , u ) : = E s e 0 τ l , u ω ( S w ) d w ( K S τ l , u ) +
over l and u > l . To find v A Put ω ( s , l , u ) , we use exit identities for spectrally negative Lévy processes containing so-called omega scale functions introduced in Li and Palmowski (2018).
As shown in (Al-Hadad and Palmowski 2021, Theorem 9), another way of finding the optimal thresholds l * < u * is to apply the classical smooth and continuous fit conditions.
Typically, a price of the option is a solution to a certain Hamiltonian–Jacobi–Bellman (HJB) system, and the optimal thresholds are identified using the smooth fit conditions. We want to underline that our approach is different, although still finding the omega scale functions is carried out via solving certain ordinary differential equations.
The paper is organised as follows. In Section 2, we introduce basic theory and notation. Section 2.4 provides the main theoretical results of this paper. In Section 3, we present some specific examples where the option price can be expressed in the explicit way. Section 4 focuses on the purely numerical analysis. We also show there that these two approaches are consistent. The last section includes our conclusions.

2. Preliminaries

2.1. Assumptions

It the beginning, we note that the analysed American put option will not be realised when its payoff is equal to 0. Hence, we can transform the form of the value function given in (1) into the following one
V A Put ω ( s ) : = sup τ T E s e 0 τ ω ( S w ) d w ( K S τ ) .
We work under the same assumptions as those formulated in (Al-Hadad and Palmowski 2021, Section 2.2). However, this time we consider slightly more specific assumptions on the ω function, namely, that
Assumption 1.
A discount function ω is concave, nondecreasing and bounded from below.
From (Al-Hadad and Palmowski 2021, Remark 3) we can then conclude that under Assumption 1, the value function V A Put ω ( s ) is convex.

2.2. Optimal Stopping Time

From (Al-Hadad and Palmowski 2021, Section 2.4), it follows that the optimal exercise time is the first entrance of the process S t into some interval, that is, it has the following form
τ l , u = inf { t 0 : S t [ l , u ] } .
Hence, we can represent value function (4) as
V A Put ω ( s ) = v A Put ω ( s , l * , u * ) ,
where
v A Put ω ( s , l * , u * ) : = sup 0 l u K v A Put ω ( s , l , u )
and
v A Put ω ( s , l , u ) : = E s e 0 τ l , u ω ( S w ) d w ( K S τ l , u ) .
Moreover, we denote the optimal stopping time by
τ * : = τ l * , u * ,
where l * and u * realise the supremum above. As shown in (Al-Hadad and Palmowski 2021, Theorem 9), another way of identifying the critical points l * and u * can be achieved via application of the smooth fit property. In that case, l * and u * satisfy
( V A Put ω ) ( u * ) = 1 and ( V A Put ω ) ( l * ) = 1 .

2.3. Scale Functions

By applying the fluctuation theory of Lévy processes, we can find a closed form of (5), and hence of (4), in terms of the so-called omega scale functions.
To introduce them formally, first let us define the Laplace exponent of X t via
ψ ( θ ) : = 1 t log E [ e θ X t X 0 = 0 ] ,
which is well-defined for θ 0 since our X t is a spectrally negative Lévy process. In the case of X t given in (3), the Laplace exponent takes the form
ψ ( θ ) = ζ θ + σ 2 2 θ 2 λ θ φ + θ .
By Φ ( q ) , we denote the right inverse of ψ ( θ ) , i.e.,
Φ ( q ) : = sup { θ 0 : ψ ( θ ) = q } ,
where q 0 .
The first scale function W ( q ) ( x ) is defined as a continuous and increasing function such that W ( q ) ( x ) = 0 for all x < 0 , while for x 0 , it is defined via the following Laplace transform
0 e θ x W ( q ) ( x ) d x = 1 ψ ( θ ) q
for θ > Φ ( q ) . We define also the related scale function Z ( q ) ( x ) by
Z ( q ) ( x ) : = 1 + q 0 x W ( q ) ( y ) d y ,
where x R . From Cohen et al. (2013), we know that for X t given in (3), we have
W ( q ) ( x ) = e γ 1 x ψ ( γ 1 ) + e γ 2 x ψ ( γ 2 ) + e Φ ( q ) x ψ ( Φ ( q ) ) ,
where { γ 1 , γ 2 , Φ ( q ) } is the set the real solutions to ψ ( θ ) = q . In turn, Z ( q ) ( x ) is as follows
Z ( q ) ( x ) = 1 + q e γ 1 x 1 γ 1 ψ ( γ 1 ) + e γ 2 x 1 γ 2 ψ ( γ 2 ) + e Φ ( q ) x 1 Φ ( q ) ψ ( Φ ( q ) ) .
If we take σ = 0 or λ = 0 in (7), then W ( q ) ( x ) and Z ( q ) ( x ) take simplified forms
W ( q ) ( x ) = e γ 1 x ψ ( γ 1 ) + e γ 2 x ψ ( γ 2 )
and
Z ( q ) ( x ) = 1 + q e γ 1 x 1 γ 1 ψ ( γ 1 ) + e γ 2 x 1 γ 2 ψ ( γ 2 )
for γ 1 and γ 2 being again the real solutions to ψ ( θ ) = q .
The generalisation of W ( q ) ( x ) and Z ( q ) ( x ) are the ξ -scale functions { W ( ξ ) ( x ) , x R } , { Z ( ξ ) ( x ) , x R } , where ξ is an arbitrary measurable function. They are defined as the unique solutions to the following equations
W ( ξ ) ( x ) = W ( x ) + 0 x W ( x y ) ξ ( y ) W ( ξ ) ( y ) d y ,
Z ( ξ ) ( x ) = 1 + 0 x W ( x y ) ξ ( y ) Z ( ξ ) ( y ) d y ,
where W ( x ) = W ( 0 ) ( x ) is a classical zero scale function.
To simplify notation, we introduce also the following S t counterparts of the scale functions (11) and (12)
W ( ξ ) ( s ) : = W ( ξ exp ) ( log s ) ,
Z ( ξ ) ( s ) : = Z ( ξ exp ) ( log s ) ,
where ξ exp ( x ) : = ξ ( e x ) .
For α , for which the Laplace exponent is well-defined, we can define a new probability measure P ( α ) via
d P s ( α ) d P s F t = e α ( X t log s ) ψ ( α ) t .
By Palmowski and Rolski (2002) and (Kyprianou 2006, Cor. 3.10), under P ( α ) , the process X t is again a spectrally negative Lévy process with the new Laplace exponent
ψ ( α ) ( θ ) : = ψ ( θ + α ) ψ ( α ) = ζ ( α ) θ + σ ( α ) 2 2 θ 2 λ ( α ) θ φ ( α ) + θ ,
where
ζ ( α ) : = ζ + σ 2 α , σ ( α ) : = σ , λ ( α ) : = λ φ φ + α and φ ( α ) : = φ + α .
For the new probability measure P ( α ) , we can define the ξ -scale functions which are denoted by the adding subscript α to the regular counterparts, i.e., W α ( ξ ) ( s ) , Z α ( ξ ) ( s ) .
Lastly, we define the following auxiliary functions
ω u ( s ) : = ω ( s u ) and ω u α ( s ) : = ω u ( s ) ψ ( α ) .

2.4. Theoretical Representation of the Price

The starting point for our entire analysis is the following results. The first one is a corollary from (Al-Hadad and Palmowski 2021, Theorem 15).
Theorem 1.
Let Assumption 1 hold, and assume that ω is non-negative. Then, the optimal stopping region is of the form ( 0 , u * ] and we have
1. 
For σ = 0 and λ > 0
V A Put ω ( s ) : = sup u > 0 v A Put ω ( s , 0 , u ) = sup u > 0 K u φ φ + 1 Z ( ω u ) s u c Z ( ω ) / W ( ω ) W ( ω u ) s u .
2. 
For λ = 0 and σ > 0
V A Put ω ( s ) : = sup u > 0 v A Put ω ( s , 0 , u ) = sup u > 0 ( K u ) lim α s u α Z α ( ω u α ) s u c Z α ( ω α ) / W α ( ω α ) W α ( ω u α ) s u .
3. 
For σ > 0 and λ > 0
V A Put ω ( s ) : = sup u > 0 v A Put ω ( s , 0 , u ) = sup u > 0 { K u φ φ + 1 Z ( ω u ) s u c Z ( ω ) / W ( ω ) W ( ω u ) s u + ( K u ) lim α s u α Z α ( ω u α ) s u c Z α ( ω α ) / W α ( ω α ) W α ( ω u α ) s u } ,
where
c Z ( ω ) / W ( ω ) = lim z Z ( ω ) ( z ) W ( ω ) ( z ) and c Z α ( ω α ) / W α ( ω α ) = lim z Z α ( ω α ) ( z ) W α ( ω α ) ( z ) .
The optimal boundary u * in (18) and (19) can be determined by the smooth fit condition
( V A Put ω ) ( u * ) = 1 ,
while the optimal boundary u * in (17) can be determined by the continuous fit condition
V A Put ω ( u * ) = K u * .
Remark 1.
Let us note that, using (13) and (14), we can interpret the value functions occurring in Theorem 1 both as the functions of s variable and x variable, where x = log s .
Remark 2.
Formula (18) corresponds to the continuous transition of S t to an interval ( 0 , u ] , while (17) describes the situation when S t jumps from ( u , ) into ( 0 , u ] . A combination of these two components makes up formula (19).
To find the option price V A Put ω ( s ) , we have to identify
W α ( ω u α ) ( s ) = W α ( ω u α exp ) ( log s ) and Z α ( ω u α ) ( s ) = Z α ( ω u α exp ) ( log s ) ,
where ω u α ( s ) is given in (16), and the case of α = 0 corresponds to W ( ω u ) ( s ) and Z ( ω u ) ( s ) .
Observe that we need to find the above ξ -scale functions for ξ = ω u α exp under measure P ( α ) , i.e., we have to identify W α ( ξ ) ( x ) and Z α ( ξ ) ( x ) . This is equivalent to taking our asset price process S t of the form of (2) but with the new parameters given in (15).
The second key result for our numerical analysis follows straightforward from (Al-Hadad and Palmowski 2021, Theorem 16) and allows us to identify the above omega scale functions using ordinary differential equations.
Theorem 2.
Assume that ξ is continuously differentiable. Then
1. 
If σ = 0 and λ > 0 or λ = 0 and σ > 0 , then W ( ξ ) ( x ) solves
W ( ξ ) ( x ) = ( ( Υ 1 + Υ 2 ) ξ ( x ) + γ 2 ) W ( ξ ) ( x ) + ( ( Υ 1 + Υ 2 ) ξ ( x ) Υ 1 γ 2 ξ ( x ) ) W ( ξ ) ( x )
with
W ( ξ ) ( 0 ) = Υ 1 + Υ 2 , W ( ξ ) ( 0 ) = Υ 2 γ 2 + ( Υ 1 + Υ 2 ) 2 ξ ( 0 ) .
Moreover, Z ( ξ ) ( x ) solves the same Equation (20) with
Z ( ξ ) ( 0 ) = 1 , Z ( ξ ) ( 0 ) = ( Υ 1 + Υ 2 ) ξ ( 0 ) .
2. 
If σ > 0 and λ > 0 , then W ( ξ ) ( x ) solves
W ( ξ ) ( x ) = γ 2 + γ 3 W ( ξ ) ( x ) + Υ 2 ( γ 2 γ 3 ) ξ ( x ) γ 2 γ 3 Υ 1 γ 3 ξ ( x ) W ( ξ ) ( x ) + Υ 2 ( γ 2 γ 3 ) ξ ( x ) + Υ 1 γ 2 γ 3 ξ ( x ) Υ 1 γ 3 ξ ( x ) W ( ξ ) ( x )
with
W ( ξ ) ( 0 ) = 0 , W ( ξ ) ( 0 ) = Υ 2 γ 2 + Υ 3 γ 3 , W ( ξ ) ( 0 ) = Υ 2 γ 2 2 + Υ 3 γ 3 2 .
Moreover, Z ( ξ ) ( x ) solves the same Equation (21) with
Z ( ξ ) ( 0 ) = 1 , Z ( ξ ) ( 0 ) = 0 , Z ( ξ ) ( 0 ) = ξ ( 0 ) ( Υ 2 ( γ 2 γ 3 ) Υ 1 γ 3 ) .

3. Option Pricing—Analytical Approach

In this section, we present some examples of discount functions for which we are able to determine the analytical form of the value function V A Put ω ( s ) .

3.1. Constant Discount Function

The case when ω function is constant, i.e., ω ( s ) = q is the standard example which appears in the literature quite extensively. However, this case is quite special, as it turns out that the second term of the sum in (19) simplifies and we do not need to deal with the measure P ( α ) (and thus to calculate the limit for α ) to find V A Put ω ( s ) . This fact is stated in the below theorem.
Theorem 3.
Assume that ω ( s ) = q . Then
lim α s u α Z α ( q ψ ( α ) ) s u c Z α ( q ψ ( α ) ) / W α ( q ψ ( α ) ) W α ( q ψ ( α ) ) s u = σ 2 2 W ( q ) s u Φ ( q ) W ( q ) s u .
Proof. 
Note that
lim α s u α Z α ( q ψ ( α ) ) s u c Z α ( q ψ ( α ) ) / W α ( q ψ ( α ) ) W α ( q ψ ( α ) ) s u
corresponds to the continuous transition of the process S t to the interval ( 0 , u ] or, in other words, continuous exit from half-line ( u , ) . We define
τ b + = inf { t 0 : X t b } , τ 0 = inf { t 0 : X t 0 } .
It turns out that formula (23) is equivalent to
E s u e q τ 0 ; τ 0 < τ b + , X τ 0 = 0 ,
for details, see (Al-Hadad and Palmowski 2021, Proof of Theorem 7).
Then, using (Loeffen et al. 2014, (13)) for x = log s u , a = 0 and v ( q ) ( x ) = W ( q ) ( x ) together with the fact that W ( q ) ( 0 ) = 2 σ 2 (see (Kyprianou 2006, 8.5 (ii), p. 235)), we obtain
E s u e q τ 0 ; τ 0 < τ b + , X τ 0 = 0 = σ 2 2 W ( q ) ( x log u ) W ( q ) ( b ) W ( q ) ( b ) W ( q ) ( x log u ) .
Lastly, taking limit b and applying L’Hospital’s Rule, complete the proof. □
Ultimately, value function (19) for the constant discount function ω ( s ) = q can be written as
V A Put ω ( s ) = sup u > 0 { K u φ φ + 1 Z ( q ) s u c Z ( q ) / W ( q ) W ( q ) s u + ( K u ) σ 2 2 W ( q ) s u Φ ( q ) W ( q ) s u } .
Remark 3.
For the case of λ = 0 , using (22), one can show that value function (18) simplifies to the well-known formula for the value function in the Black–Scholes model, i.e.,
V A Put ω ( s ) = sup u > 0 ( K u ) s u 2 r σ 2 ,
where we substituted q = r and ζ = r σ 2 2 . Therefore, we are not forced to apply the smooth fit condition in order to find the optimal value of u. We can do this analytically by finding the maximum of V A Put ω ( s ) with respect to u.
Figure 1 presents value function (19) for three different values of q, i.e., q { 0.3 , 0.6 , 0.9 } .
Based on Figure 1, we can simply note that a higher value of the discount function ω results in a smaller value of V A Put ω ( s ) , which is in line with the financial intuition.
In turn, Figure 2 shows a comparison of (17)–(19) for the same value of q = 0.5 .
The resulting relation between these functions is again consistent with the economical expectations.

3.2. Linear Discount Function

In this subsection, we consider a linear discount function of the form ω ( s ) = C s for some positive constant C.

3.2.1. σ = 0

Let us consider the case of σ = 0 . Then, the asset price process S t jumps into the interval ( 0 , u ] , which means from Theorem 1 that
V A Put ω ( s ) = sup u > 0 K u φ φ + 1 Z ( ω u ) s u c Z ( ω ) / W ( ω ) W ( ω u ) s u ,
where ω u s u = ω ( s ) = C s . Equivalently, (24) can be rewritten as
V A Put η ( x ) = sup u > 0 K u φ φ + 1 Z ( η u ) x log u c Z ( η ) / W ( η ) W ( η u ) x log u ,
where x = log s and η u ( x log u ) = η ( x ) = C e x .
To find a closed form of value function (25), we need to identify the scale functions W ( η u ) ( x log u ) and Z ( η u ) ( x log u ) . From Theorem 2, it follows that both W ( η ) ( x ) and Z ( η ) ( x ) solve the following ordinary differential equation
f ( x ) = ( A e x + B ) f ( x ) + D e x f ( x )
with A = C ζ , B = λ φ ζ ζ and D = C 1 + φ ζ , while the initial conditions are as follows
W ( η ) ( 0 ) = 1 ζ , W ( η ) ( 0 ) = C + λ ζ 2
and
Z ( η ) ( 0 ) = 1 , Z ( η ) ( 0 ) = C ζ .
Substituting t = A e x and F ( t ) = f ( x ) to (26), we obtain the Kummer’s equation of the form
t F ( t ) + ( b t ) F ( t ) a F ( t ) = 0 ,
where b = 1 B and a = D A .
If b is not an integer, then the general solution to (29) has the form
F ( t ) = K 1 1 F 1 ( a 1 , b 1 ; t ) + K 2 t 1 b 1 F 1 ( a 2 , b 2 ; t ) ,
where K 1 and K 2 are the constants that can be found based on the initial conditions, a 1 = a , b 1 = b , a 2 = a b + 1 , b 2 = 2 b , while 1 F 1 ( · , · ; · ) is the Kummer confluent hypergeometric function.
We denote by K 1 W , K 2 W and K 1 Z , K 2 Z the constants corresponding to W ( η ) ( x ) and Z ( η ) ( x ) , respectively. Using initial conditions (27) and (28), we can simply calculate these constants for both W ( η ) ( x ) and Z ( η ) ( x ) . By shifting these functions by log u , we simply produce W ( η u ) ( x log u ) and Z ( η u ) ( x log u ) .
The asymptotic behaviour of 1 F 1 a , b ; t for t is as follows
1 F 1 ( a , b ; t ) = Γ ( b ) Γ ( a ) e t t a b 1 + O 1 t .
Based on (31), we calculate the constant c Z ( η ) / W ( η ) (or equivalenly c Z ( ω ) / W ( ω ) ) occuring in (25). It has the following form
c Z ( η ) / W ( η ) = K 1 Z Γ ( b 1 ) Γ ( a 1 ) A a 1 b 1 + K 2 Z Γ ( b 2 ) Γ ( a 2 ) A a 2 1 K 1 W Γ ( b 1 ) Γ ( a 1 ) A a 1 b 1 + K 2 W Γ ( b 2 ) Γ ( a 2 ) A a 2 1 .
Combining all the obtained results and substituting them into (25), we can present value function (25) graphically for some sample parameter values. This is carried out in Section 4.

3.2.2. λ = 0

Let us consider the case of λ = 0 . In this case, the asset price process S t enters the interval ( 0 , u ] in a continuous way only. Therefore, from Theorem 1
V A Put ω ( s ) = sup u > 0 ( K u ) lim α s u α Z α ( ω u α ) s u c Z α ( ω α ) / W α ( ω α ) W α ( ω u α ) s u
which is equivalent to
V A Put η ( x ) = sup u > 0 ( K u ) lim α e α ( x log u ) Z α ( η u α ) ( x log u ) c Z α ( η α ) / W α ( η α ) W α ( η u α ) ( x log u ) .
It suffices to find now W α ( η u α ) ( x log u ) and Z α ( η u α ) ( x log u ) . From Theorem 2, it follows that W α ( η α ) ( x ) and Z α ( η α ) ( x ) solve
f ( x ) = B α f ( x ) + ( D α e x + E α ) f ( x )
with B α = 2 σ 2 ( ζ + σ 2 α ) , D α = 2 C σ 2 and E α = 2 σ 2 ζ α + σ 2 2 α 2 . The initial conditions have the following form
W α ( η α ) ( 0 ) = 0 , W α ( η α ) ( 0 ) = 2 σ 2
and
Z α ( η u α ) ( 0 ) = 1 , Z α ( η u α ) ( 0 ) = 0 .
Substituting t = 2 D α e x and F ( t ) = e B α x 2 f ( x ) to (35), we obtain the Bessel differential equation of the form
t 2 F ( t ) + t F ( t ) + ( t 2 v 2 ) F ( t ) = 0 ,
where v = B α 2 + 4 E α = 2 ζ σ 2 . The general solution to (38) is equal to
F ( t ) = K 1 J v ( t ) + K 2 Y v ( t )
and therefore
f ( x ) = e B α x 2 K 1 J v ( 2 D α e x ) + K 2 Y v ( 2 D α e x ) .
Based on the form of (39) and the fact that D α does not depend on α , we can simply note that value function (34) is also independent of α . Therefore, we can take an arbitrary value of α in (35). Thanks to this key observation, its solution (39) can take a simplified form. Indeed, for α = 0 , Equation (35) is equal to
f ( x ) = B 0 f ( x ) + D 0 e x f ( x ) ,
where B 0 = 2 ζ σ 2 and D 0 = 2 C σ 2 . Hence, the general solution to (40) takes the following form
f ( x ) = e B 0 x 2 K 1 J v ( 2 D 0 e x ) + K 2 Y v ( 2 D 0 e x ) .
For B 0 = 1 2 n , where n N 0 and D 0 t > 0 , Equation (41) reduces to
f ( x ) = K 1 cosh ( 4 D 0 e x ) n + K 2 sinh ( 4 D 0 e x ) n .
If we take the following sample parameters r = 0.05 and σ = 0.2 , we obtain n = 2 and therefore
f ( x ) = K 1 3 sinh ( 2 e x ) 4 e 5 2 x + sinh ( 2 e x ) e 3 2 x 3 cosh ( 2 e x ) 2 e 2 x + K 2 3 cosh ( 2 e x ) 4 e 5 2 x + cosh ( 2 e x ) e 3 2 x 3 sinh ( 2 e x ) 2 e 2 x .
Applying initial conditions (36) and (37), we can simply obtain K 1 W , K 2 W and K 1 Z , K 2 Z . Using equality (42), which holds for both W α ( η α ) ( x ) and Z α ( η α ) ( x ) , we can calculate that
c Z α ( η α ) / W α ( η α ) = c Z ( η ) / W ( η ) = K 1 Z + K 2 Z K 1 W + K 2 W .
Taking into account all the obtained results, we can obtain value function (34) for sample data. This is performed in Section 4 as well.

3.3. Power Discount Function

This time, we take into account a power function of the form ω ( s ) = C s n for n ( 0 , 1 ] and C being some positive contant. This case is a generalisation of a linear discount function.

3.3.1. σ = 0

Similar to the case of a linear discount function, the scale functions W ( η ) ( x ) and Z ( η ) ( x ) solve
f ( x ) = ( A e n x + B ) f ( x ) + D e n x f ( x )
with A = C ζ , B = λ φ ζ ζ and D = C n + φ ζ , while the initial conditions are the same as those provided in (27) and (28). Applying a substitution t = A n e n x and F ( t ) = f ( x ) , we transform (44) into
t F ( t ) + ( b t ) F ( t ) a F ( t ) = 0 ,
where b = 1 B n and a = D A n . The general solution to (45) has the same form as was provided in (30).
Therefore, for both the linear and the power discount function ω , the form of the value function V A Put η ( x ) is identical.

3.3.2. λ = 0

As in the above case, the idea of finding a closed form of the value function can be borrowed from the linear case. This time, the scale functions W α ( η α ) ( x ) and Z α ( η α ) ( x ) satisfy the equation
f ( x ) = B α f ( x ) + ( D α e n x + E α ) f ( x )
with B α = 2 σ 2 ( ζ + σ 2 α ) , D α = 2 C σ 2 and E α = 2 σ 2 ζ α + σ 2 2 α 2 , while the initial conditions are of the form (36) and (37). If we substitute t = 2 n D α e n x and F ( t ) = e B α x 2 f ( x ) , we receive the Bessel differential equation for F ( t ) with the solution
F ( t ) = K 1 J v ( t ) + K 2 Y v ( t ) .
Therefore, we have
f ( x ) = e B α x 2 K 1 J v 2 n D α e n x + K 2 Y v 2 n D α e n x ,
where v = B α 2 + 4 E α n = 2 ζ n σ 2 . We can show, as in the prevoius section, that the value function which arises in this scenario does not depend on α . Thus, for α = 0 , (46) takes the form
f ( x ) = e B 0 x 2 K 1 J v 2 n D 0 e n x + K 2 Y v 2 n D 0 e n x .
Again, having exact formulas for the scale functions, we can easily represent the form of the value function. All results are presented in Section 4.

4. Option Pricing—Numerical Approach

In this section, we show how to numerically identify the value function V A Put ω ( s ) for arbitrary discount function ω . We present some figures corresponding to the various discount functions as well.

4.1. Different Discount Functions

For some discount functions ω , we are unable to find the analytical forms of W ( η ) ( x ) , Z ( η ) ( x ) , W α ( η α ) ( x ) and Z α ( η α ) ( x ) , which are the solutions to the ordinary differential equations occurring in Theorem 2. That is, formally we cannot identify explicitly the value function either. In such a situation, we can proceed a numerical analysis of these equations.
In general, solving a high-order ordinary differential equation consists of transforming it into first-order vector form and then applying an appropriate algorithm that returns the numerical solution of the n + 1 -dimensional system of first-order ordinary differential equations. For practical purposes, however—such as in financial engineering—numeric approximations to the solutions of ordinary differential equations are often sufficient. In this paper, we focus on the Higher-Order Taylor Method. This method employs the Taylor polynomial of the solution to the equation. It approximates the 0-th order term by using the previous step’s value (which is the initial condition for the first step), and the subsequent terms of the Taylor expansion by using the differential equation.

4.1.1. σ = 0

For the case of σ = 0 , we showed in the previous section how to derive the value function for the linear discount function ω . Figure 3 presents comparison of the value function given in (24) when the scale functions were calculated analytically (as shown in Section 3.2.1) and numerically by solving differential ordinary Equation (26). In this example, a linear discount function was chosen, i.e., ω ( s ) = C s . The difference between these functions is so small that this is negligible.
Moreover, Figure 4 illustrates the constant c Z ( η ) / W ( η ) obtained in (32) together with the quotient of the functions Z ( η ) ( x ) and W ( η ) ( x ) .
As we mentioned at the beginning of this section, we are not always able to get an analytical solution to an ordinary differential equation. This is the case, for example, when the discount function is of the form ω ( s ) = C arctan ( s ) for some positive C. Then, we can only obtain the scale functions numerically. Figure 5 shows two value functions, for ω ( s ) = C s and ω ( s ) = C arctan ( s ) , respectively. Since for all positive s we have s > arctan ( s ) , then we expect that the value function corresponding to ω ( s ) = C arctan ( s ) takes greater values rather than for ω ( s ) = C s . We can also note that the difference between these functions becomes greater for higher values of s, which is in line with economical intuition since the difference between ω ( s ) = s and ω ( s ) = arctan ( s ) also extends as s increases.

4.1.2. λ = 0

The case of λ = 0 corresponds to the situation when the stock price process (2) does not have any jumps. Then, the value function takes the form (18). From the numerical point of view, the problem lies in choosing a sufficiently large value of α in (18) to obtain the final and right form of the value function. In this section, we avoid this problem by selecting discount functions for which the value function is independent of the α parameter.
Figure 6 presents a comparison of the value function given in (33) for ω ( s ) = C s and for the scale functions obtained analytically and numerically. As shown in Section 3.2.2, in this case, we can take an arbitrary value of α and obtain a simplified form of the value function and ordinary differential equation that the scale functions solve. Similar to the previous example, we again can observe a negligible difference between these two value functions.
In turn, Figure 7 shows the constant c Z ( η ) / W ( η ) given in (43) with the ratio of Z ( η ) ( x ) and W ( η ) ( x ) .
Lastly, in Figure 8, we can observe the value functions for both ω ( s ) = C s and ω ( s ) = C s + Z for some positive Z, i.e., we compare two discount functions that differ in shift. This time, we can see that the value functions obtained in this way also differ only in shift, which confirms financial intuition.

4.1.3. σ > 0 and λ > 0

The most general case is when σ > 0 and λ > 0 . Then, the considered value function V A Put ω ( s ) is given by (19). It can be also represented as the function of x variable:
V A Put η ( x ) = sup u > 0 { K u φ φ + 1 Z ( η u ) ( x log u ) c Z ( η ) / W ( η ) W ( η u ) ( x log u ) + ( K u ) lim α e α ( x log u ) Z α ( η u α ) ( x log u ) c Z α ( η α ) / W α ( η α ) W α ( η u α ) ( x log u ) } .
For the linear discount function ω ( s ) = C s , the scale functions W ( η ) ( x ) and Z ( η ) ( x ) occurring in (47) are the solutions to the following ordinary differential equation
f ( x ) = A f ( x ) + ( B e x + D ) f ( x ) + E e x f ( x )
with A = γ 2 + γ 3 , B = C Υ 2 ( γ 2 γ 3 ) Υ 1 γ 3 , D = γ 2 γ 3 , E = C [ Υ 2 ( γ 2 γ 3 ) + Υ 1 γ 2 γ 3 Υ 1 γ 3 ] .
The initial conditions for W ( η ) ( x ) and Z ( η ) ( x ) are as follows
W ( η ) ( 0 ) = 0 , W ( η ) ( 0 ) = Υ 2 γ 2 + Υ 3 γ 3 , W ( η ) ( 0 ) = Υ 2 γ 2 2 + Υ 3 γ 3 2
and
Z ( η ) ( 0 ) = 1 , Z ( η ) ( 0 ) = 0 , Z ( η ) ( 0 ) = C Υ 2 ( γ 2 γ 3 ) Υ 1 γ 3 .
Thus, W α ( η α ) ( x ) and Z α ( η α ) ( x ) solve
f ( x ) = A α f ( x ) + ( B α e x + D α ) f ( x ) + ( E α e x + F α ) f ( x )
with A α = γ α 2 + γ α 3 , B α = C Υ α 2 ( γ α 2 γ α 3 ) Υ α 1 γ α 3 , D α = Υ α 2 ( γ α 2 γ α 3 ) ψ ( α ) γ 2 γ 3 + Υ α 1 γ α 3 ψ ( α ) , E α = C Υ α 2 ( γ α 2 γ α 3 ) + Υ α 1 γ α 2 γ α 3 Υ α 1 γ α 3 , F α = Υ α 1 γ α 2 γ α 3 ψ ( α ) . The initial conditions for W α ( η α ) ( x ) and Z α ( η α ) ( x ) are as follows
W α ( η α ) ( 0 ) = 0 , W α ( η α ) ( 0 ) = Υ α 2 γ α 2 + Υ α 3 γ α 3 , W α ( η α ) ( 0 ) = Υ α 2 γ α 2 2 + Υ α 3 γ α 3 2
and
Z α ( η α ) ( 0 ) = 1 , Z α ( η α ) ( 0 ) = 0 , Z α ( η α ) ( 0 ) = C Υ α 2 ( γ α 2 γ α 3 ) Υ α 1 γ α 3 .
In this case, we are not able to identify explicit solutions to third-order ordinary differential Equations (48) and (49). Therefore, we are forced to use only the numerical method to find the scale functions and, hence, the value function.
Figure 9 shows several graphs of the value function (19) for different values of α parameter together with the first and second component occurring in (19).

5. Conclusions

In this paper, we have presented a novel approach to pricing the perpetual American put options with asset-dependent discounting. For the asset price process S t being the geometric spectrally negative Lévy process, we have shown that the value function (1) is the solution of the ordinary differential equations given in Theorem 2. We have used these theoretical results to perform extended numerical analysis for some key financial examples. In particular, for some cases, we have managed to produce some explicit formulas for the value function. In the cases where it was impossible to do so, we have used the numerical analysis of the above-mentioned ordinary differential equations based on the Higher-Order Taylor Method. We have presented many figures of the value functions that arise in various scenarios.
One can think of further generalisations; for example, when the discount factor is randomised. Thia can be achieved in different ways, e.g., by introducing additional Markov economical environment. One can also think about other very realistic models for financial assets, for example, based on generalised hyperbolic distributions and their subclasses; see Prause (1999), for more details. This type of research is left for future investigations.

Author Contributions

Methodology, J.A.-H. and Z.P.; Formal Analysis, J.A.-H. and Z.P.; Investigation, J.A.-H. and Z.P.; Writing—Original Draft Preparation, J.A.-H.; Writing—Review and Editing, Z.P.; Visualization, J.A.-H.; Supervision, Z.P.; Project Administration, Z.P.; Funding Acquisition, Z.P. All authors have read and agreed to the published version of the manuscript.

Funding

This paper is supported by the National Science Centre under the grant 2016/23/B/HS4/00566 (2017–2020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Both authors thank the anonymous referees for their helpful comments.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Al-Hadad, Jonas, and Zbigniew Palmowski. 2021. Perpetual American Options with Asset-Dependent Discounting. Submitted for Publication. Available online: https://arxiv.org/abs/2007.09419 (accessed on 18 July 2020).
  2. Cohen, Serge, Alexey Kuznetsov, Andreas Kyprianou, and Victor Rivero. 2013. Lévy Matters II. Berlin/Heidelberg: Springer. [Google Scholar]
  3. Cont, Rama, and Peter Tankov. 2004. Financial Modelling with Jump Processes. Boca Raton: Chapman & Hall. [Google Scholar]
  4. De Donno, Marzia, Zbigniew Palmowski, and Joanna Tumilewicz. 2020. Double continuation regions for American and Swing options with negative discount rate in Lévy models. Mathematical Finance 30: 196–227. [Google Scholar] [CrossRef] [Green Version]
  5. Kou, Steven G. 2002. A Jump-Diffusion Model for Option Pricing. Management Science 48: 1086–1101. [Google Scholar] [CrossRef] [Green Version]
  6. Kyprianou, Andreas E. 2006. Introductory Lectures on Fluctuations of Lévy Processes with Applications. Berlin/Heidelberg: Springer. [Google Scholar]
  7. Li, Bo, and Zbigniew Palmowski. 2018. Fluctuations of omega–killed spectrally negative Lévy processes. Stochastic Processes and their Applications 128: 3273–99. [Google Scholar] [CrossRef] [Green Version]
  8. Linetsky, Vadim. 1999. Step options. Mathematical Finance 9: 55–96. [Google Scholar] [CrossRef]
  9. Loeffen, Ronnie L., Jean-François Renaud, and Xiaowen Zhou. 2014. Occupation times of intervals until first passage times for spectrally negative Lévy processes. Stochastic Processes and their Applications 124: 1408–143. [Google Scholar] [CrossRef] [Green Version]
  10. Palmowski, Zbigniew, and Tomasz Rolski. 2002. A technique for exponential change of measure for Markov processes. Bernoulli 8: 767–85. [Google Scholar]
  11. Peskir, Goran, and Albert Shiryaev. 2006. Optimal Stopping and Free–Boundary Problems. Basel: Birkhäuser. [Google Scholar]
  12. Prause, Karsten. 1999. The Generalized Hyperbolic Model: Estimation, Financial Derivatives, and Risk Measures. Ph.D. thesis, University of Freiburg, Breisgau, Germany. Available online: https://d-nb.info/961152192/34 (accessed on October 1999).
  13. Rodosthenous, Neofytos, and Hongzhong Zhang. 2018. Beating the omega clock: An optimal stopping problem with random time-horizon under spectrally negative Lévy models. Annals of Applied Probability 28: 2105–40. [Google Scholar] [CrossRef]
Figure 1. The value and payoff functions for the given set of parameters: K = 20 , r = 0.05 , σ = 0.2 , λ = 6 , φ = 2 and q { 0.3 , 0.6 , 0.9 } .
Figure 1. The value and payoff functions for the given set of parameters: K = 20 , r = 0.05 , σ = 0.2 , λ = 6 , φ = 2 and q { 0.3 , 0.6 , 0.9 } .
Jrfm 14 00130 g001
Figure 2. The value and payoff functions for the given set of parameters: K = 20 , r = 0.05 , σ = 0.4 , λ = 6 , φ = 2 and q = 0.5 .
Figure 2. The value and payoff functions for the given set of parameters: K = 20 , r = 0.05 , σ = 0.4 , λ = 6 , φ = 2 and q = 0.5 .
Jrfm 14 00130 g002
Figure 3. Comparison of value function (24) for ω ( s ) = C s for both methods of determining the scale functions—analytical and numerical one. The chosen set of parameters is as follows: K = 20 , C = 0.1 , r = 0.05 , λ = 6 , φ = 2 .
Figure 3. Comparison of value function (24) for ω ( s ) = C s for both methods of determining the scale functions—analytical and numerical one. The chosen set of parameters is as follows: K = 20 , C = 0.1 , r = 0.05 , λ = 6 , φ = 2 .
Jrfm 14 00130 g003
Figure 4. Comparison of the constant c Z ( η ) / W ( η ) and the ratio of Z ( η ) ( x ) and W ( η ) ( x ) for a linear discount function ω ( s ) = C s and K = 20 , C = 0.1 , r = 0.05 , λ = 6 , φ = 2 .
Figure 4. Comparison of the constant c Z ( η ) / W ( η ) and the ratio of Z ( η ) ( x ) and W ( η ) ( x ) for a linear discount function ω ( s ) = C s and K = 20 , C = 0.1 , r = 0.05 , λ = 6 , φ = 2 .
Jrfm 14 00130 g004
Figure 5. Comparison of value function (24) for both ω ( s ) = C s and ω = C arctan ( s ) and for the given set of parameters: K = 20 , C = 0.5 , r = 0.05 , λ = 6 , φ = 2 .
Figure 5. Comparison of value function (24) for both ω ( s ) = C s and ω = C arctan ( s ) and for the given set of parameters: K = 20 , C = 0.5 , r = 0.05 , λ = 6 , φ = 2 .
Jrfm 14 00130 g005
Figure 6. Comparison of value function (33) for ω ( s ) = C s for both methods of determining the scale functions–analytical and numerical one. The chosen set of parameters is as follows: K = 20 , C = 0.1 , r = 0.05 , σ = 0.2 .
Figure 6. Comparison of value function (33) for ω ( s ) = C s for both methods of determining the scale functions–analytical and numerical one. The chosen set of parameters is as follows: K = 20 , C = 0.1 , r = 0.05 , σ = 0.2 .
Jrfm 14 00130 g006
Figure 7. Comparison of the constant c Z ( η ) / W ( η ) and the ratio of Z ( η ) ( x ) and W ( η ) ( x ) for a linear discount function ω ( s ) = C s and K = 20 , C = 0.1 , r = 0.05 , σ = 0.2 .
Figure 7. Comparison of the constant c Z ( η ) / W ( η ) and the ratio of Z ( η ) ( x ) and W ( η ) ( x ) for a linear discount function ω ( s ) = C s and K = 20 , C = 0.1 , r = 0.05 , σ = 0.2 .
Jrfm 14 00130 g007
Figure 8. Comparison of value function (33) for both ω ( s ) = C s and ω = C s + Z and for the given set of parameters: K = 20 , C = 0.005 , Z = 0.1 , r = 0.05 , σ = 0.2 .
Figure 8. Comparison of value function (33) for both ω ( s ) = C s and ω = C s + Z and for the given set of parameters: K = 20 , C = 0.005 , Z = 0.1 , r = 0.05 , σ = 0.2 .
Jrfm 14 00130 g008
Figure 9. Comparison of value function (19) for the particular choice of α , ω ( s ) = C s and for the given set of parameters: K = 20 , C = 0.1 , r = 0.05 , σ = 0.2 , λ = 6 , φ = 2 .
Figure 9. Comparison of value function (19) for the particular choice of α , ω ( s ) = C s and for the given set of parameters: K = 20 , C = 0.1 , r = 0.05 , σ = 0.2 , λ = 6 , φ = 2 .
Jrfm 14 00130 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Al-Hadad, J.; Palmowski, Z. Pricing Perpetual American Put Options with Asset-Dependent Discounting. J. Risk Financial Manag. 2021, 14, 130. https://doi.org/10.3390/jrfm14030130

AMA Style

Al-Hadad J, Palmowski Z. Pricing Perpetual American Put Options with Asset-Dependent Discounting. Journal of Risk and Financial Management. 2021; 14(3):130. https://doi.org/10.3390/jrfm14030130

Chicago/Turabian Style

Al-Hadad, Jonas, and Zbigniew Palmowski. 2021. "Pricing Perpetual American Put Options with Asset-Dependent Discounting" Journal of Risk and Financial Management 14, no. 3: 130. https://doi.org/10.3390/jrfm14030130

APA Style

Al-Hadad, J., & Palmowski, Z. (2021). Pricing Perpetual American Put Options with Asset-Dependent Discounting. Journal of Risk and Financial Management, 14(3), 130. https://doi.org/10.3390/jrfm14030130

Article Metrics

Back to TopTop