American put options with regime-switching volatility

Jang, Bong-Gyu; Koo, Hyeng Keun

doi:10.1108/JDQS-12-2023-0043

We present an approach for pricing American put options with a regime-switching volatility. Our method reveals that the option price can be expressed as the sum of two components: the price of a European put option and the premium associated with the early exercise privilege. Our analysis demonstrates that, under these conditions, the perpetual put option consistently commands a higher price during periods of high volatility compared to those of low volatility. Moreover, we establish that the optimal exercise boundary is lower in high-volatility regimes than in low-volatility regimes. Additionally, we develop an analytical framework to describe American puts with an Erlang-distributed random-time horizon, which allows us to propose a numerical technique for approximating the value of American puts with finite expiry. We also show that a combined approach involving randomization and Richardson extrapolation can be a robust numerical algorithm for estimating American put prices with finite expiry.

1. Introduction

The coming of modern option pricing theory, pioneered by Black and Scholes (1973) and Merton (1973), has had a profound impact on both the academic discipline of finance and the practical operations of the financial services industry. The foundation of this theory, rooted in the concept of the absence of arbitrage opportunities, has evolved into a dominant paradigm in the field of finance, as exemplified by seminal works such as Harrison and Kreps (1979) and Dybvig and Ross (2003). This transformative technology for pricing and risk management has empowered financial institutions to develop a diverse array of customized contracts and securities, tailored precisely to meet the unique requirements of their clientele. In this light, option pricing stands as a fundamental cornerstone of modern finance.

In the early stages of option pricing analysis, exemplified by the renowned Black–Scholes formula, a fundamental assumption prevailed: that the volatility of the underlying asset remained constant. Nevertheless, the veracity of this assumption has been cast into doubt, primarily due to compelling empirical evidence revealing that implied volatility, as observed across various strike prices and maturities, does not exhibit constancy.

In response to the challenges posed by this empirical evidence, extensive research endeavors have been undertaken to expand the scope of the Black–Scholes–Merton framework to encompass scenarios featuring randomly fluctuating volatility. Substantial progress has been made in this pursuit, with notable contributions from scholars such as Hull and White (1987), Wiggins (1987), Stein and Stein (1991), Heston (1993), Bates (1996) and Duffie et al. (2000). These efforts have significantly advanced our understanding and modeling capabilities in the context of dynamic and variable volatility in financial markets.

We contribute to the field of option pricing by conducting an investigation into the pricing of American options in the context of stochastic volatility. Our specific focus on stochastic volatility entails a unique characteristic: the volatility of the underlying asset exhibits a binary behavior. It transitions between two distinct states, denoted as the “high volatility regime” (H) and the “low volatility regime” (L), at the occurrence of jumps dictated by an independent Poisson process. The regime-switching volatility model can be considered as the simplest version of the stochastic volatility model. This model is useful when the underlying asset undergoes regime shifts over a period or when the maturity of the underlying asset is long enough so that the option lifetime includes abrupt environmental changes. Compared to parameter-rich volatility models like the Heston (1993) model, it has relatively fewer parameters, making it considerably easier to calibrate from options market data.

However, it is important to note that our model does not encompass the full spectrum of randomness associated with volatility. For instance, it does not account for potential correlations between changes in volatility and alterations in the underlying asset’s price. Nevertheless, our exploration of American option pricing within the context of regime-switching volatility represents an initial step toward comprehending the valuation of American options in the presence of stochastic volatility. This research sets the foundation for further, more intricate models that can capture additional nuances in the relationship between option pricing and the dynamic nature of volatility.

Option pricing in the context of regime switching has attracted the attention of several researchers, including Naik (1993), Chourdakis and Tzavalis (2000), Chourdakis (2001), Campbell and Li (2002) and Edwards (2005). Empirical findings from their work have shown promise, indicating that regime-switching models offer a fitting representation of market data comparable to more intricate models like stochastic volatility models with jumps. The papers, however, has primarily concentrated on European options.

In a distinctive departure, Bollen (1998) introduced a novel approach known as the “pentanomial tree” to calculate the price of American options within a regime-switching framework. While Bollen presented numerical results, analytical solutions were not derived. Building upon this foundation, Bollen et al. (2000) extended the application of Bollen’s pentanomial trees to investigate currency option prices. Duan et al. (2002) have contributed to this field by proposing a numerical methodology for pricing American options in a discrete-time regime-switching volatility model. Similarly, Driffill et al. (2002) have introduced a closed-form solution for perpetual American call options in scenarios where the dividend process follows a regime-switching pattern. Lastly, Guo and Zhang (2004) have made significant strides by deriving analytic formulas for perpetual American put options within the context of regime-switching volatility. These collective efforts have enriched our understanding of option pricing under regime switching and have paved the way for further exploration in this intriguing area of finance.

We establish fundamental theoretical insights within the framework of regime-switching volatility. Specifically, we demonstrate that the price of an American put option can be represented as the combination of a European put option price and the premium attributed to the early exercise privilege. This characterization extends findings established by Kim (1990), Jacka (1991) and Carr et al. (1992) within the context of constant volatility. Moreover, we provide a rigorous proof that the valuation of an American put option constitutes a unique solution to the free boundary value problem. This result builds upon earlier work by Mckean (1965) and Van Moerbeke (1976), expanding its applicability to the regime-switching volatility setting. Furthermore, we delve into the realm of perpetual American put options, demonstrating that their pricing exhibits distinct characteristics across the volatility regimes. Specifically, we establish that perpetual put options command a higher price during periods of high volatility compared to low volatility phases. Additionally, we reveal that the optimal exercise boundary is situated at a lower threshold in high-volatility regimes, as opposed to low-volatility regimes. These findings contribute valuable insights into the dynamics of option pricing in the context of fluctuating volatility regimes (Theorems 3.7 and 3.8).

To facilitate the numerical determination of the American put option’s value, we employ the randomization method introduced by Carr (1998). This method involves substituting the fixed time horizon with a random time, effectively transforming the problem into one that resembles the pricing of a perpetual put option. Consequently, the pricing of an American put option with a specified time horizon, denoted as T, can be approximated through a backward induction process. This process considers a substantial number, N, of randomly generated time points, whose summation is expected to equal the original time horizon T, with a variance equal to T/N. We derive an analytical expression for the value of an American put option with such a random time horizon. Subsequently, we present numerical examples that showcase the robustness of the randomization method in conjunction with Richardson extrapolation. This combined approach serves as a robust algorithm for approximating the price of American put options with finite maturities in the regime-switching model. Our analysis includes a comparative evaluation of the numerical results obtained through this method with those generated using the pentanomial tree proposed by Bollen (1998). Furthermore, we establish a noteworthy observation: the N-th approximation of an American put option with finite maturity consistently commands a higher price during periods of high-volatility compared to low-volatility phases. Additionally, we demonstrate that the optimal exercise boundary is situated at a lower threshold in high-volatility regimes relative to low-volatility regimes (Theorem 4.4) ^[1].

A paper closely related to ours is Boyarchenko and Levendorskii (2009), which presents an American option pricing framework within the context of Markov-modulated Lévy models, utilizing Carr’s randomization procedure. While our paper shares certain similarities with their work, it distinguishes itself by offering a unique contribution—analytical comparisons of the optional exercise boundaries and option prices, contingent on the volatility regime.

After Boyarchenko and Levendorskii (2009), Huang et al. (2011) introduce iterative procedures to value American options under regime switching. Zhang et al. (2014) exploit a penalty method to solve a system of complementarity problems arising from pricing American options. Yousuf et al. (2015) develop a second-order method based on an exponential time-differencing approach for solving American options under multi-state regime switching. Chiarella et al. (2016) solve the American option pricing problem under regime-switching by using the method-of-lines scheme. Lu and Putri (2020) use the Laplace transform method to solve the system of the partial differential equations for American option pricing. Including these papers, there are numerous papers dealing with the valuation of American options under regime switching. However, most of them do not show analytical comparisons based on the regimes covered in this paper.

The paper is organized as follows. Section 2 explains the model with regime-switching volatility. Section 3 studies a perpetual American put. Section 4 investigates pricing of an American put with finite expiry by randomization of time horizon and Section 5 discusses its implementation and shows numerical examples. Finally, Section 6 concludes. All the proofs are contained in the Appendix.

2. Model

In this section, we explain a model of a financial market. In the model, all activity occurs on a filtered probability space $(Ω, F, {F_{t}}_{t \geq 0}, P)$ ⁠, where ${F_{t}}$ is the augmentation of the filtration generated by a (one-dimensional) Brownian motion B and two Poisson processes φ_H(⋅) and φ_L(⋅) on $[0, \infty)$ ⁠, where the two Poisson processes are independent of each other and independent of the Brownian motion B.

In the financial market, there are two (underlying) assets. The first asset (the bond) is a money market account and is instantly risk-free. We assume that the return on the riskless asset is constant and equal to r. The second asset (the stock) is risky and American options are written on this asset. The price S_t of the stock evolves according to the following equation

\frac{d S (t)}{S (t)} = μ (t) d t + σ (t) d B (t), S (0) = S_{0},

where μ(t), σ(t) are ${F_{t}}$ -predictable processes. The stock pays no cash dividends ^[2].

The volatility of the second asset takes two values σ_H > σ_L > 0; namely, there are two regimes, say, “high volatility regime” H and “low volatility regime” L and volatility σ_H in regime H is larger than volatility σ_L in regime L, i.e. σ_H > σ_L > 0. This assumption implies, in particular, $σ (0) = σ_{0} = σ_{i_{0}}$ for some i₀ ∈ {H, L}, that is, the initial volatility regime is i at time 0. The volatility changes according to a Markov regime-switching process; regime i switches to regime j(≠ i) at the next jump time of Poisson process φ_i for i, j ∈ {H, L}, i.e. suppose at time t, the market is in regime i and φ_i(t) = n (n = 0, 1, 2, ⋯) then the volatility regime changes to j ≠ i at time inf{s: s ≥ t φ_i(s) = n + 1}. Under P, the intensity of the Poisson process φ_i is assumed to be ${\tilde{λ}}_{i} (t)$ at time t.

We assume that there is no market frictions, i.e. no taxes, no transaction costs and no short-sale constraints. Portfolio strategies, consumption processes and self-financing strategies with an initial wealth are defined as in a standard financial model (see, e.g. Chapter 1, Karatzas and Shreve, 1998).

We employ a standard assumption in modern finance, that is, there is no arbitrage opportunity in the financial market. Then, under a suitable regularity condition, there exists a probability measure Q such that for any financial contract $\tilde{S} (t)$ is a Q-martingale, where $\tilde{S} (t)$ is the time-t market price of the financial contract (see, e.g. Harrison and Pliska, 1981) ^[3]. The martingale measure may not be unique if the market is incomplete (see Chapter 5, Karatzas and Shreve, 1998). We will not assume the completeness of the financial market in this paper. However, we will assume the existence of a unique martingale measure that is used for pricing of financial contracts by market participants.

We will assume that the following is valid throughout the paper:

Standing assumption. There is a unique probability measure under which all the price of every financial contract is a martingale. The measure is used for pricing of all securities that will be introduced to the market. Under Q, the intensity of the Poisson process φ_i is λ_i for i ∈ {H, L}.

We will call Q the risk-neutral probability measure and the stock price process S(t) satisfies

\frac{d S (t)}{S (t)} = r d t + σ_{i} d \tilde{B} (t), in regime i,

where $\tilde{B}$ is a standard Brownian motion under Q. Note that a jump occurs only to the volatility of the stock, not to its price.

We now consider an American put option with maturity at T < ∞ and strike price K written on the stock. Suppose that the stock price is equal to S and the volatility regime is i ∈ {H, L} at time 0, then we define

P_{i} (T - t, S) \equiv ess \sup_{t \leq τ \leq T} E^{Q} [e^{- r (τ - t)} {(K - S_{τ})}^{+} | F_{t}] (ω), for ω \in {S (t) = S} \cap {σ (t) = σ_{i}},

where τ is an ${F_{s}}_{t \leq s \leq T}$ -stopping time. By the Markov property of the Brownian motion and Poisson processes, P_i(T − t, S) is uniquely determined. P_i(T − t, S) will be called the fair price of the American put at time t when the stock price and volatility regime at t are S and i, respectively.

The next proposition shows that our definition of the fair price is equal to the smallest wealth to hedge a short position of the American put if the market is complete. In order to state the proposition, we will consider (possibly non-existent) financial assets Θ_i for i = H or L whose price is perfectly correlated with the Poisson process φ_i. The financial market is complete if both Θ_H and Θ_L are tradable, i.e. the price processes exist and investors can form their portfolios using the assets.

Proposition 2.1.

Suppose that the financial market is complete, i.e. for j = H and L the financial contract Θ_j is tradable in the market. Suppose also that S₀ = S, σ(0) = σ_i(i ∈ {H, L}) Then, P_i(t, S) is equal to the smallest initial wealth required such that there exists a self-financing strategy whose wealth process X(u) satisfies

(i).
X(u) ≥ (K − S(u))⁺ for all 0 ≤ u < T.
(ii).
X(T) = (K − S(T))⁺ .

Now we state a lemma that will be useful to define the optimal exercise boundaries.

Lemma 2.2.

For i ∈ {H, L}, P_i(T, S) is Lipschitz continuous in S, uniformly continuous in T, nondecreasing in T, nonincreasing in S and convex in S.

Let

\begin{matrix} C_{i} \equiv {(u, S) \in {[0, \infty)}^{2} : P_{i} (u, S) > {(K - S)}^{+}} \\ S_{i} \equiv {(u, S) \in {[0, \infty)}^{2} \times [0, T) : P_{i} (u, S) = {(K - S)}^{+}}, \end{matrix}

for i ∈ {H, L}. By Lemma 2.2 $S_{i}$ is closed and $C_{i}$ is open and $S_{i, u} \equiv {S : (u, S) \in S_{i}}$ and $C_{i, u} \equiv {S : (u, S) \in C_{i}}$ are connected for every u ∈ [0, ∞)

For i ∈ {H, L}, let us define the optimal exercise boundary ${\underline{S}}_{i}$ ⁠: for $u \in [0, \infty)$

{\underline{S}}_{i} (u) \equiv \sup {S : S \in S_{i, u}} .

Then, the optimal stopping timeτ* for exercising the American put can be characterized as follows:

τ^{*} = \min_{i \in {H, L}} {\inf_{0 \leq u \leq T} {u : σ (u) = σ_{i} & S_{u} \leq {\underline{S}}_{i} (u)}} .

The following theorem states that the American put option price is the sum of the price of a European put option and the premium for the early exercise privilege. The results have been obtained by Kim (1990), Jacka (1991) and Carr et al. (1992) in the constant volatility case. For a measurable set, $A \in F_{t}$ ⁠, let χ_A denotes its characteristic function, i.e. χ_A(x) = 1, if x ∈ A and χ_A(x) = 0 if x ∈ Ω − A

Theorem 2.3.

For i ∈ {H, L}, let p_i(T, S) be the fair price of a European option with a strike price equal to K, time to maturity equal to T and the underlying state S at time T, when the current regime is i. Then,

P_{i} (T, S) = p_{i} (T, S) + \int_{0}^{T} e^{- r t} r K χ_{{S (t) < {\underline{S}}_{H} & σ (t) = σ_{H}} \cup {S (t) < {\underline{S}}_{L} & σ (t) = σ_{L}}} d t

In the constant volatility case, Mckean (1965) and Van Moerbeke (1976) have shown that the fair value of the American put is a solution to a free boundary value problem. We establish a similar result for the case with regime-switching volatility in the following. Let us define a differential operator $L$ for a pair (f_H(u, S), f_L(u, S)) of functions which are C² on a subset of $R^{2}$ ⁠: for i ∈ {H, L}

L f_{i} \equiv \frac{σ_{i}^{2}}{2} S^{2} f_{i, S S} + r S f_{i, S} - r f_{i} + λ_{i} (f_{j} - f_{i}) - f_{i, u},

where j ∈ {H, L}, j ≠ i. Here, subscripts denote partial derivatives, i.e. $f_{i, S S} = \frac{\partial^{2} f_{i}}{\partial S^{2}}, f_{i, S} = \frac{\partial f_{i}}{\partial S}$ and $f_{i, u} = \frac{\partial f_{i}}{\partial u}$ ⁠.

Then, we consider the following free boundary value problem.

Problem 2.4.

Find a pair (f_H(u, S), f_L(u, S)) of continuous functions defined on ${[0, \infty)}^{2}$ and a pair $(D_{H} (u), D_{L} (u))$ of boundaries defined on [0, T] such that

L f_{i} (u, S) = 0 for S > D_{i} (u), i \in {H, L}

subject to the following six boundary conditions

\lim_{S ↑ \infty} \max_{0 \leq u \leq T} | f_{i} (u, S) | = 0, \forall T < \infty \lim_{S ↓ D_{i} (u)} f_{i} (u, S) = K - D_{i} (u), \lim_{S ↓ D_{i} (u)} f_{i, S} (u, S) = - 1,

and four terminal conditions

D_{i} (T) = K, f_{i} (T, S) = {(K - S)}^{+}, for S > 0 .

Furthermore,

f_{i} (u, S) = K - S, for S \leq D_{i} (u), u \in [0, T],

and

f_{i} (u, S) \geq {(K - S)}^{+}, for u \in [0, T], S > 0 .

Proposition 2.5.

(P_H(u, S), P_L(u, S)) and $({\underline{S}}_{H} (u), {\underline{S}}_{L} (u))$ is the unique solution to Problem 2.4.

3. Perpetual American puts

In this section, we consider the case where the option’s lifetime is infinite, that is, T = ∞. McKean (1965) and Merton (1973) derived an analytic form for the option price for the case where the volatility is constant. In this section, we derive an analytic form for the value of a perpetual American put as in Guo and Zhang (2004). Furthermore, we provide a result concerning the existence and uniqueness of a solution to the free boundary value problem (Theorem 3.5) and show that the price of the perpetual put is always higher in the high-volatility regime than in the low-volatility regime and the exercise boundary is lower in the high-volatility regime than in the low-volatility regime (Theorems 3.7 and 3.8).

We will use the notation P_i(S) for the value of American put and ${\underline{S}}_{i}$ for the free boundary in regime i ∈ {H, L}, omitting time variable t. The notation is justified because the value and optimal exercise boundary of the perpetual American put option are time-homogeneous, i.e. independent of time t, given the volatility regime and stock price.

Now we consider the following free boundary value problem that will characterize the value of a perpetual American put option.

Problem 3.1.

Find a pair (f_H(S), f_L(S)) of continuous functions defined on (0, ∞) and a pair $(D_{H}, D_{L})$ of positive real numbers such that

L^{P} f_{i} (S) = 0 for S > D_{i}, i \in {H, L}

(3.1)

where

L^{P} f_{i} \equiv \frac{σ_{i}^{2}}{2} S^{2} f_{i, S S} + r S f_{i, S} - r f_{i} + λ_{i} (f_{j} - f_{i})

withj ≠ i, subject to the following six boundary conditions

\lim_{S ↑ \infty} | f_{i} (S) | = 0, \lim_{S ↓ D_{i}} f_{i} (S) = K - D_{i}, \lim_{S ↓ D_{i}} f_{i, S} (S) = - 1,

(3.2)

Furthermore,

f_{i} (S) = K - S, for S \leq D_{i},

and

f_{i} (S) \geq {(K - S)}^{+}, for S > 0 .

Proposition 3.2.

(P_H(S), P_L(S)) and $({\underline{S}}_{H}, {\underline{S}}_{L})$ is the unique solution to Problem 3.1.

We will later show that ${\underline{S}}_{L} > {\underline{S}}_{H}$ ⁠. However, we still do not know which of ${\underline{S}}_{L}$ and ${\underline{S}}_{H}$ is larger and will temporarily let ${\underline{S}}_{M} \equiv \max {{\underline{S}}_{H}, {\underline{S}}_{L}}$ and ${\underline{S}}_{m} \equiv \min {{\underline{S}}_{H}, {\underline{S}}_{L}}$ ⁠. Namely, M and m will denote H or L such that ${\underline{S}}_{M} \geq {\underline{S}}_{m}$ (When ${\underline{S}}_{M} = {\underline{S}}_{m}$ ⁠, we will choose notation such that M = L and m = H).

In the following, we will derive the unique solution (P_H(S), P_L(S)) and $({\underline{S}}_{H}, {\underline{S}}_{L})$ to Problem 3.1 in explicit form. An alternative derivation can be found in Guo and Zhang (2004).

Lemma 3.3.

There exist two real negative numbers n₁ and n₂(n₁ > n₂) such that a general solution to equations (3.1) for $S > {\underline{S}}_{M}$ subject to boundary conditions $\lim_{S ↑ \infty} P_{i} (S) = 0$ for i, j ∈ {H, L} are given by

P_{H} (S) = C_{1} S^{n_{1}} + C_{2} S^{n_{2}}, P_{L} (S) = C_{1} ξ_{1} S^{n_{1}} + C_{2} ξ_{2} S^{n_{2}} .

(3.3)

for some

C_{1}, C_{2} \in R

⁠. Furthermore,n₁, n₂, ξ₁andξ₂satisfy

\begin{matrix} n_{2} < - 2 r / σ_{L}^{2} < n_{1} < - 2 r / σ_{H}^{2}, \\ 1 > ξ_{1} > 0 > ξ_{2} . \end{matrix}

(3.4)

Lemma 3.4.

Suppose that ${\underline{S}}_{m} < S < {\underline{S}}_{M}$ ⁠. Then, the value of the perpetual American put P_m(S) in regime m takes the following form:

P_{m} (S) = D_{1} S^{m_{1}} + D_{2} S^{m_{2}} - S + \frac{K λ_{m}}{r + λ_{m}}

(3.5)

for some constantsD₁andD₂, where

m_{1} = \frac{σ_{m}^{2} - 2 r + \sqrt{a_{m}}}{2 σ_{m}^{2}}, m_{2} = \frac{σ_{m}^{2} - 2 r - \sqrt{a_{m}}}{2 σ_{m}^{2}}, a_{m} = {(2 r - σ_{m}^{2})}^{2} + 8 σ_{m}^{2} (r + λ_{m}) > 0 .

Furthermore,m₁, m₂satisfy

m_{1} > 1, m_{2} < 0 .

(3.6)

Now we proceed to determine constants ${\underline{S}}_{m}$ ⁠, ${\underline{S}}_{M}$ ⁠, C₁, C₂, D₁, D₂. Let us define

\begin{matrix} f (z, s) \equiv \frac{r + λ_{m}}{r z K} (\frac{(z - n_{1}) (n_{2} (K - s) + s)}{(n_{2} - n_{1}) ξ_{1}} + \frac{(z - n_{2}) (n_{1} (K - s) + s)}{(n_{1} - n_{2}) ξ_{2}} + (z - 1) s - \frac{z λ_{m} K}{r + λ_{m}}) . \end{matrix}

The following result shows how to determine the solution to the free boundary problem.

Proposition 3.5.

${\underline{S}}_{M}, {\underline{S}}_{m}, C_{1}, C_{2}$ satisfy the followings:

{[f (m_{1}, {\underline{S}}_{M})]}^{m_{1}} = {[f (m_{2}, {\underline{S}}_{M})]}^{m_{2}},

({\frac{{\underline{S}}_{M}}{{\underline{S}}_{m}})}^{m_{2}} = f (m_{1}, {\underline{S}}_{M}), ({\frac{{\underline{S}}_{M}}{{\underline{S}}_{m}})}^{m_{1}} = f (m_{2}, {\underline{S}}_{M}),

(3.7)

and

C_{1} = \frac{n_{2} (K - {\underline{S}}_{M}) + {\underline{S}}_{M}}{(n_{2} - n_{1}) ξ_{1} {\underline{S}}_{M}^{n_{1}}}, C_{2} = \frac{n_{1} (K - {\underline{S}}_{M}) + {\underline{S}}_{M}}{(n_{1} - n_{2}) ξ_{2} {\underline{S}}_{M}^{n_{2}}},

(3.8)

D_{1} = \frac{r m_{2} K}{(r + λ_{m}) (m_{2} - m_{1}) {\underline{S}}_{m}^{m_{1}}}, D_{2} = \frac{r m_{1} K}{(r + λ_{m}) (m_{1} - m_{2}) {\underline{S}}_{m}^{m_{2}}} .

(3.9)

Corollary 3.6.

(i) Both D₁ and D₂ are positive and (ii) C₁ ≠ 0 and C₂ ≠ 0.

We now show in the following theorems that the price of the perpetual put is always higher in the high-volatility regime than in the low-volatility regime and the optimal exercise boundary is lower in the high-volatility regime than in the low-volatility regime.

Theorem 3.7.

${\underline{S}}_{L} > {\underline{S}}_{H}$ , i.e. m = H and M = L.

Theorem 3.8.

P_H(S) ≥ P_L(S) for all S > 0 and the inequality is strict for $S > {\underline{S}}_{H}$ ⁠.

In Figure 1, we draw a graph of the value of a perpetual American put option where the parameters are given as r = 0.1, K = 1, T = 1, σ_H = 0.4, σ_L = 0.2, λ_H = 1.0 and λ_L = 0.5. The intensity parameters in the figure are chosen such as the average time it takes a regime to change from high volatility to low volatility is a year, while the average time it takes a regime to change from low volatility to high volatility is half a year. The higher intensity λ_H in the high volatility regime H is supported by empirical studies, for example, Ang and Bekaert (2002). The result shown in the figure is consistent with Theorems 3.7 and 3.8.

Figure 1

View large Download slide

Value of a perpetual American put option

These results are also consistent with our intuition: High volatility yields a higher time value, so the price P_H is higher and the corresponding optimal exercise boundary is lower. In the continuation region (the upper region of the optimal exercise boundary), a lower optimal exercise boundary gives the put holder a higher time value because the current stock price is above it.

4. American puts with finite expiry

In this section, we will study the valuation of an American put which has a finite expiration time T. In the case where volatility is constant Carr (1998) has studied the same problem by considering an American put with a random expiration time which is the sum of N independent random times that are identically distributed according to an exponential distribution with a mean equal to T/N. The random time is distributed according to an Erlang distribution with a mean equal to T and a variance equal to T/N, and therefore, the random-time approaches to T in probability as N → ∞. He has suggested that the value of an American put with finite time T to expiration can be approximated by the value of an American put with the previously explained random time if N is sufficiently large. We will apply his technique to study the value of an American put with finite expiry.

4.1 Exponentially distributed time

In this section, we will consider an American put with a random time τ until expiration, where

P [τ \in d t] = β e^{- β t} d t

with β ≡ 1/T. Namely, there is an independent Poisson process with an intensity β such that the option expires at the first jump time of the process ^[4].

Note that the expected time to expiration of the American put is β. Analysis of the American put with a random expiration time is very similar to that of the perpetual option in the previous section. As in the case of the perpetual option, the value and early exercise boundary of the American put option are time-homogeneous, i.e. independent of time t, given the volatility regime and stock price. Therefore, the equations describing the fair value of the American put are ordinary differential equations without having time derivatives. We will use notation $P_{i}^{1} (S)$ and ${\underline{S}}_{i}^{1}$ to denote the value and early exercise boundary of the American put, respectively, in regime i ∈ {H, L}.

As in Carr (1998), the free boundary problem for the fair price of the American put takes the following form: for i, j ∈ {H, L}, i ≠ j,

\{\begin{matrix} \frac{1}{2} σ_{i}^{2} S^{2} P_{i, S S}^{1} (S) + r S P_{i, S}^{1} (S) - (r + λ_{i} + β) P_{i}^{1} (S) + λ_{i} P_{j}^{1} (S) = 0 & if S > K \\ \frac{1}{2} σ_{i}^{2} S^{2} P_{i, S S}^{1} (S) + r S P_{i, S}^{1} (S) - (r + λ_{i} + β) P_{i}^{1} (S) + λ_{i} P_{j}^{1} (S) = - β (K - S) & if {\underline{S}}_{i}^{1} < S < K \\ P_{i}^{1} (S) = K - S & if S < {\underline{S}}_{i}^{1}, \end{matrix}

with boundary conditions

\lim_{S ↑ \infty} P_{i}^{1} (S) = 0, \lim_{S ↓ {\underline{S}}_{i}^{1}} P_{i}^{1} (S) = K - {\underline{S}}_{i}^{1}, \lim_{S ↓ {\underline{S}}_{i}^{1}} P_{i, S}^{1} (S) = - 1 .

As in the previous section, let ${\underline{S}}_{M}^{1} \equiv \max {{\underline{S}}_{H}^{1}, {\underline{S}}_{L}^{1}}$ and ${\underline{S}}_{m}^{1} \equiv \min {{\underline{S}}_{H}^{1}, {\underline{S}}_{L}^{1}}$ ⁠.

We will make the following assumption in this section.

Assumption I.

There exist four distinct real roots n₁, n₂, n₃, n₄ of the algebraic equation I_F(n) = 0, where

I_{F} (n) = (\frac{1}{2} σ_{H}^{2} n^{2} + (r - \frac{1}{2} σ_{H}^{2}) n - (r + λ_{H} + β)) (\frac{1}{2} σ_{L}^{2} n^{2} + (r - \frac{1}{2} σ_{L}^{2}) n - (r + λ_{L} + β)) - λ_{H} λ_{L},

and among them the two roots n₁, n₂(n₁ > n₂) are negative and others n₃, n₄(n₃ > n₄) are greater than 1.

Note the fact that I_F(0) > 0 and I_F(1) > 0. Thus, if the roots of equation I_F(n) = 0 exist in one of the open intervals (−∞, 0), (0, 1) and (1, ∞), then there are exactly two or four. Assumption I tells us that we assume there exist exactly two roots in (−∞, 0) and the other two roots in (1, ∞), respectively.

Proposition 4.1.

Suppose that Assumption I is valid and let ξ₁, ξ₂, ξ₃ and ξ₄ be unique solutions to the following equation:

\{\begin{matrix} \frac{1}{2} σ_{H}^{2} n (n - 1) + r n - r - λ_{H} - β + λ_{H} ξ = 0 \\ ξ (\frac{1}{2} σ_{L}^{2} n (n - 1) + r n - r - λ_{L} - β) + λ_{L} = 0 . \end{matrix}

Then the fair value of the American put with random expiration time τ takes the following form.

(i) In the region S > K, for some constants $C_{1}^{1}, C_{2}^{1}$

P_{m}^{1} (S) = C_{1}^{1} S^{n_{1}} + C_{2}^{1} S^{n_{2}}, P_{M}^{1} (S) = C_{1}^{1} ξ_{1} S^{n_{1}} + C_{2}^{1} ξ_{2} S^{n_{2}} .

(ii) In the region ${\underline{S}}_{M}^{1} < S < K$ ⁠, for some constants ${\tilde{C}}_{1}^{1}, {\tilde{C}}_{2}^{1}, {\tilde{C}}_{3}^{1}$ and ${\tilde{C}}_{4}^{1}$ ⁠,

\{\begin{matrix} P_{m}^{1} (S) = {\tilde{C}}_{1}^{1} S^{n_{1}} + {\tilde{C}}_{2}^{1} S^{n_{2}} + {\tilde{C}}_{3}^{1} S^{n_{3}} + {\tilde{C}}_{4}^{1} S^{n_{4}} - S + \frac{β K}{r + β}, \\ P_{M}^{1} (S) = {\tilde{C}}_{1}^{1} ξ_{1} S^{n_{1}} + {\tilde{C}}_{2}^{1} ξ_{2} S^{n_{2}} + {\tilde{C}}_{3}^{1} ξ_{3} S^{n_{3}} + {\tilde{C}}_{4}^{1} ξ_{4} S^{n_{4}} - S + \frac{β K}{r + β} . \end{matrix}

(iii) In the region ${\underline{S}}_{m}^{1} < S < {\underline{S}}_{M}^{1}$ ⁠, for some constants $D_{1}^{1}$ and $D_{2}^{1}$ ⁠,

P_{m}^{1} (S) = D_{1}^{1} S^{m_{1}} + D_{2}^{1} S^{m_{2}} - S + \frac{K (λ_{m} + β)}{r + λ_{m} + β},

where

m_{1} = \frac{σ_{m}^{2} - 2 r + \sqrt{a_{m}}}{2 σ_{m}^{2}}

and

m_{2} = \frac{σ_{m}^{2} - 2 r - \sqrt{a_{m}}}{2 σ_{m}^{2}}, a_{m} = {(2 r - σ_{m}^{2})}^{2} + 8 σ_{m}^{2} (r + λ_{m} + β) > 0

and

P_{M}^{1} (S) = K - S .

There are 10 equations that can be derived from the smoothness of $P_{M}^{1} (\cdot)$ and $P_{m}^{1} (\cdot)$ ⁠. And by some heavy calculation, we can reduce the 10 equations to two similar to the case of a perpetual put ^[5].

Theorem 4.2.

Suppose that Assumption I is valid. Then, ${\underline{S}}_{L}^{1} > {\underline{S}}_{H}^{1}$ and $P_{H}^{1} (S) \geq P_{L}^{1} (S)$ for all S > 0. The inequality is strict for $S > {\underline{S}}_{H}^{1}$ ⁠.

4.2 Erlang-distributed time

In this section, we will consider an American put with a random expiration time τ which is the sum of N independent random-times which are identically and exponentially distributed with a mean equal to T/N. That is, the American put will expire after the N-th jump time of an independent Poisson process with intensity β_N = N/T. Then, τ is distributed according to an Erlang distribution, that is,

P [τ \in d t] = \frac{{(β_{N})}^{N}}{(N - 1)!} t^{N - 1} e^{- β_{N} t} d t .

Let us consider, N put options with the same strike price K. The first option is the put with an exponentially distributed time to expiration with β = β_N explained in the previous section. The k-th option for 1 < k ≤ N is a put which can be exercised for (K − S)⁺ at any time up to and including the k-th jump time of a Poisson Process with intensity β_N, and if the option is not exercised after the k-th jump time of the Poisson Process it becomes the (k − 1)-th put after the jump. Let $P_{i}^{k} (S) (k = 2,3, \dots, N)$ denote the fair value of the k-th American put in regime i ∈ {H, L}. We conjecture that the N-th put value $P_{i}^{N} (S)$ converges to the exact value of the American put P_i(0, S) as N goes to infinity, as Carr (1998) did. Numerical results that will be shown later support this conjecture, but the proof of the conjecture is still unknown.

If we follow the argument in Carr (1998), we can easily derive the following equation for $P_{i}^{k} (\cdot)$ ⁠:

\begin{matrix} \frac{1}{2} σ_{i}^{2} S^{2} P_{i, S S}^{k} (S) + r S P_{i, S}^{k} (S) - r P_{i}^{k} (S) + λ_{i} (P_{j}^{k} (S) & - P_{i}^{k} (S)) = β_{N} (P_{i}^{k} (S) - P_{i}^{k - 1} (S)), \\ i \neq j, j \in {H, L}, for S > {\underline{S}}_{i}^{k} . \end{matrix}

Let ${\underline{S}}_{M}^{k} \equiv \max {{\underline{S}}_{H}^{k}, {\underline{S}}_{L}^{k}}$ and ${\underline{S}}_{m}^{k} \equiv \min {{\underline{S}}_{H}^{k}, {\underline{S}}_{L}^{k}}$ ⁠. We use the same notation n_l’s and ξ_l’s(l = 1, 2, 3, 4) as in the previous section with β = β_N.

The problem is to find $P_{i}^{k} (S)$ satisfying

\frac{1}{2} σ_{i}^{2} S^{2} P_{i, S S}^{k} (S) + r S P_{i, S}^{k} (S) - r P_{i}^{k} (S) + λ_{i} (P_{j}^{k} (S) - P_{i}^{k} (S)) = β_{N} (P_{i}^{k} (S) - P_{i}^{k - 1} (S)), i \neq j

and six boundary conditions

\lim_{S ↑ \infty} P_{i}^{k} (S) = 0, \lim_{S ↓ {\underline{S}}_{i}^{k}} P_{i}^{k} (S) = K - {\underline{S}}_{i}^{k}, \lim_{S ↓ {\underline{S}}_{i}^{k}} P_{i, S}^{k} (S) = - 1 .

Theorem 4.3.

Suppose that Assumption I is valid and the determinant of Φ(1), defined in (6.13), is not zero. Then $P_{m}^{k} (S)$ and $P_{M}^{k} (S)$ for k = 1, 2, …, N take the following form:

(i) In the region $S > {\underline{S}}_{M}^{k}$ ⁠, for two constants $C_{h}^{k}$ ’s(h = 1, 2)

\{\begin{cases} P_{m}^{k} (S) = \sum_{h = 1}^{4} (C_{h}^{k} + u_{h}^{k} (S)) S^{n_{h}} \\ P_{M}^{k} (S) = \sum_{h = 1}^{4} ξ_{h} (C_{h}^{k} + u_{h}^{k} (S)) S^{n_{h}} \end{cases}

(4.1)

where

u_{h}^{k} (\cdot)

’s are defined by (6.15) in Appendix and

C_{h}^{k} = - \lim_{S \to \infty} u_{h}^{k} (S) for h = 3,4 .

(4.2)

(ii) In the region ${\underline{S}}_{m}^{k} < S < {\underline{S}}_{M}^{k}$ ⁠, for some constants $D_{1}^{k}$ and $D_{2}^{k}$ and m₁, m₂ defined in the previous section,

\begin{matrix} P_{m}^{k} (S) = & D_{1}^{k} S^{m_{1}} + D_{2}^{k} S^{m_{2}} \\ - S^{m_{1}} \int_{{\underline{S}}_{m}^{k}}^{S} \frac{2 (- λ_{m} K + λ_{m} t - β P_{m}^{k - 1} (t))}{(m_{2} - m_{1}) σ_{m}^{2} t^{m_{1} + 1}} d t + S^{m_{2}} \int_{{\underline{S}}_{m}^{k}}^{S} \frac{2 (- λ_{m} K + λ_{m} t - β P_{m}^{k - 1} (t))}{(m_{2} - m_{1}) σ_{m}^{2} t^{m_{2} + 1}} d t, \end{matrix}

and

P_{M}^{k} (S) = K - S .

Theorem 4.4.

Suppose that Assumption I is valid. ${\underline{S}}_{L}^{k} > {\underline{S}}_{H}^{k}$ and $P_{H}^{k} (S) \geq P_{L}^{k} (S)$ for all S > 0. The inequality is strict for $S > {\underline{S}}_{H}^{k}$ ⁠.

5. Implementation and numerical results

In this section, we will show numerical results. First, we will propose the Richardson extrapolation together with the values of American puts with a random expiration time studied in the previous section as an approximation to the American option with finite expiry. We will compare the solutions obtained from this method with the solutions obtained by pentanomial trees in Bollen (1998) for a certain set of parameter values. Second, we will illustrate behavior of American put values in a regime-switching model with a parameter change by using numerical solutions obtained from the Richardson extrapolation method.

5.1 Richardson extrapolation with randomized put values

We first explain the Richardson extrapolation method. As in the previous section, let P_i(0, S) denote the exact American Put value with time to expiration equal to T < ∞ at time 0 when the stock price is equal to S and regime is i ∈ {H, L}. And let $P_{i}^{N} (S; N)$ the value of an American put with an Erlang-distributed random time to expiration with β = N/T. Then, as in Geske and Johnson (1984), Broadie and Detemple (1996) and Carr (1998), the Richardson extrapolation for M ≥ 2, ${\tilde{P}}_{i}^{M} (0, S)$ ⁠, is defined as

{\tilde{P}}_{i}^{M} (0, S) \equiv \sum_{N = 1}^{M} \frac{{(- 1)}^{M - N} N^{M}}{N! (M - N)!} P_{i}^{N} (S; N) .

For example, a three-point Richardson extrapolation is represented by

{\tilde{P}}_{i}^{3} (0, S) = \frac{1}{2} P_{i}^{1} (S; 1) - 4 P_{i}^{2} (S; 2) + \frac{9}{2} P_{i}^{3} (S; 3) .

We conjecture that ${\tilde{P}}_{i}^{M} (0, S)$ ⁠, for a sufficiently large M, is a good approximation to P_i(0, S), namely,

P_{i} (0, S) \approx {\tilde{P}}_{i}^{M} (0, S),

for a sufficiently large M.

Table 1 exhibits the value of an American put calculated by the pentanomial-tree method with 10³ time steps and those calculated by the method in this paper with three-point Richardson extrapolation. The default parameters are r = 0.1, K = 1, T = 1, σ_L = 0.2, λ_L = 0.5. It clearly shows that we can approximate the value of the American put quite accurately by the three-point Richardson extrapolation together with randomized put values (mostly errors are less than 0.0003).

Table 1

Comparison of the value of an American put option for a set of parameter values with default parameters r = 0.1, K = 1, T = 1, σ_L = 0.2, λ_L = 0.5

Parameters			Pentanomial tree		3-Point Richardson
S	σ_H	λ_H	P_H(0, S)	P_L(0, S)	${\tilde{P}}_{H}^{3} (0, S)$	${\tilde{P}}_{L}^{3} (0, S)$
0.9	0.4	1.0	0.1483	0.1106	0.1483	0.1106
0.9	0.4	2.0	0.1390	0.1093	0.1393	0.1094
0.9	0.5	1.0	0.1738	0.1150	0.1737	0.1149
0.9	0.5	2.0	0.1594	0.1128	0.1597	0.1127
1.0	0.4	1.0	0.1015	0.0594	0.1014	0.0592
1.0	0.4	2.0	0.0904	0.0574	0.0905	0.0572
1.0	0.5	1.0	0.1293	0.0660	0.1292	0.0658
1.0	0.5	2.0	0.1128	0.0629	0.1131	0.0626

Parameters			Pentanomial tree		3-Point Richardson
S	σ_H	λ_H	P_H(0, S)	P_L(0, S)	${\tilde{P}}_{H}^{3} (0, S)$	${\tilde{P}}_{L}^{3} (0, S)$
0.9	0.4	1.0	0.1483	0.1106	0.1483	0.1106
0.9	0.4	2.0	0.1390	0.1093	0.1393	0.1094
0.9	0.5	1.0	0.1738	0.1150	0.1737	0.1149
0.9	0.5	2.0	0.1594	0.1128	0.1597	0.1127
1.0	0.4	1.0	0.1015	0.0594	0.1014	0.0592
1.0	0.4	2.0	0.0904	0.0574	0.0905	0.0572
1.0	0.5	1.0	0.1293	0.0660	0.1292	0.0658
1.0	0.5	2.0	0.1128	0.0629	0.1131	0.0626

Note(s): The values in Columns 4 and 5 are obtained by the pentanomial-tree method in Bollen (1998) with 10³ time steps and those in Columns 6 and 7 are obtained by our accelerated method with a three-point Richardson extrapolation

Source(s): Authors' own work

5.2 Characteristics of American put values

Figure 2 shows the value of an American put option as a function of volatility σ_H of the stock in state H for the following parameter values: r = 0.1, K = 1, T = 1, σ_L = 0.2, λ_H = 2.0, λ_L = 0.5 and S = 0.9. It shows the values in both the high-volatility state and low-volatility state increase as the high-state volatility σ_H increases.

Figure 2

View large Download slide

Value of an American put option as a function of volatility σ_H of the stock in state H

Figure 3 shows the value of an American put option as a function of volatility σ_L of the stock in state L for the following parameter values: r = 0.1, K = 1, T = 1, σ_H = 0.4, λ_H = 2.0, λ_L = 0.5 and S = 0.9. It shows the values in both the high-volatility state and low-volatility state increase as the low-state volatility σ_L increases. Put prices tend to increase as either the high-state volatility or the low-state volatility increases.

Figure 3

View large Download slide

Value of an American put option as a function of volatility σ_L of the stock in state L

The results in Figures 2 and 3 are because there is a possibility that the future regime may differ from the current regime through regime shifts, thus the option price in a regime might be highly influenced by the volatility of the opposite regime. This implies, if a period incorporating regime shifts is considered as the lifetime of an American option, practitioners should assess the market value of options taking into account the volatility of a regime different from the present one.

Figure 4 shows the value of an American put option as a function of intensity parameter λ_H for the following parameter values r = 0.1, K = 1, T = 1, σ_H = 0.4, σ_L = 0.2, λ_L = 0.5 and S = 0.9. It shows the values in both the high-volatility state and low-volatility state decline as the intensity λ_H increases. Namely, as it is more likely for the high-volatility regime to change to the low-volatility regime, the option values in both the high-volatility state and low-volatility state tend to decline. Figure 5 shows the value of an American put option as a function of intensity parameter λ_L for the following parameter values r = 0.1, K = 1, T = 1, σ_H = 0.4, σ_L = 0.2, λ_H = 2.0 and S = 0.9. It shows the values in both the high-volatility state and low-volatility state increase as the intensity λ_L increases.

Figure 4

View large Download slide

Value of an American put option as a function of intensity parameter λ_H

Figure 5

View large Download slide

Value of an American put option as a function of intensity parameter λ_L

We conjecture the above comparative statics results illustrated by numerical examples are generally valid and leave their analytic proof as an open question.

6. Conclusion

We introduce a novel approach to find the value of American puts with regime-switching volatility. We demonstrate that the value of an American put can be represented as the combination of the corresponding European put price and the premium attributed to the early exercise privilege. We also provide rigorous proof that the American put value constitutes a unique solution to the free boundary value problem. Furthermore, we show that perpetual American put options can exhibit distinct characteristics depending on the prevailing volatility regime. Specifically, we establish that perpetual put options command a higher price during periods of high-volatility compared to low-volatility phases. Additionally, we reveal that the optimal exercise boundary is situated at a lower threshold in high-volatility regimes, as opposed to low-volatility regimes. These findings contribute valuable insights into the dynamics of option pricing in the context of fluctuating volatility regimes.

Moreover, we develop an analytical framework to describe American puts with an Erlang-distributed random-time horizon, which allows us to propose a numerical technique for approximating the value of American puts with finite expiry. We also show that a combined approach involving randomization and Richardson extrapolation can be a numerical algorithm for estimating American put prices with finite expiry.

This work was supported by the Ministry of Education of the Republic of Korea and the National Research Foundation of Korea (NRF- 2021S1A5A2A03065678).

Notes

1.

Yi (2008) shows the same result with ours.

2.

Analysis in this paper can be easily extended to the case where the stock pays cash dividends at a constant rate.

3.

The probability measure Q is equivalent to P on $(Ω, F_{T}^{B, φ_{H}, φ_{L}})$ for T < ∞, where $F_{T}^{B, φ_{H}, φ_{L}}$ is a σ-algebra generated by B(s), φH(s), φL(s) for s ∈ [0, T]. But in an infinite horizon P and Q are not equivalent. See Section 1.7, Karatzas and Shreve (1998).

4.

We need to enlarge the filtration ${F_{t}}$ to make the Poisson process progressively measurable. We will assume such a modification has been done to the filtration.

5.

The two equations are very lengthy and we do not report them in the paper. They are available from the authors upon request.

References

Ang

,

A.

and

Bekaert

,

G.

(

2002

), “

International asset allocation with regime shift

”,

Review of Financial Studies

, Vol.

15

No.

4

, pp.

1137

-

1187

, doi:

https://doi.org/10.1093/rfs/15.4.1137

.

Google Scholar

Crossref

Bates

,

D.S.

(

1996

), “

Jumps and stochastic volatility: exchange rate processes implicit in deutsche mark options

”,

Review of Financial Studies

, Vol.

9

No.

1

, pp.

69

-

107

, doi:

https://doi.org/10.1093/rfs/9.1.69

.

Google Scholar

Crossref

Black

,

F.

and

Scholes

,

M.

(

1973

), “

The pricing of options and corporate liabilities

”,

Journal of Political Economy

, Vol.

81

No.

3

, pp.

637

-

654

, doi:

https://doi.org/10.1086/260062

.

Google Scholar

Crossref

Bollen

,

N.P.B.

(

1998

), “

Valuing options in regime-switching models

”,

Journal of Derivatives

, Vol.

6

No.

1

, pp.

38

-

49

, doi:

https://doi.org/10.3905/jod.1998.408011

.

Google Scholar

Crossref

Bollen

,

P.B.

,

Gray

,

S.

and

Whaley

,

R.

(

2000

), “

Regime-switching in foreign exchange rates: evidence from currency option prices

”,

Journal of Econometrics

, Vol.

94

Nos

1-2

, pp.

239

-

276

, doi:

https://doi.org/10.1016/s0304-4076(99)00022-6

.

Google Scholar

Crossref

Boyarchenko

,

S.

and

Levendorskii

,

S.

(

2009

), “

American options in regime-switching models

”,

SIAM Journal on Control and Optimization

, Vol.

48

No.

3

, pp.

1353

-

1376

, doi:

https://doi.org/10.1137/070682897

.

Google Scholar

Crossref

Broadie

,

M.

and

Detemple

,

J.

(

1996

), “

American option valuation: new bounds, approximations, and a comparison of existing methods

”,

Review of Financial Studies

, Vol.

9

No.

4

, pp.

1211

-

1250

, doi:

https://doi.org/10.1093/rfs/9.4.1211

.

Google Scholar

Crossref

Campbell

,

S.

and

Li

,

C.

(

2002

), “

Option pricing in regime-switching environments

”,

Working Paper, Brown University

.

Google Scholar

Carr

,

P.

(

1998

), “

Randomization and the American put

”,

Review of Financial Studies

, Vol.

11

No.

3

, pp.

597

-

626

, doi:

https://doi.org/10.1093/rfs/11.3.597

.

Google Scholar

Crossref

Carr

,

P.

,

Jarrow

,

R.

and

Myneni

,

R.

(

1992

), “

Alternative characterizations of American put options

”,

Mathematical Finance

, Vol.

2

, pp.

87

-

106

, doi:

https://doi.org/10.1111/j.1467-9965.1992.tb00040.x

.

Google Scholar

Crossref

Chiarella

,

C.

,

Nikitopoulos

,

C.S.

and

Schlögl

,

E.

(

2016

), “

Pricing American options under regime switching using method of lines

”,

Working Paper, University of Technology Sydney

.

Google Scholar

Chourdakis

,

K.

(

2001

), “

Volatility persistence, regime switches, jumps and option pricing

”,

Working Paper, University of London

.

Google Scholar

Chourdakis

,

K.

and

Tzavalis

,

E.

(

2000

), “

Option pricing under discrete shifts in stock returns

”,

Working Paper, University of London

.

Google Scholar

Driffill

,

J.

,

Kenc

,

T.

and

Sola

,

M.

(

2002

), “

Merton-style option pricing under regime switching

”,

Working Paper, University of London

.

Google Scholar

Duan

,

J.C.

,

Popova

,

I.

and

Ritchken

,

P.

(

1999

), “

Option pricing under regime switching

”,

Working Paper, Hong Kong University of Science and Technology

.

Google Scholar

Duffie

,

D.

,

Pan

,

J.

and

Singleton

,

K.

(

2000

), “

Transform analysis and asset pricing for affine jump-diffusions

”,

Econometrica

, Vol.

68

No.

6

, pp.

1343

-

1376

, doi:

https://doi.org/10.1111/1468-0262.00164

.

Google Scholar

Crossref

Dybvig

,

P.

and

Ross

,

S.A.

(

2003

),

Arbitrage, State Prices and Portfolio Theory: Handbook of the Economics of Finance

,

Washington University in St. Louis

.

Google Scholar

Edwards

,

C.

(

2005

), “

Derivative pricing models with regime switching: a general approach

”,

Journal of Derivatives

, Vol.

13

No.

1

, pp.

41

-

47

, doi:

https://doi.org/10.3905/jod.2005.580516

.

Google Scholar

Crossref

Elliot

,

R.J.

(

1980

),

Stochastic Calculus and Applications

,

Springer-Verlag

,

New York

.

Google Scholar

Friedman

,

A.

(

1964

),

Partial Differential Equations of Parabolic Type

,

Prentice-Hall

,

London

.

Google Scholar

Geske

,

R.

and

Johnson

,

H.

(

1984

), “

The American put option valued analytically

”,

Journal of Finance

, Vol.

39

No.

5

, pp.

1511

-

1524

, doi:

https://doi.org/10.1111/j.1540-6261.1984.tb04921.x

.

Google Scholar

Crossref

Guo

,

S.

and

Zhang

,

Q.

(

2004

), “

Closed-form solutions for perpetual American put options with regime switching

”,

SIAM Journal on Applied Mathematics

, Vol.

64

No.

6

, pp.

2034

-

2049

, doi:

https://doi.org/10.1137/s0036139903426083

.

Google Scholar

Crossref

Harrison

,

J.M.

and

Kreps

,

D.

(

1979

), “

Martingales and arbitrage in multiperiod securities markets

”,

Journal of Economic Theory

, Vol.

20

No.

3

, pp.

381

-

408

, doi:

https://doi.org/10.1016/0022-0531(79)90043-7

.

Google Scholar

Crossref

Harrison

,

J.M.

and

Pliska

,

S.R.

(

1981

), “

Martingales and stochastic integrals in the theory of continuous trading

”,

Stochastic Processes and Applications

, Vol.

11

No.

3

, pp.

215

-

260

, doi:

https://doi.org/10.1016/0304-4149(81)90026-0

.

Google Scholar

Crossref

Heston

,

S.L.

(

1993

), “

A closed-form solution for options with stochastic volatility with applications to bond and currency options

”,

Review of Financial Studies

, Vol.

6

No.

2

, pp.

327

-

343

, doi:

https://doi.org/10.1093/rfs/6.2.327

.

Google Scholar

Crossref

Huang

,

Y.

,

Forsyth

,

P.A.

and

Labahn

,

G.

(

2011

), “

Methods for pricing American options under regime switching

”,

SIAM Journal on Scientific Computing

, Vol.

33

No.

5

, pp.

2144

-

2168

, doi:

https://doi.org/10.1137/110820920

.

Google Scholar

Crossref

Hull

,

J.

and

White

,

A.

(

1987

), “

The pricing of options on assets with stochastic volatilities

”,

Journal of Finance

, Vol.

42

No.

2

, pp.

281

-

300

, doi:

https://doi.org/10.2307/2328253

.

Google Scholar

Crossref

Jacka

,

S.D.

(

1991

), “

Optimal stopping and the American put

”,

Mathematical Finance

, Vol.

1

No.

2

, pp.

1

-

14

, doi:

https://doi.org/10.1111/j.1467-9965.1991.tb00007.x

.

Google Scholar

Crossref

Karatzas

,

I.

and

Shreve

,

S.E.

(

1998

),

Methods of Mathematical Finance

,

Springer-Verlag

,

New York

.

Google Scholar

Kim

,

I.J.

(

1990

), “

The analytic valuation of American options

”,

The Review of Financial Studies

, Vol.

3

No.

4

, pp.

547

-

572

, doi:

https://doi.org/10.1093/rfs/3.4.547

.

Google Scholar

Crossref

Lu

,

X.

and

Putri

,

E.

(

2020

),

Physica A

, Vol.

538

, 122968, doi:

https://doi.org/10.1016/j.physa.2019.122968

.

Google Scholar

McKean

,

H.P.

, Jr (

1965

), “

Appendix: a free boundary problem for the heat equation arising from a problem in mathematical economics

”,

Industrial Management Review

, Vol.

6

, pp.

32

-

39

.

Google Scholar

Merton

,

R.

(

1973

), “

The theory of rational option pricing

”,

Bell Journal of Economics and Management Science

, Vol.

4

No.

1

, pp.

141

-

183

, doi:

https://doi.org/10.2307/3003143

.

Google Scholar

Crossref

Naik

,

V.

(

1993

), “

Option valuation and hedging strategies with jumps in the volatility of asset returns

”,

Journal of Finance

, Vol.

48

No.

5

, pp.

1969

-

1984

, doi:

https://doi.org/10.1111/j.1540-6261.1993.tb05137.x

.

Google Scholar

Crossref

Stein

,

E.M.

and

Stein

,

J.C.

(

1991

), “

Stock price distributions with stochastic volatility: an analytic approach

”,

Review of Financial Studies

, Vol.

4

, pp.

727

-

752

, doi:

https://doi.org/10.1093/rfs/4.4.727

.

Google Scholar

Crossref

Van Moerbeke

,

P.

(

1976

), “

On optimal stopping and free boundary problems

”,

Archive for Rational Mechanics and Analysis

, Vol.

60

No.

2

, pp.

101

-

148

, doi:

https://doi.org/10.1007/bf00250676

.

Google Scholar

Crossref

Wiggins

,

J.B.

(

1987

), “

Option values under stochastic volatility: theory and empirical estimates

”,

Journal of Financial Economics

, Vol.

19

No.

2

, pp.

351

-

372

, doi:

https://doi.org/10.1016/0304-405x(87)90009-2

.

Google Scholar

Crossref

Yi

,

F.

(

2008

), “

American put option with regime-switching volatility (finite time horizon) — variational inequality approach

”,

Mathematical Methods in the Applied Sciences

, Vol.

31

No.

12

, pp.

1461

-

1477

, doi:

https://doi.org/10.1002/mma.984

.

Google Scholar

Crossref

Yousuf

,

M.

,

Khaliq

,

A.Q.M.

and

Liu

,

R.H.

(

2015

), “

Pricing American options under multi-state regime switching with an efficient L − stable method

”,

International Journal of Computer Mathematics

, Vol.

92

No.

12

, pp.

2530

-

2550

, doi:

https://doi.org/10.1080/00207160.2015.1071799

.

Google Scholar

Crossref

Zhang

,

K.

,

Teo

,

K.L.

and

Swartz

,

M.

(

2014

), “

A robust numerical scheme for pricing American options under regime switching based on penalty method

”,

Computational Economics

, Vol.

43

No.

4

, pp.

463

-

483

, doi:

https://doi.org/10.1007/s10614-013-9361-3

.

Google Scholar

Crossref

Appendix

We provide proof of the results in this appendix.

The Snell Envelope. Let us define

Y (t) \equiv e^{- r t} {(K - S (t))}^{+}

By Theorem D.7 in Appendix of Karatzas and Shreve (1998) (Hereafter we will refer to it as “KS”), there is a Q-supermartingale {ξ(t): 0 ≤ t ≤ T} with right continuous and left limit (RCLL) paths such that

ξ (t) \geq Y (t) for all t \in [0, T]

almost surely and for a stopping time v ∈ [0, T]

ξ (v) = ess \sup_{τ \in S_{v, T}} E^{Q} [Y (τ) | F_{v}] a . s .

where $S_{v, T}$ is the set of stopping times between v and T. We know, in particular, $ξ (0) = \sup_{τ \in S_{0, T}} E^{Q} [Y (τ)] = P_{i} (T, S)$ ⁠. By Theorem D.12 in KS, with the stopping time

τ^{*} \equiv {t \in [0, T) : ξ (t) = Y (t)} \land T,

we have ξ(0) = E^QY(τ*).

Proof of Proposition 2.1. The proof is a modification of the proof of Theorem 5.3, Chapter 2 in KS. A self-financing process dominating the American option payoff satisfies the following inequality: there exists some portfolio process (π(t), π_H(t), π_L(t)), t ∈ [0, T] and γ ∈ R such that for every $τ \in S_{0, T}$

Y (τ) \leq γ + \int_{0}^{τ} \frac{1}{S (u)} π (u) σ (u) d \tilde{B} (u) + \int_{0}^{T} \frac{1}{S_{H} (u)} π_{H} (u) d {\tilde{φ}}_{H} (u) + \int_{0}^{T} \frac{1}{S_{L} (u)} π_{L} (u) d {\tilde{φ}}_{L} (u),

(6.1)

where ${\tilde{φ}}_{i} (t) \equiv φ_{i} (t) - λ_{i} t$ and S_i(t) is the price of Θ_i at time i for i ∈ {H, L} (by market completeness assumption Θ_i is tradable in the market).

Let

V_{\min} \equiv \inf {γ \in R : there exists a portfolio process (π (t), π_{H} (t), π_{L} (t)), t \in [0, T], satisfying (6.1)} .

That is, V_min is the smallest initial wealth required such that there exists self-financing strategy whose wealth process satisfies (i) and (ii) in the proposition.

Theorem D.13 in KS asserts that ξ(⋅) = M(⋅) − Λ(⋅), where M(⋅) is a uniformly integrable RCLL martingale under Q and Λ(⋅) is an adapted, continuous, nonincreasing process with Λ(0) = Λ(τ*) = 0 a.s. Because of the assumption of market completeness, the $F_{T}$ -measurable random variable B ≡ S(T)M(T) can be self-financed, i.e. by the martingale representation theorem for a process driven by a Brownian motion and independent Poisson processes (see Theorem 12.33 and 12.35 of Elliot (1980)), there is a triplet $(\hat{π} (\cdot), {\hat{π}}_{H} (\cdot), {\hat{π}}_{L} (\cdot))$ satisfying

M (T) = ξ (0) + \int_{0}^{T} \frac{1}{S (u)} \hat{π} (u) σ (u) d \tilde{B} (u) + \int_{0}^{T} \frac{1}{S_{H} (u)} {\hat{π}}_{H} (u) d {\tilde{φ}}_{H} (u) + \int_{0}^{T} \frac{1}{S_{L} (u)} {\hat{π}}_{L} (u) d {\tilde{φ}}_{L} (u),

where ${\tilde{φ}}_{i} (t) \equiv φ_{i} (t) - λ_{i} t$ and S_i(t) is the price of Θ_i at time i for i ∈ {H, L}.

Taking conditional expectations with respect to $F_{t}$ in the above equation, we have for t ∈ [0, T]

\begin{matrix} Y (t) \leq ξ (t) & = & M (t) - Λ (t) \\ = & ξ (0) - Λ (t) + \int_{0}^{t} \frac{1}{S (u)} \hat{π} (u) σ (u) d \tilde{B} (u) + \int_{0}^{t} \frac{1}{S_{H} (u)} {\hat{π}}_{H} (u) d {\tilde{φ}}_{H} (u) + \int_{0}^{t} \frac{1}{S_{L} (u)} {\hat{π}}_{L} (u) d {\tilde{φ}}_{L} (u) \\ \leq & ξ (0) + \int_{0}^{t} \frac{1}{S (u)} \hat{π} (u) σ (u) d \tilde{B} (u) + \int_{0}^{t} \frac{1}{S_{H} (u)} {\hat{π}}_{H} (u) d {\tilde{φ}}_{H} (u) + \int_{0}^{t} \frac{1}{S_{L} (u)} {\hat{π}}_{L} (u) d {\tilde{φ}}_{L} (u) . \end{matrix}

Therefore,

V_{\min} \leq ξ (0) = P_{i} (0, S),

Suppose that there exists a self-financing strategy whose wealth process X(t) satisfies (i) and (ii) in the statements of the proposition. For any stopping time $τ \in S_{0, T}$

Y (τ) \leq X (0) + \int_{0}^{τ} \frac{1}{S (u)} {\hat{π}}^{'} (u) σ (u) d \tilde{B} (u) + \int_{0}^{τ} \frac{1}{P_{H} (u)} {\hat{π}}_{H}^{'} (u) d {\tilde{φ}}_{H} (u) + \int_{0}^{τ} \frac{1}{P_{L} (u)} {\hat{π}}_{L}^{'} (u) d {\tilde{φ}}_{L} (u),

where $({\hat{π}}^{'} (\cdot), {\hat{π}}_{H}^{'} (\cdot), {\hat{π}}_{L}^{'} (\cdot))$ is a self-financing portfolio for the process X(⋅). Taking expectation in the above equation we know

E^{Q} Y (τ) \leq X (0), \forall τ \in S_{0, T} .

Therefore, P_i(0, S) = ξ(0) ≤ X(0) and this implies that

P_{i} (0, S) \leq V_{\min} .

This completes the proof. ♠

Proof of Lemma 2.2. The proof is essentially the same as that of Proposition 7.1 and Lemma 7.4, Chapter 2 in KS. Therefore, we establish only the proof of the uniform continuity in T here.

Let us define, for i ∈ {H, L},

\begin{matrix} ψ_{i} (t) \equiv E^{Q} [\max_{0 \leq s \leq t} {1 - e^{- r s} H_{i} (s)}^{+}], \\ H_{i} (s) \equiv \exp (r s - \frac{1}{2} \int_{0}^{s} σ^{2} (u) d u + \int_{0}^{s} σ (u) d {\tilde{B}}_{u}) when σ (0) = σ_{i}, \end{matrix}

and

ψ (t) \equiv ψ_{H} (t) \lor ψ_{L} (t) .

Then by the bounded convergence theorem, $\lim_{t ↓ 0} ψ (t) = 0$ ⁠. Suppose that σ(0) = σ_i. For given $S_{0} = S \in [0, \infty)$ ⁠, if we let 0 ≤ T₁ ≤ T₂ and

τ_{2} = \inf {t \in [0, T_{2}) | ξ (t) = e^{- r t} {(K - S H_{i} (t))}^{+}} \land T_{2}, τ_{1} = τ_{2} \land T_{1},

then by the fact S(t) = S H(t), we obtain

\begin{matrix} 0 & \leq P_{i} (T_{2}, S) - P_{i} (T_{1}, S) \leq E^{Q} [e^{- r τ_{2}} {(K - S (τ_{2}))}^{+} - e^{- r τ_{1}} {(K - S (τ_{1}))}^{+}] \\ \leq E^{Q} [{(e^{- r τ_{1}} S (τ_{1}) - e^{- r τ_{2}} S (τ_{2}))}^{+}] \\ \leq E^{Q} [e^{- r τ_{1}} S (τ_{1}) \cdot E^{Q} [{(1 - \min_{T_{1} \leq t \leq T_{2}} \{\exp (- \frac{1}{2} \int_{T_{1}}^{t} σ^{2} (u) d u + \int_{T_{1}}^{t} σ (u) d {\tilde{B}}_{u})\})}^{+}| F_{T_{1}}]] \\ \leq E^{Q} [e^{- r τ_{1}} S (τ_{1})] \cdot ψ (T_{2} - T_{1}) \leq S ψ (T_{2} - T_{1}) . \end{matrix}

Thus, P_i(T, S) is uniformly continuous in T for all $S \in [0, \infty)$ ⁠. ♠

We need the following lemmas in order to prove Theorem 2.3 and Proposition 2.5.

Lemma 6.1.

(P_H(u, S), P_L(u, S)) is a solution to the initial-boundary value problem

\begin{matrix} L f_{i} (u, S) = 0, S > {\underline{S}}_{i} (u), u \in [0, \infty) \\ f_{i} (u, {\underline{S}}_{i} (u)) = K - {\underline{S}}_{i} (u), u \in [0, \infty), \\ f_{i} (0, S) = {(K - S)}^{+}, S \geq {\underline{S}}_{i} (0), \\ \lim_{S \to \infty} \max_{0 \leq u \leq T} | f_{i} (u, S) | = 0 for all T < \infty, \end{matrix}

fori ∈ {H, L}. In particular, the partial derivativesP_i,SS,P_i,S,P_i,texist andarecontinuous in

C_{i}

fori ∈ {H, L} .

Proof. The proof is a modification of the proof of Theorem 7.6 Chapter 2 in KS.

Consider the following initial-boundary value problem: for i ∈ {H, L}

\frac{σ_{i}^{2}}{2} S^{2} f_{i, S S} (u, S) + r S f_{i, S} (u, S) - r f_{i} (u, S) - λ_{i} f_{i} (u, S) - f_{i, u} (u, S) = - λ_{i} P_{j} (u, S), in R

f_{i} = P_{i}, on \partial_{0} R,

where $R = (t_{1}, t_{2}) \times (S_{1}, S_{2}) \subset C_{i}$ and $\partial_{0} R = \partial_{0} R - {t_{2}} \times (S_{1}, S_{2})$ ⁠.

We know by Lemma 2.2 P_i(i ∈ {H, L}) is uniformly Hölder continuous and Theorem 7, Chapter 3 in Friedman (1964) implies that there exists a unique classical solution to the above initial-boundary value problem.

We show that f_i = P_i on $R$ ⁠. Suppose that σ(0) = σ_i and j ∈ {H, L}, j ≠ i,. Let (t₀, S₀) be an element in $R$ and consider a stopping time $τ \in S_{0, t_{0} - t_{1}}$ defined by

τ \equiv \inf {θ \in [0, t_{0} - t_{1}) : (t_{0} - θ, S (θ)) \in \partial_{0} R} \land \inf {θ \in [0, t_{0} - t_{1}) : σ (τ) = σ_{j}} \land (t_{0} - t_{1})

where S(θ) is the price of the stock at time θ with S(0) = S₀ and the process

N (θ) \equiv e^{- r θ} f_{i} (t_{0} - θ, S (θ)) χ_{{σ (θ) = σ_{i}}} + e^{- r θ} P_{j} (t_{0} - θ, S (θ)) χ_{{σ (θ) = σ_{j}}}, θ \in [t_{0} - t_{1}] .

By Ito’s rule N(⋅∧ τ) is a bounded P₀-martingale, and thus

f (t_{0}, S_{0}) = N (0) = E^{Q} N (τ) = \sum_{k \in {H, L}} E^{Q} [e^{- r τ} P_{k} (t_{0} - τ, S (τ)) χ_{{σ (τ) = σ_{k}}}] .

But $(t_{0} - τ, S (t)) \in C_{i}$ implies

τ \leq τ_{S} \equiv \inf {θ \in [0, t_{0}) : ξ (t_{0} - θ) = {(K - S (θ))}^{+}} \land t_{0}

Now the proof of Proposition 2.1 implies that the stopped process

{(e^{- r (t \land τ_{S})} ξ (T - (t \land τ_{S})), F_{t}) : 0 \leq t \leq T}

(6.2)

is a Q-martingale. Now the optional sampling theorem and equation (6.2) yield

\sum_{k \in {H, L}} E^{Q} [e^{- r τ} P_{k} (t_{0} - τ, S (τ)) χ_{{σ (τ) = σ_{k}}}] = P_{i} (t_{0}, S_{0}) .

Thus, f_i and P_i agree on $R$ and hence P_i,SS, P_i,S and P_i,t exist and continuous in $C_{i}$ ⁠.

We now show the last limit in the lemma. Let T ∈ (0, ∞) and for (t, S) ∈ [0, T] × (0, ∞) define

τ_{S} \equiv \inf {θ \in [0, t) : P_{i} (t - θ, S (θ)) = {(K - S (θ))}^{+}} \land t

where S(θ) is the price of the stock at time θ with S(0) = S. Notice that $P_{i} (t, S) = E^{Q} [e^{- r τ_{S}} {(K - S (τ_{S}))}^{+}]$ ⁠. Set

ρ_{S} \equiv {θ \in [0, \infty) : S (θ) \leq K} .

Then,

\begin{matrix} 0 \leq P_{i} (t, S) & \leq & K E^{Q} [χ_{{ρ_{S} \leq t}} e^{- r ρ_{S}}] + E^{Q} [χ_{{ρ_{S} > t}} e^{- r t} {(K - S (t))}^{+}] \\ \leq & K Q [ρ_{S} \leq T] . \end{matrix}

Since $\lim_{S \to \infty} Q [ρ_{S} \leq T] = 0$ ⁠, the last limit in the lemma is valid. ♠

Lemma 6.2.

The Snell envelope ξ(t) has the following decomposition

ξ (t) = M (t) - Λ (t), t \in [0, T],

(6.3)

where

\begin{matrix} M (t) & \equiv & P_{i} (T, S) + \sum_{k \in {H, L}} \int_{0}^{t} e^{- r u} S (u) P_{k, S} (T - u, S (u)) σ_{k} χ_{{σ (u) = σ_{k}}} d \tilde{B} (u) \\ + \sum_{\begin{matrix} k, j \in {H, L}, \\ j \neq k \end{matrix}} \int_{0}^{t} e^{- r u} (P_{j} (T - u, S (u)) - P_{k} (T - u, S (u)) χ_{{σ (u) = σ_{k}}} d {\tilde{φ}}_{k} (u), \end{matrix}

and

Λ (t) \equiv \int_{0}^{t} r K e^{- r u} χ_{{S (u) < {\underline{S}}_{H} & σ (u) = σ_{H}} \cup {S (u) < {\underline{S}}_{L} & σ (u) = σ_{L}}} d u .

Proof. The proof is a modification of the proof of Theorem 7.9, Chapter 2 in KS. Let $ζ : R^{2} \to [0, \infty)$ be a C^∞ function integrating to 1 and having support in [0,1]². For k ∈ {H, L} and ϵ > 0, define

P_{k}^{ϵ} (t, S) \equiv \int_{0}^{\infty} \int_{0}^{\infty} P_{k} (t + ϵ u, S + ϵ v) ζ (u, v) d u d v .

Then, $P_{k}^{ϵ} (t, S)$ is of class C^∞ on (0,∞)². And by the same calculation as in the proof of Theorem 7.9 in KS, we know

\begin{matrix} P_{k, S}^{ϵ} (t, S) & = & \int_{0}^{\infty} \int_{0}^{\infty} P_{k, S} ((t + ϵ u, S + ϵ v) ζ (u, v) d u d v \\ P_{k, S S}^{ϵ} (t, S) & = & \int_{0}^{\infty} \int_{0}^{\infty} P_{k, S S} ((t + ϵ u, S + ϵ v) ζ (u, v) d u d v \\ P_{k, t}^{ϵ} (t, S) & = & \int_{0}^{\infty} \int_{0}^{\infty} P_{k, t} ((t + ϵ u, S + ϵ v) ζ (u, v) d u d v \end{matrix}

These formulas show that $P_{k}^{ϵ}$ and $L P_{k}^{ϵ}$ are bounded on compact subsets of (0,∞)² and

\begin{matrix} P_{k, S} (t, S) & = & \lim_{ϵ ↓ 0} P_{k, S}^{ϵ} (t, S), \\ L P_{k} (t, S) & = & \lim_{ϵ ↓ 0} L P_{k}^{ϵ} (t, S), \forall (t, S) \in {(0, \infty)}^{2}, S \neq {\underline{S}}_{i} (t) . \end{matrix}

According to generalized Ito’s rule,

\begin{matrix} \sum_{k \in {H, L}} e^{- r t} P_{k}^{ϵ} (T - t, S (t)) χ_{{σ (u) = σ_{k}}} \\ = P_{i} (T, S) + \sum_{k \in {H, L}} \int_{0}^{t} e^{- r u} L P_{k}^{ϵ} (T - u, S (u)) χ_{{σ (u) = σ_{k}}} d u \\ + \sum_{k \in {H, L}} \int_{0}^{t} e^{- r u} S (u) P_{k, S}^{ϵ} (T - u, S (u)) σ_{k} χ_{{σ (u) = σ_{k}}} d \tilde{B} (u) \\ + \sum_{\begin{matrix} k, j \in {H, L}, \\ j \neq k \end{matrix}} \int_{0}^{t} e^{- r u} (P_{j}^{ϵ} (T - u, S (u)) - P_{k}^{ϵ} (T - u, S (u)) χ_{{σ (u) = σ_{k}}}) d {\tilde{φ}}_{k} (u), \end{matrix}

(6.4)

for t ∈ [0, T].

For each u ∈ (0, T), we have

Q [{S (u) = {\underline{S}}_{H} & σ (u) = σ_{H}} \cup {S (u) = {\underline{S}}_{L} & σ (u) = σ_{L}}] = 0 .

Therefore, $\sum_{k \in {H, L}} \int_{0}^{t} e^{- r u} L P_{S} (T - u, S (u)) χ_{{σ (u) = σ_{k}}} d u$ is defined and equal to

\int_{0}^{t} e^{- r u} r K χ_{{S (u) < {\underline{S}}_{H} & σ (u) = σ_{H}} \cup {S (u) < {\underline{S}}_{L} & σ (u) = σ_{L}}} d u,

a.s.

Letting ϵ ↓ 0 in equation (6.4), we obtain (6.3), for $t \in [0, T)$ and then for t = T by letting t ↑ T. The process M(⋅) is a martingale because −1 ≤ P_k,S(t, S) ≤ 0 and $E^{Q} [\int_{0}^{T} S {(u)}^{2} d u] < \infty$ ⁠. ♠

Proof of Theorem 2.3. Suppose that σ(0) = σ_i(i ∈ {H, L}). By Lemma 6.2, we have

\begin{matrix} p_{i} (T, S) & = E^{Q} [e^{- r T} {(K - S (T))}^{+}] = E^{Q} [ξ (T)] \\ = E^{Q} [M (T)] - E^{Q} [Λ (T)] = ξ (0) - E^{Q} [Λ (T)] \\ = P_{i} (T, S) - E^{Q} [Λ (T)] . \end{matrix}

Therefore, the decomposition in the theorem follows. ♠

Proof of Proposition 2.5. The proof is a modification of the proof of Theorem 7.12, Chapter 2 in KS. (P_H(u, S), P_L(u, S)) and $({\underline{S}}_{H} (u), {\underline{S}}_{L} (u))$ solve Problem 2.4 by Lemma 6.1.

Suppose (f_H(u, S), f_L(u, S)) is a solution to Problem 2.4 together with $(D_{H} (u), D_{L} (u))$ ⁠.

Fix $(T, S) \in {[0, \infty)}^{2}$ and assume that S(0) = S and σ(0) = σ_i(i ∈ {H, L}). Let us define a process ξ^f(t) by letting

ξ^{f} (t) \equiv \sum_{k \in {H, L}} e^{- r t} f_{k} (T - t, S (t)) χ_{{σ (t) = σ_{k}}} .

By using the same argument as in the proof of Lemma 6.2, we obtain

ξ^{f} (t) = M^{f} (t) - Λ^{f} (t)

where M^f(⋅) is a local martingale defined according as

\begin{matrix} M^{f} (t) & \equiv & f_{i} (T, S) + \sum_{k \in {H, L}} \int_{0}^{t} e^{- r u} S (u) f_{k, S} (T - u, S (u)) σ_{k} χ_{{σ (u) = σ_{k}}} d \tilde{B} (u) \\ + \sum_{\begin{matrix} k, j \in {H, L}, \\ j \neq k \end{matrix}} \int_{0}^{t} e^{- r u} (f_{j} (T - u, S (u)) - f_{k} (T - u, S (u)) χ_{{σ (u) = σ_{k}}} d {\tilde{φ}}_{k} (u), \end{matrix}

and Λ^f(⋅) is nondecreasing and defined as

Λ^{f} (t) \equiv \int_{0}^{t} e^{- r u} r K χ_{{S (u) < D_{H} & σ (u) = σ_{H}} \cup {S (u) < D_{L} & σ (u) = σ_{L}}} d u .

Let ${(τ_{n})}_{n = 1}^{\infty}$ be a sequence of stopping times with τ_n↑ T almost surely and such that {M^f(t ∧ τ_n): 0 ≤ t ≤ T} is a Q-martingale. For any stopping time $τ \in S_{0, T}$ we have

E^{Q} [ξ^{f} (τ \land τ_{n})] = f_{i} (T, S) - E^{Q} Λ^{f} (τ \land τ_{n}) .

By the fact that $\lim_{S ↑ \infty} \max_{0 \leq t \leq T} | f_{k} (u, S) | = 0$ ⁠, k ∈ {H, L}, the function f_k is bounded on $[0, T] \times [0, \infty)$ ⁠, so by taking limit n → ∞ in the above equation the dominated convergence theorem implies

E^{Q} [ξ^{f} (τ)] = f_{i} (T, S) - E^{Q} Λ^{f} (τ), \forall τ \in S_{0, T} .

(6.5)

From the fact f_k(u, S) ≥ (K − S)⁺ for k ∈ {H, L} and $(u, S) \in {[0, \infty)}^{2}$ and the nonnegativity of Λ^f, we have

E^{Q} [e^{- r τ} {(K - S (τ))}^{+}] \leq f_{i} (T, S), \forall τ \in S_{0, T} .

Therefore, f_i(T, S) ≥ P_i(T, S).

Let us now show the reverse inequality. Let

τ_{S}^{'} \equiv T \land \inf {t \geq 0 : (S (t) \leq D_{H} & σ (t) = σ_{H}) or (S (t) \leq D_{L} & σ (t) = σ_{L})}

Then, we have $ξ^{f} (T - τ_{S}^{'}) = {(K - S (τ_{S}^{'}))}^{+}$ and $Λ^{f} (τ_{S}^{'}) = 0$ ⁠, so equation (6.5) implies $E^{Q} [e^{- r τ_{S}^{'}} {(K - S (τ_{S}^{'}))}^{+}] = f_{i} (T, S)$ ⁠. This implies f_i(T, S) ≤ P_i(T, S), and therefore, f_i(T, S) = P_i(T, S).

To show that $D_{i} = {\underline{S}}_{i}$ ⁠, it is enough to show the two open sets $C_{i}^{0}$ and $D_{i}$ are equal, where $C_{i}^{0} \equiv {(t, S) \in {(0, \infty)}^{2} : S > {\underline{S}}_{i} (t)}$ and $D_{i} \equiv {(t, S) \in {(0, \infty)}^{2} : S > D_{i} (t)}$ ⁠. For $(t, S) \in C_{i}^{0}$ ⁠, we have $L P_{i} (t, S) = 0$ which means that $(t, S) \in {\bar{D}}_{i}$ ⁠. Therefore $C_{i}^{0} \subseteq D_{i}$ since $D_{i}$ is open. The other inclusion $D_{i} \subseteq C_{i}^{0}$ can be proved by the same argument. ♠

Proof of Proposition 3.2. We can extend Lemma 6.1 to the case T = ∞ by considering a system of the second-order ordinary differential equations instead of parabolic partial differential equations. Namely, for i ∈ {H, L}

\frac{σ_{i}^{2}}{2} S^{2} f_{i, S S} (S) + r S f_{i, S} (S) - r f_{i} (S) - λ_{i} f_{i} (S) = - λ_{i} P_{j} (S), in R

f_{i} = P_{i}, on \partial R,

where $R = (S_{1}, S_{2}) \subset C_{i}$ and $\partial R = {S_{1}, S_{2}}$ ⁠. We know by the theory of ordinary differential equations that there exists a unique classical solution to the above boundary value problem. Then, the proof of the extended lemma follows in exactly the same way as in the proof of Lemma 6.1.

Now (P_H(S), P_L(S)) and $({\underline{S}}_{H}, {\underline{S}}_{L})$ is a solution to Problem 3.1 by the extended lemma.

Conversely, a solution (f_H(S), f_L(S)) and $(D_{H}, D_{L})$ to Problem 3.1 satisfies f_i(S) = P_i(S) and $D_{i} = {\underline{S}}_{i}$ for i ∈ {H, L} by Theorem 3.1 in Guo and Zhang (2004). ♠

Proof of Lemma 3.3. We first try to find a solution of the form P_H(S) = Sⁿ, P_L(S) = ξSⁿ to (3.1), where n and ξ are some complex numbers. A straightforward calculation shows that if n and ξ satisfy the following system of algebraic equations then P_H(S) = Sⁿ, P_L(S) = ξSⁿ are solutions to (3.1):

\{\begin{cases} \frac{1}{2} σ_{H}^{2} n (n - 1) + r n - r - λ_{H} + λ_{H} ξ = 0 \\ ξ (\frac{1}{2} σ_{L}^{2} n (n - 1) + r n - r - λ_{L}) + λ_{L} = 0 . \end{cases}

(6.6)

Now n satisfies equations (6.6) if and only if it satisfies the fourth-order equation:

(n - 1) [σ_{H}^{2} σ_{L}^{2} n^{3} + (2 r σ_{H}^{2} + 2 r σ_{L}^{2} - σ_{H}^{2} σ_{L}^{2}) n^{2} + 2 (2 r^{2} - r σ_{H}^{2} - r σ_{L}^{2} - λ_{H} σ_{L}^{2} - λ_{L} σ_{H}^{2}) n - 4 r (r + λ_{H} + λ_{L})] = 0 .

(6.7)

We will now show that all the roots of equation (6.7) are real and exactly two of them are negative. Obviously, n = 1 is a one positive root of equation (6.7). Let

I (n) \equiv σ_{H}^{2} σ_{L}^{2} n^{3} + (2 r σ_{H}^{2} + 2 r σ_{L}^{2} - σ_{H}^{2} σ_{L}^{2}) n^{2} + 2 (2 r^{2} - r σ_{H}^{2} - r σ_{L}^{2} - λ_{H} σ_{L}^{2} - λ_{L} σ_{H}^{2}) n - 4 r (r + λ_{H} + λ_{L}) .

Then,

I (0) = - 4 r (r + λ_{1} + λ_{2}) < 0, I (- \frac{2 r}{σ_{L}^{2}}) = \frac{4 r λ_{2} (σ_{H}^{2} - σ_{L}^{2})}{σ_{L}^{2}} > 0 .

Therefore, there exist two distinct real negative roots and two (possibly equal) real positive roots of equation (6.7).

Note that among the solutions of the form P_H(S) = Sⁿ, P_L(S) = ξSⁿ, only those with negative n satisfy boundary conditions $\lim_{S ↑ \infty} P_{i} (S) = 0$ for i ∈ {H, L}.

We know that the homogeneous system of equations (3.1) has a solution space of dimension 4. Suppose that the positive real roots to (6.7) are distinct. Then the basis to the solution space consists of pairs of functions P_H(S) = Sⁿ, P_L(S) = ξSⁿ with n, ξ satisfying (6.6). Solutions to (3.1) satisfying the boundary conditions $\lim_{S ↑ \infty} P_{i} (S) = 0$ for i = H, L form a two-dimensional subspace whose basis consists of pairs of functions P_H(S) = Sⁿ, P_L(S) = ξSⁿ with n < 0 and ξ satisfying (6.6). Therefore, a general solution takes the form in equation (3.3). A straightforward extension of the proof shows that the result is also valid for the case where the positive real roots to (6.7) are both equal to 1.

It is easy to show that $I (- \frac{2 r}{σ_{H}^{2}}) < 0$ ⁠. This proves inequalities for n₁ and n₂ in (3.4).

Let $J (n) \equiv \frac{1}{2} σ_{H}^{2} n (n - 1) + r n - r$ ⁠. Then, it is easy to show that J′(n) < 0 for $n \leq - \frac{2 r}{σ_{H}^{2}}$ and $J (- \frac{2 r}{σ_{H}^{2}}) = 0$ ⁠. This proves inequalities 1 > ξ₁ > ξ₂. Also, the first and second equations in (3.1) can be converted into

1 - ξ = (n - 1) (\frac{1}{2} σ_{H}^{2} n + r) / λ_{H}, 1 - \frac{1}{ξ} = (n - 1) (\frac{1}{2} σ_{L}^{2} n + r) / λ_{L},

respectively. Dividing the first equation by the second, we obtain

ξ = - \frac{λ_{H} (\frac{1}{2} σ_{L}^{2} n + r)}{λ_{L} (\frac{1}{2} σ_{H}^{2} n + r)},

and this fact together with $n_{2} < - 2 r / σ_{L}^{2} < n_{1} < - 2 r / σ_{H}^{2}$ implies that ξ₁ > 0 > ξ₂.

The proof is now complete. ♠

Proof of Lemma 3.4.Let us assume ${\underline{S}}_{m} < S < {\underline{S}}_{M}$ ⁠. Then, by (3.2) we know that P_M(S) = K − S. Therefore, P_m(S) satisfies

\frac{1}{2} σ_{m}^{2} S^{2} P_{m, S S} + r S P_{m, S} - (r + λ_{m}) P_{m} + λ_{m} (K - S) = 0 .

Now a solution to the above equation takes the form in equation (3.5), where m₁, m₂ are solutions to the characteristic equation

\frac{1}{2} σ_{m}^{2} z^{2} + (r - \frac{1}{2} σ_{m}^{2}) z - (r + λ_{m}) = 0 .

Therefore, m₁, m₂ satisfy inequalities in (3.6). This completes the proof. ♠

Proof of Proposition 3.5. From the boundary conditions for P_m(S) at ${\underline{S}}_{m}$ ⁠, we obtain equation (3.9), and from the boundary conditions for P_M(S) at ${\underline{S}}_{M}$ ⁠, we get equation (3.8). Finally, by (3.9) and (3.8) the C²-property of P_m(S) at ${\underline{S}}_{M}$ is equivalent to equation (3.7). ♠

Proof of Corollary 3.6. (i) (3.6) and (3.9) yield the property for D₁ and D₂ immediately.

(ii) Equation (3.8) implies that both C₁ and C₂ are non-zero. ♠

Proof of Theorem 3.7. Assume that the converse is true, i.e. ${\underline{S}}_{H} \geq {\underline{S}}_{L}$ ⁠.

By Lemma 2.2, the value of the American put is a convex function of the underlying stock price S (The lemma is stated for a put option with finite expiry, but the proof is valid also for a perpetual put).

Therefore, P_i,SS(S) ≥ 0 for all $S > {\underline{S}}_{i}$ for i = H, L. Convexity of P_L(S) and the boundary condition $P_{L, S} ({\underline{S}}_{L}) = - 1$ implies that

P_{H} (S) = K - S \leq P_{L} (S) for {\underline{S}}_{L} \leq S \leq {\underline{S}}_{H},

and

P_{L, S} ({\underline{S}}_{H}) - P_{H, S} ({\underline{S}}_{H}) \geq 0 .

(6.8)

Therefore,

P_{L} ({\underline{S}}_{H}) - P_{H} ({\underline{S}}_{H}) \equiv l_{1} \geq 0 .

(6.9)

From (3.3), we know $\lim_{S ↑ \infty} P_{H, S S} (S) = 0$ ⁠, and hence there exists l₂ > 0 such that 0 ≤ P_H,SS(S) ≤ l₂.

Let g(S) ≡ P_L(S) − P_H(S). Then, for $S \geq {\underline{S}}_{H}$ ⁠, we derive the following equation from (3.1), (3.2), (6.8) and (6.9)

\frac{1}{2} σ_{L}^{2} S^{2} g^{″} (S) + r S g^{'} (S) - (r + λ_{H} + λ_{L}) g (S) = \frac{1}{2} (σ_{H}^{2} - σ_{L}^{2}) S^{2} P_{H, S S} (S),

with boundary conditions $g ({\underline{S}}_{H}) = l_{1}, g^{'} ({\underline{S}}_{H}) \geq 0, \lim_{S ↑ \infty} g (S) = 0$ ⁠. Let $L (S) \equiv \frac{σ_{H}^{2} - σ_{L}^{2}}{σ_{L}^{2}} P_{H, S S} (S)$ ⁠. Then 0 ≤ L(S) ≤ l₃, for some constant l₃, since 0 ≤ P_H,SS(S) ≤ l₂. Applying a standard argument for an Euler equation, a general solution takes the following form

g (S) = c_{1} S^{k_{1}} + c_{2} S^{k_{2}} - S^{k_{1}} \int_{{\underline{S}}_{H}}^{S} \frac{L (t)}{(k_{2} - k_{1}) σ_{L}^{2} t^{k_{1} - 1}} d t + S^{k_{2}} \int_{{\underline{S}}_{H}}^{S} \frac{L (t)}{(k_{2} - k_{1}) σ_{L}^{2} t^{k_{2} - 1}} d t

for some constants c₁ and c₂, where.

$k_{1} \equiv \frac{r - σ_{L}^{2} / 2 + \sqrt{{(r - σ_{L}^{2} / 2)}^{2} + 2 σ_{L}^{2} (r + λ_{H} + λ_{L})}}{σ_{L}^{2}}, k_{2} \equiv \frac{r - σ_{L}^{2} / 2 - \sqrt{{(r - σ_{L}^{2} / 2)}^{2} + 2 σ_{L}^{2} (r + λ_{H} + λ_{L})}}{σ_{L}^{2}} .$

Note that k₁ > 1 and k₂ < 0.

By (3.3), we know that

|\int_{{\underline{S}}_{H}}^{\infty} \frac{L (t)}{(k_{2} - k_{1}) σ_{L}^{2} t^{k_{1} - 1}} d t| < \infty .

From (3.3) and the fact that n₁, n₂, k₂ < 0, we know

\lim_{S ↑ \infty} S^{k_{2}} \int_{{\underline{S}}_{H}}^{S} \frac{L (t)}{(k_{2} - k_{1}) σ_{L}^{2} t^{k_{2} - 1}} d t = 0 .

(6.10)

(6.10), the fact that k₁ > 1 and the boundary condition $\lim_{S ↑ \infty} g (S) = 0$ force c₁ to satisfy the following

c_{1} = \int_{{\underline{S}}_{H}}^{\infty} \frac{L (t)}{(k_{2} - k_{1}) σ_{L}^{2} t^{k_{1} - 1}} d t .

Here, the improper integral in the right-hand side is finite because of equation (3.3) and the fact that n₁, n₂ < 0, k₁ > 0.

Therefore,

l_{1} = g ({\underline{S}}_{H}) = {\underline{S}}_{H}^{k_{1}} \int_{{\underline{S}}_{H}}^{\infty} \frac{L (t)}{(k_{2} - k_{1}) σ_{L}^{2} t^{k_{1} - 1}} d t + c_{2} {\underline{S}}_{H}^{k_{2}},

or equivalently,

c_{2} = l_{1} {\underline{S}}_{H}^{- k_{2}} - {\underline{S}}_{H}^{k_{1} - k_{2}} \int_{{\underline{S}}_{H}}^{\infty} \frac{L (t)}{(k_{2} - k_{1}) σ_{L}^{2} t^{k_{1} - 1}} .

But,

g^{'} ({\underline{S}}_{H}) = k_{1} c_{1} {\underline{S}}_{H}^{k_{1} - 1} + k_{2} c_{2} {\underline{S}}_{H}^{k_{2} - 1} = - {\underline{S}}_{H}^{k_{1} - 1} \int_{{\underline{S}}_{H}}^{\infty} \frac{L (t)}{σ_{L}^{2} t^{k_{1} - 1}} d t + l_{1} k_{2} {\underline{S}}_{H}^{- 1} < 0

where the last inequality follows from the fact that L(t) is not identically zero for $S \geq {\underline{S}}_{H}$ because of Corollary 3.6 (ii). This is contradictory to the boundary condition $g^{'} ({\underline{S}}_{H}) \geq 0$ ⁠.

The proof is complete. ♠

Proof of Theorem 3.8. Let h(S) ≡ P_H(S) − P_L(S). Then, by an argument similar to the proof of Theorem 3.7, we know that h(S) satisfies the following equation for $S > {\underline{S}}_{L}$

\frac{1}{2} σ_{H}^{2} S^{2} h^{″} (S) + r S h^{'} (S) - (r + λ_{H} + λ_{L}) h (S) = \frac{1}{2} (σ_{L}^{2} - σ_{H}^{2}) S^{2} P_{L, S S} (S),

with boundary conditions $h ({\underline{S}}_{L}) \equiv l \geq 0, h^{'} ({\underline{S}}_{L}) \geq 0, \lim_{S ↑ \infty} h (S) = 0$ ⁠. Let $L (S) \equiv \frac{σ_{L}^{2} - σ_{H}^{2}}{σ_{H}^{2}} P_{L, S S} (S)$ ⁠. Then, $- \tilde{l} \leq L (S) \leq 0$ for some positive constant $\tilde{l}$ ⁠.

Applying a standard argument, we know

h (S) = c_{1} S^{k_{1}} + c_{2} S^{k_{2}} - S^{k_{1}} \int_{{\underline{S}}_{L}}^{S} \frac{L (t)}{(k_{2} - k_{1}) σ_{H}^{2} t^{k_{1} - 1}} d t + S^{k_{2}} \int_{{\underline{S}}_{L}}^{S} \frac{L (t)}{(k_{2} - k_{1}) σ_{H}^{2} t^{k_{2} - 1}} d t

(6.11)

for some constants c₁ and c₂, where

\begin{matrix} k_{1} = \frac{r - σ_{H}^{2} / 2 + \sqrt{{(r - σ_{H}^{2} / 2)}^{2} + 2 σ_{H}^{2} (r + λ_{H} + λ_{L})}}{σ_{H}^{2}}, \\ k_{2} = \frac{r - σ_{H}^{2} / 2 - \sqrt{{(r - σ_{H}^{2} / 2)}^{2} + 2 σ_{H}^{2} (r + λ_{H} + λ_{L})}}{σ_{H}^{2}} . \end{matrix}

Here, k₁ > 1 and k₂ < 0.

By an argument similar to the proof of Theorem 3.7, we can show that

c_{1} = \int_{{\underline{S}}_{L}}^{\infty} \frac{L (t)}{(k_{2} - k_{1}) σ_{H}^{2} t^{k_{1} - 1}} d t \geq 0,

and

c_{2} = l {\underline{S}}_{L}^{- k_{2}} - {\underline{S}}_{L}^{k_{1} - k_{2}} \int_{{\underline{S}}_{L}}^{\infty} \frac{L (t)}{(k_{2} - k_{1}) σ_{H}^{2} t^{k_{1} - 1}} .

Plugging these into (6.11), we know that h must be of the form

\begin{matrix} h (S) = & [S^{k_{1}} - {(\frac{S}{{\underline{S}}_{L}})}^{k_{2}} {\underline{S}}_{L}^{k_{1}}] \int_{S}^{\infty} \frac{L (t)}{(k_{2} - k_{1}) σ_{H}^{2} t^{k_{1} - 1}} d t \\ + l {(\frac{S}{{\underline{S}}_{L}})}^{k_{2}} + {(\frac{S}{{\underline{S}}_{L}})}^{k_{2}} {\underline{S}}_{L} \int_{{\underline{S}}_{L}}^{S} \frac{L (t)}{(k_{2} - k_{1}) σ_{H}^{2}} [{(\frac{{\underline{S}}_{L}}{t})}^{k_{2} - 1} - {(\frac{{\underline{S}}_{L}}{t})}^{k_{1} - 1}] d t . \end{matrix}

Since $S^{k_{1}} - {(\frac{S}{{\underline{S}}_{L}})}^{k_{2}} {\underline{S}}_{L}^{k_{1}} \geq 0$ for $S \geq {\underline{S}}_{L}$ and ${(\frac{{\underline{S}}_{L}}{t})}^{k_{2} - 1} \geq {(\frac{{\underline{S}}_{L}}{t})}^{k_{1} - 1}$ for all $t \geq {\underline{S}}_{L}$ ⁠, h(S) ≥ 0 for all $S > {\underline{S}}_{L}$ ⁠. Strict inequality for $S > {\underline{S}}_{L}$ also follows from the above equation.

For $S \leq {\underline{S}}_{L}$ ⁠, P_L(S) = K − S ≤ P_H(S) by convexity of P_H(S) and the boundary condition $P_{H, S} ({\underline{S}}_{H}) = - 1$ ⁠. The strict inequality for ${\underline{S}}_{H} < S \leq {\underline{S}}_{L}$ follows from equation (3.5) and Corollary 3.6.

The proof is complete. ♠

Proof of Proposition 4.1. The proof is a slight modification of the proofs of Lemmas 3.3 and 3.4. ♠

Proof of Theorem 4.2. The proof proceeds in the same way as that of 3.7 and 3.8. ♠

Proof of Theorem 4.3. (i) Equations in (4.1) are derived in the next in this appendix. Since $P_{m}^{k - 1} (t)$ is bounded and n₂ < n₁ < 0 < 1 < n₄ < n₃, from (6.16) we can show that

\lim_{S \to \infty} | S^{n_{h}} u_{h}^{k} (S) | < \infty for h = 1,2,

and

\lim_{S \to \infty} | u_{h}^{k} (S) | < \infty for h = 3,4 .

Therefore, (4.2) follows immediately.

(ii) This result follows from a standard argument for a second-order ordinary differential equations. ♠

Proof of Theorem 4.4.

Use mathematical induction on k. Theorem 4.2 establishes the proof for k = 1. The rest of the proof is essentially the same as the proofs of Theorem 3.7 and Theorem 3.8. ♠

Derivation of Equation (4.1). Note that the homogeneous solution of the system of ordinary differential equations is of the following form with four undetermined constants $C_{h}^{k}$ ’s:

φ_{i}^{k} (S) = C_{1}^{k} S^{n_{1}} + C_{2}^{k} S^{n_{2}} + C_{3}^{k} S^{n_{3}} + C_{4}^{k} S^{n_{4}}, φ_{j}^{k} (S) = C_{1}^{k} ξ_{1} S^{n_{1}} + C_{2}^{k} ξ_{2} S^{n_{2}} + C_{3}^{k} ξ_{3} S^{n_{3}} + C_{4}^{k} ξ_{4} S^{n_{4}} .

In order to find the particular solution, we use the method of variation of parameters. We derive the matrix equation

X^{'} (S) = P (S) X (S) + G (S),

(6.12)

where, for particular solutions $ϕ_{i}^{k}$ and $ϕ_{j}^{k}$ ⁠,

X (z) = {(ϕ_{i}^{k} (z), ϕ_{j}^{k} (z), {ϕ_{i}^{k}}^{'} (z), {ϕ_{j}^{k}}^{'} (z))}^{⊤}, G (z) = {(0, 0, - \frac{2 β P_{i}^{k - 1} (z)}{σ_{i}^{2} z^{2}}, - \frac{2 β P_{j}^{k - 1} (z)}{σ_{j}^{2} z^{2}})}^{⊤}

and

P (z) = (\begin{matrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ \frac{2 (r + λ_{i} + β)}{σ_{i}^{2} z^{2}} & \frac{- 2 λ_{i}}{σ_{i}^{2} z^{2}} & \frac{- 2 r}{σ_{i}^{2} z} & 0 \\ \frac{2 (r + λ_{j} + β)}{σ_{j}^{2} z^{2}} & \frac{- 2 λ_{j}}{σ_{j}^{2} z^{2}} & 0 & \frac{- 2 r}{σ_{j}^{2} z} \end{matrix}) .

Also set the fundamental solution Φ to be

Φ (z) = (\begin{matrix} z^{n_{1}} & z^{n_{2}} & z^{n_{3}} & z^{n_{4}} \\ ξ_{1} z^{n_{1}} & ξ_{2} z^{n_{2}} & ξ_{3} z^{n_{3}} & ξ_{4} z^{n_{4}} \\ n_{1} z^{n_{1} - 1} & n_{2} z^{n_{2} - 1} & n_{3} z^{n_{3} - 1} & n_{4} z^{n_{4} - 1} \\ ξ_{1} n_{1} z^{n_{1} - 1} & ξ_{2} n_{2} z^{n_{2} - 1} & ξ_{3} n_{3} z^{n_{3} - 1} & ξ_{4} n_{4} z^{n_{4} - 1} \end{matrix}) .

(6.13)

Note that the fundamental solution satisfies the following equation

Φ^{'} (S) = P (S) Φ (S) .

(6.14)

Also note that Φ(z) is invertible for all z > 0 if and only if the determinant of Φ(1) is not equal to zero.

The method of variation of parameters is started by letting

X (S) \equiv Φ (S) U (S) \equiv Φ (S) {(u_{1}^{k} (S), u_{2}^{k} (S), u_{3}^{k} (S), u_{4}^{k} (S))}^{⊤}

(6.15)

From (6.12) and (6.14), we obtain

Φ (S) U^{'} (S) = G (S),

or equivalently,

U (S) = \int_{{\underline{S}}_{j}^{k}}^{S} Φ^{- 1} (t) G (t) d t .

Thus the particular solution can be represented as

(\begin{matrix} ϕ_{i}^{k} (S) \\ ϕ_{j}^{k} (S) \end{matrix}) = (\begin{matrix} \sum_{h = 1}^{4} S^{n_{h}} \cdot u_{h}^{k} (S) \\ \sum_{h = 1}^{4} ξ_{h} S^{n_{h}} \cdot u_{h}^{k} (S) \end{matrix}) .

Consequently, the solution is of the form

(\begin{matrix} P_{i}^{k} (S) \\ P_{j}^{k} (S) \end{matrix}) = (\begin{matrix} φ_{i}^{k} (S) \\ φ_{j}^{k} (S) \end{matrix}) + (\begin{matrix} ϕ_{i}^{k} (S) \\ ϕ_{j}^{k} (S) \end{matrix}) = (\begin{matrix} \sum_{h = 1}^{4} (C_{h}^{k} + u_{h}^{k} (S)) S^{n_{h}} \\ \sum_{h = 1}^{4} ξ_{h} (C_{h}^{k} + u_{h}^{k} (S)) S^{n_{h}} \end{matrix}) .

After immediate calculation, we obtain

u_{h}^{k} (S) = \int_{{\underline{S}}_{j}^{k}}^{S} t^{- n_{h} - 1} (A_{h} P_{i}^{k - 1} (t) + B_{h} P_{j}^{k - 1} (t)) d t,

(6.16)

for some constants A_h and B_h. ♠

2024

Bong-Gyu Jang and Hyeng Keun Koo

Published in Journal of Derivatives and Quantitative Studies: 선물연구. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence maybe seen at http://creativecommons.org/licences/by/4.0/legalcode

American put options with regime-switching volatility

1. Introduction

2. Model

3. Perpetual American puts

4. American puts with finite expiry

4.1 Exponentially distributed time

4.2 Erlang-distributed time

5. Implementation and numerical results

5.1 Richardson extrapolation with randomized put values

5.2 Characteristics of American put values

6. Conclusion

Notes

References

Appendix

Email Alerts

Cited By

American put options with regime-switching volatility Open Access

1. Introduction

2. Model

3. Perpetual American puts

4. American puts with finite expiry

4.1 Exponentially distributed time

4.2 Erlang-distributed time

5. Implementation and numerical results

5.1 Richardson extrapolation with randomized put values

5.2 Characteristics of American put values

6. Conclusion

Notes

References

Appendix

Email Alerts

Suggested Reading

Related Chapters

Recommended for you

Cited By

American put options with regime-switching volatility