Open AccessCommunication

Pricing American Options with a Non-Constant Penalty Parameter

Applied Mathematics & Numerical Analysis Group, University of Wuppertal, 42119 Wuppertal , Germany

Department of Applied Mathematics and Statistics, Division of Applied Mathematics, Comenius University, 842 48 Bratislava, Slovakia

Author to whom correspondence should be addressed.

J. Risk Financial Manag. 2020, 13(6), 124; https://doi.org/10.3390/jrfm13060124

Submission received: 11 May 2020 / Revised: 3 June 2020 / Accepted: 10 June 2020 / Published: 13 June 2020

(This article belongs to the Special Issue Option Pricing)

Download Versions Notes

Abstract

As the American early exercise results in a free boundary problem, in this article we add a penalty term to obtain a partial differential equation, and we also focus on an improved definition of the penalty term for American options. We replace the constant penalty parameter with a time-dependent function. The novelty and advantage of our approach consists in introducing a bounded, time-dependent penalty function, enabling us to construct an efficient, stable, and adaptive numerical approximation scheme, while in contrast, the existing standard approach to the penalisation of the American put option-free boundary problem involves a constant penalty parameter. To gain insight into the accuracy of our proposed extension, we compare the solution of the extension to standard reference solutions from the literature. This illustrates the improvement of using a penalty function instead of a penalising constant.

Keywords:

American Options; PDE option pricing; Penalty term; projected SOR; penalization strategy

1. Introduction

As American options give the holder the right to exercise the option before and at maturity, this leads to free boundary value problems which have to be solved numerically. Several schemes for solving American option problems have been proposed, such as the projected SOR scheme (Cryer 1971; Wilmott et al. 1993), the binomial method, front-fixing schemes (Nielsen et al. 2002), the power penalty method (Wang et al. 2006), and Monte Carlo simulation techniques. These schemes compute the free boundary value implicitly. Other researchers focused on an explicit representation of the free boundary value (Lauko and Ševčovič 2010; Stamicar et al. 1999; Zhu 2006).

Another approach to solving American option problems is to add a penalty term to the problem (Günther and Jüngel 2010; Nielsen et al. 2002; Wang et al. 2006). A penalty term forces the problem to fulfill the free boundary constraint asymptotically. If the free boundary constraint is fulfilled, the penalty term is zero, otherwise it penalizes the problem with a factor. Until now, the penalty term included a penalisation constant being roughly estimated by an educated guess or an optimization. The recent approach of the power penalty function (Wang et al. 2006) improves the convergence and accuracy by introducing a non-constant penalization function of a given power form. Here, we further improve the performance of the penalization by replacing the penalty constant with a general penalty function, allowing for an adaptive penalization to be determined in a consistent way while solving the partial differential equation (PDE). An explicit formulation for the derivation of the parameters of the penalty function is part of our further research.

Once we add a penalty term to the American option problem, the problem reduces to a PDE on a fixed domain and we can apply standard numerical methods, such as finite difference methods or finite element methods. In front-fixing schemes, a change of variables is used to obtain a PDE on a fixed domain. Through the change of variables, the free boundary value is tracked implicitly (Han and Wu 2003; Schwartz 1977). In our approach, the domain is fixed by the initial guess of the free boundary value obtained by the intersection point of the payoff and the exact solution of the Black–Scholes equation. The penalisation of a time-step depends on the solution itself in comparison with the payoff (Günther and Jüngel 2010; Nielsen et al. 2002). Our presented penalty term is bounded by an initial guess of the free boundary value, and is thus independent from the solution of the last time-step. The usual approach to penalisation of the American put option-free boundary problem involves a small parameter, making the numerical analysis harder. The novelty and advantage of our approach consists in introducing a bounded penalty function, enabling us to construct an efficient and stable numerical approximation scheme.

The outline of this article is as follows. Section 2 reviews the mathematical modeling for an American option with and without a penalty term. As a benchmark, we chose the classical Black–Scholes equation and use a variable transformation to simplify our computations and add the penalty function. In Section 3 the model is discretized, and the numerical results of the different test cases are presented in Section 4. In Section 5 we conclude this work with a brief outlook.

2. Mathematical Modeling

American options are more expensive than European options, as American options give the holder the right to exercise the option before the maturity T. In the following, we will focus on American put options for clearness of the idea, but as it can be seen in the numerical results, all assumptions also hold analogously for the American call options. For pricing an American put option P with the Black–Scholes model, we seek a pair of functions

(P (S, t), S_{f} (t))

such that

\begin{matrix} L_{B S M} [P] \equiv \frac{\partial P}{\partial t} & + \frac{σ^{2}}{2} S^{2} \frac{\partial^{2} P}{\partial S^{2}} + (r - q) S \frac{\partial P}{\partial S} - r P \leq 0, 0 \leq t \leq T, \\ {(K - S)}^{+} & = P (S, t) for S \leq S_{f} (t), \\ {(K - S)}^{+} & < P (S, t) for S > S_{f} (t), \end{matrix}

(1)

where K denotes the predefined strike price, r is the risk-free interest rate, q is the dividend rate,

σ

is the volatility, S is the price of an asset, and

S_{f} (t)

is the free boundary value at time t with

0 \leq t \leq T

. In (1), we used the standard notation

{(f)}^{+} : = max (f, 0)

. The differential operator appearing in the Black–Scholes PDE (1) is abbreviated by

L_{B S M}

The terminal condition at maturity

t = T

reads

P (S, T) = {(K - S)}^{+},

(2)

and the "spatial" boundary conditions at

S = S_{f} (t)

S \to \infty

, are given by

P (S_{f} (t), t) = {(K - S_{f} (t))}^{+}, \frac{\partial P}{\partial S} (S_{f} (t), t) = - 1, lim_{S \to \infty} P (S, t) = 0, 0 \leq t \leq T .

The American Put option problem can be split into two regions: the exercise region

0 \leq S \leq S_{f} (t)

, and the holding region

S > S_{f} (t)

. In the holding region, the American Put option problem follows the Black–Scholes equation for the European Put option

L_{B S M} [P] (S, t) = 0, S > S_{f} (t) .

For the exercise region,

P (S, t) = {(K - S)}^{+}

has to hold. By inserting this condition into

L_{B S M} [P] (S, t)

, we obtain

L_{B S M} [{(K - S)}^{+}] = - r S - r (K - S) = - r K < 0,

(3)

fulfilling the condition (1). We obtain a case system for the American put option problem

L_{B S M} [P] (S, t) = g (S, t) \{\begin{matrix} - r K, & 0 < S \leq S_{f} (t), \\ 0, & S > S_{f} (t) . \end{matrix}

(4)

To simplify the computational effort, we transform the problem (4) into the result that we obtain on the left-hand side of the heat equation. This transformation is well-known for the European options, since the heat equation has an analytical solution and is analysed very often. Since the right-hand side is in comparison to the European Put options unequal to zero, we sketch the transformation in more detail.

2.1. Transformation

For the rewriting of the problem (4), we use the standard transformations

x = ln (\frac{S}{K}), τ = \frac{σ^{2}}{2} (T - t), v (x, τ) = \frac{P (S, t)}{K}

(5)

and use the chain rule

- \frac{σ^{2}}{2} K v_{τ} + \frac{σ^{2}}{2} S^{2} \frac{K}{S} (v_{x x} - v_{x}) + r S \frac{K}{S} v_{x} - r K v = \{\begin{matrix} - r K, & x \leq x_{f} (τ), \\ 0, & x > x_{f} (τ), \end{matrix}

where

x_{f} (τ) = ln (S_{f} (τ) / K)

denotes the free boundary value for x. We simplify this equation and divide by

- σ / 2

v_{τ} - v_{x x} + (1 - k) v_{x} + k v = \{\begin{matrix} k, & x \leq x_{f} (τ), \\ 0, & x > x_{f} (τ), \end{matrix}

(6)

with

k = \frac{2 r}{σ^{2}}

. In the last step, we use the transformation

u (x, τ) = exp (α x + β τ) v (x, τ) = e (x, τ) v (x, τ)

, where

α = - \frac{1}{2} (k - 1) and β = - \frac{1}{4} {(k + 1)}^{2} .

(7)

We obtain the transformed problem

u_{τ} - u_{x x} = \tilde{g} (x, τ) = exp (- α x - β τ) \{\begin{matrix} k & x \leq x_{f} (τ), \\ 0, & x > x_{f} (τ) . \end{matrix}

(8)

with the transformed terminal condition (2)

f (x, τ) = exp (- α x - β τ) {(1 - exp (x))}^{+},

(9)

supplied with the initial and boundary conditions

u (x, 0) = f (x, 0), x \in R, lim_{x \to \pm \infty} (u (x, τ) - f (x, τ)) = 0, 0 \leq τ \leq T .

(10)

2.2. The Penalty Term

For the penalty term, we chose an affine function,

δ (t) = a t + b,

(11)

in a way that the penalisation function is given by

p (S, t) = δ (t) g (S, t)

. This choice preserves the novelty of this approach, as the known penalty terms are neither bounded nor independent from the solution itself (Günther and Jüngel 2010; Nielsen et al. 2002). The penalised case system for (4) reads

L_{B S M} [P] (S, t) = δ (t) \cdot g (S, t) .

(12)

t = T

, the terminal condition (2) and (3) has to hold, and we set

p (t) = a (T - t) + 1

. Similarly, we chose an affine function

\tilde{δ} (τ) = \tilde{a} τ + b

as a penalisation term for the penalisation function

\tilde{p} = \tilde{δ} (τ) \tilde{g} (x, τ)

for the transformed case system (8). Since the relation between

u (x, τ)

and

P (S, t)

is given by

u (x, τ) = \frac{1}{K} exp (- α x - β τ) P (S, t),

the same relation has to hold for the penalised right-hand side. We focus only on the exercise region since the right-hand side of the holding region is always zero. As with t, we step forward in time, and with

τ

backward,

\frac{1}{K} exp (- α x - β τ) p (S, t) = - \tilde{p} (x, τ)

can be assumed. We thus obtain

\begin{matrix} \frac{- r K}{K} exp (- α x - β τ) (a (T - t) + 1) & = & - (\tilde{a} τ + \tilde{b}) exp (- α x - β τ) \frac{2 r}{σ^{2}} \\ a (T - t) + 1 & = & (\tilde{a} \frac{σ^{2}}{2} (T - t) + b) \frac{2}{σ^{2}} \\ \frac{σ^{2}}{2} & = & b . \end{matrix}

The transformed case system with the bounded penalisation function is given by

u_{τ} - u_{x x} = \tilde{p} (x, τ) = (\tilde{a} + \frac{σ^{2}}{2}) \{\begin{matrix} k exp (- α x - β τ), & x \leq x_{f} (τ) \\ 0, & x > x_{f} (τ) \end{matrix},

(13)

supplied with the initial and boundary conditions

u (x, 0) = f (x, 0), x \in R, lim_{x \to \pm \infty} (u (x, τ) - f (x, τ)) = 0, 0 \leq τ \leq T .

(14)

Our approach requires an initial guess for the free boundary value. Since the price of an American option is larger than or equal to the price of the corresponding European option, we can use the intersection point

{\bar{x}}_{f} (τ)

between the payoff and the solution of the European Put option as the initial guess for the free boundary value

x_{f} (τ)

, with the result that the bounded penalty term forces the PDE (12) to fulfill the conditions (1) asymptotically. As the approximation of the free boundary and the analysis for the approximation is one of the main research areas in the field of American options, well-known approximation formulas can be found in the literature, such as the approximations from Zhu (2006), Stamicar et al. (1999), and Evans et al. (2002). Each of those formulations can be considered to gain the initial values for the free boundary. Since the penalisation is based on the initial choice of the free boundary, the choice of the formula for the computation of the initial free boundary has a large effect on the accuracy of the method. A detailed analysis of the choice of the initial free boundary computed by different approximation formulas from the literature is part of our future research. After computing the penalty term with the initial free boundary value, we solve the penalised heat Equation (13). The obtained solution is the solution of the American option problem.

3. Discretization

Let us introduce a temporal discretization

τ_{j} = T - j Δ τ

Δ τ = T / M

j = 0, \dots, M

, and a spatial grid between the points

x_{m i n}

and

x_{max}

x_{i} = x_{min} + i Δ x, Δ x = \frac{x_{max} - x_{min}}{N}, i = 0, \dots, N .

We use the finite difference

θ

-scheme for discretization and simplify the notation by

α_{1} = θ \frac{Δ τ}{{(Δ x)}^{2}} and α_{2} = (1 - θ) \frac{Δ τ}{{(Δ x)}^{2}} with 0 \leq θ \leq 1 .

We obtain

w^{j} = {(w_{1}^{j}, \dots, w_{N - 1}^{j})}^{⊤}

with

w_{i}^{j}

as the approximation for

u (x_{i}, τ_{j})

f^{j} = {(f_{1}^{j}, \dots, f_{N - 1}^{j})}^{⊤}

with

f_{i}^{j} \sim f (x_{i}, τ_{j})

and the diagonal matrices A and B

A = diag (- α_{1}, 2 α_{1} + 1, - α_{1}), B = diag (α_{2}, - 2 α_{2} + 1, α_{2}),

as well as the vector

d^{j}

containing the boundary values

{(d^{j})}^{⊤} = (α_{1} w_{0}^{j + 1} + α_{2} w_{0}^{j}, 0, \dots, 0, α_{1} w_{N + 1}^{j + 1} + α_{2} w_{N + 1}^{j}) .

The discretized penalty term (13) is given by

\begin{matrix} p^{j} = δ^{j} \cdot g^{j} with δ^{j} = (a τ_{j} + \frac{σ^{2}}{2}) and with \\ g_{i}^{j} = \{\begin{matrix} k exp (- α x_{i} - β τ_{j}), & x_{i} < {\bar{x}}_{f}^{j}, \\ 0, & otherwise, \end{matrix} \end{matrix}

where

{\bar{x}}_{f}^{j}

is the unique solution to the European put option problem

\bar{u} ({\bar{x}}_{f} (τ), τ) = f ({\bar{x}}_{f} (τ), τ),

(15)

where

\bar{u}

is a solution to the Cauchy problem:

{\bar{u}}_{τ} - {\bar{u}}_{x x} = 0

\bar{u} (x, 0) = f (x, 0)

, which is given in a closed explicit form. Including all the components, we obtain the

θ

-scheme discretized formulation for the penalised heat Equation (12)

A w^{j + 1} - B w^{j} - d^{j} = Δ τ \cdot p^{j},

where the multiplication of

Δ τ

results from the discretization.

4. Numerical Results

In this section, we consider the example for pricing American put options from Nielsen et al. (2002). All results were computed on an Intel^® Core™ i7-5557U CPU with 3.10 GHz. We chose

x_{\min} = - 4

x_{\max} = 4

M = 5000

, and used the parameter sets from Table 1. To facilitate the optimization, we summarised

Δ τ p^{j}

Δ τ p^{j} = (\bar{a} τ_{j} + Δ τ \cdot r) \cdot \{\begin{matrix} exp (- α x_{i} - β τ_{j}), & x_{i} < {\bar{x}}_{f}^{j}, \\ 0, & otherwise, \end{matrix}

(16)

where

\bar{a} = Δ τ \cdot k \cdot a

. Through the initial guess, the only unknown parameter is

\tilde{a}

. Since a deterministic expression for

\tilde{a}

is a goal of our future research, the penalty parameter

\tilde{a}

was obtained by an optimization. The optimization was done by minimizing the mean square error (MSE) of the solution corresponding to

\tilde{a}

. The MSE is given by

M S E = \frac{1}{N} \sum_{i = 1}^{N} {(P^{PSOR} (S_{i}, 0) - P^{Pen} (S_{i}, 0))}^{2}, S_{i} = K exp (x_{i}),

where

P^{PSOR}

is the solution obtained by the projected SOR algorithm and

P^{Pen}

is the maximum of the solution of the penalised system and the payoff function. The maximum was used to gain comparable results to the PSOR algorithm. Since the parameter set 2 is widely used in research, we compared the free boundary value of the parameter set 2 with the free boundary solution of Nielsen et al. (2002), Fazio et al. (2019), and Company et al. (2014). The value for the free boundary solution obtained by Nielson is 0.8622, Fazio obtained 0.86274, and the free boundary value of Company et al. (2014) is 0.8628, where our approach gives 0.86269 with a finer optimization and

\tilde{a} = 10.11 \times 10^{- 4}

. This comparison illustrates the high accuracy of this method in comparison with related research.

Since the PSOR method does not consider a penalty term, we also computed the MSE of the solution of the penalised PDE

u_{τ} - u_{x x} = \frac{1}{\hat{δ}} (f - u),

(17)

using a well-known penalty term (Günther and Jüngel 2010) with

\hat{δ} = 1 \times 10^{- 4}

and

N = 1000

. The results are presented in Table 2, and the comparison of the numerical results illustrate the accuracy of the method. As we included

Δ τ

\bar{a}

, the obtained values for

\bar{a}

are small, and the corresponding values for

\tilde{a}

are larger than 1. The best results were obtained by the sample sets with little volatility and short time maturity. The observation of the short time maturity is based on the fact that the number of points were different. The dependence on the volatility was caused by the simplification of the term p, since we cancelled out

σ^{2} / 2

and included

2 / σ^{2}

\tilde{a}

. We observed that the differences were in the range between an estimated free boundary value and the final free boundary value. They were caused by the time-dependent movement of the free boundary position. There are several ways to analyse this approach in detail. Since a sensitive point of the presented method is the choice of the initial free boundary value, an interesting approach for future research is a detailed analysis of the effect of the choice of the approximation formulas to the solution. As approximation formulas for the initial guess, one can choose the formulas in (Evans et al. 2002; Stamicar et al. 1999; Zhu 2006). Another idea is the consideration of an iterative scheme. In this idea, the free boundary value of the obtained solution is used as an initial guess for a second iteration of solving the penalised American Put problem.

5. Conclusions

The numerical results give clear evidence that using a non-constant penalty parameter

δ

is both feasible and beneficial. Future work will focus on the inclusion of the free boundary movement, a deterministic penalty function, and on the extension to multi-asset American options.

Author Contributions

All authors contributed equally to this research article. All authors have read and agree to the published version of the manuscript.

Funding

This research was funded by by the German-Slovakian Project ENANEFA–Efficient Numerical Approximation of Nonlinear Equations in Financial Applications financed by the DAAD and the Slovakian Ministry of Education.

Conflicts of Interest

The authors declare no conflict of interest.

References

Company, Rafael, Vera N. Egorova, and Lucas Jódar. 2014. Solving American option pricing models by the front-fixing method: Numerical analysis and computing. Abstract and Applied Analysis 2014: 146745. [Google Scholar] [CrossRef]
Cryer, Colin W. 1971. The solution of a quadratic programming problem using systematic overrelaxation. SIAM Journal on Control 9: 385–92. [Google Scholar] [CrossRef]
Evans, Jonathan D., Rachel Kuske, and Joseph Bishop Keller. 2002. American options on assets with dividends near expiry. Mathematical Finance 12: 219–37. [Google Scholar] [CrossRef]
Fazio, Riccardo, Alessandra Insana, and Alessandra Jannelli. 2019. Front-fixing finite difference schemes for American put options model. AIP Conference Proceedings 1738: 480123. [Google Scholar] [CrossRef]
Günther, Michael, and Ansgar Jüngel. 2010. Finanzderivate mit MATLAB: Mathematische Modellierung und Numerische Simulation, 2nd ed. Wiesbaden: Vieweg. [Google Scholar]
Han, Houde, and Xiaonan Wu. 2003. A fast numerical method for the Black-Scholes equation of American options. SIAM Journal on Numerical Analysis 41: 2081–95. [Google Scholar] [CrossRef]
Lauko, Martin, and Daniel Ševčovič. 2010. Comparison of numerical and analytical approximation of the early boundary of American put options. The ANZIAM Journal 51: 430–48. [Google Scholar] [CrossRef] [Green Version]
Nielsen, Bjørn Fredrik, Ola Skavhaug, and Aslak Tveito. 2002. Penalty and front-fixing methods for the numerical solution of American option problems. Journal of Computational Finance 5: 69–97. [Google Scholar] [CrossRef]
Schwartz, Eduardo S. 1977. The valuation of warrants: Implementing a new approach. Journal of Financial Economics 4: 79–93. [Google Scholar] [CrossRef]
Stamicar, Robert, Daniel Ševčovič, and John Chadam. 1999. The early exercise boundary for the American put near expiry: Numerical approximation. Canadian Applied Mathematics Quaterly 7: 427–44. [Google Scholar]
Wang, Song, Xiao-qi Yang, and Kok Lay Teo. 2006. Power penalty method for a linear complementarity problem arising from American option valuation. Journal of Optimization Theory and Applications 129: 227–54. [Google Scholar] [CrossRef]
Wilmott, Paul, Jeff Dewynne, and Sam Howison. 1993. Option Pricing: Mathematical Models and Computation. Oxford: Oxford Financial Press. [Google Scholar]
Zhu, Song-Ping. 2006. A new analytical approximation formula for the optimal exercise boundary of American put options. International Journal of Theoretical and Applied Finance 9: 1141–77. [Google Scholar] [CrossRef]

Table 1. Numerical results of the corresponding parameter sets.

Example	T	K	r	$σ$	N	$\bar{a}$ (×10⁻⁴)	$S_{f}$	MSE
1	3	100	0.08	0.2	1000	7.5	81.87	$9.6 \times 10^{- 3}$
1	3	100	0.08	0.2	2500	7.4	82.00	$6.2 \times 10^{- 3}$
2	1	1	0.1	0.2	1000	10.2	0.862	$5.2 \times 10^{- 5}$
2	1	1	0.1	0.2	2500	10.0	0.863	$3.3 \times 10^{- 5}$
3	0.05	10	0.1	0.25	1000	8.0	9.158	$3.5 \times 10^{- 5}$
3	0.05	10	0.1	0.25	2500	8.5	9.142	$3.5 \times 10^{- 5}$
4	0.1	100	0.1	0.3	1000	5.5	86.59	$1.2 \times 10^{- 3}$
4	0.1	100	0.1	0.3	2500	5.4	86.87	$7.6 \times 10^{- 4}$
5	1	100	0.1	0.4	1000	2.6	66.49	$1.4 \times 10^{- 2}$
5	1	100	0.1	0.4	2500	2.63	66.60	$9.2 \times 10^{- 3}$
6	0.05	50	0.1	0.4	1000	3.0	42.61	$3.8 \times 10^{- 4}$
6	0.05	50	0.1	0.4	2500	3.1	42.61	$2.5 \times 10^{- 4}$

Table 2. MSE for the parameter sets from the Table 1 computed with the penalty term (17).

1	2	3	4	5	6
$3.3 \times 10^{- 1}$	$1.8 \times 10^{- 3}$	$1.1 \times 10^{- 3}$	$2.2 \times 10^{- 2}$	$1.8 \times 10^{- 1}$	$5.5 \times 10^{- 4}$

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Clevenhaus, A.; Ehrhardt, M.; Günther, M.; Ševčovič, D. Pricing American Options with a Non-Constant Penalty Parameter. J. Risk Financial Manag. 2020, 13, 124. https://doi.org/10.3390/jrfm13060124

AMA Style

Clevenhaus A, Ehrhardt M, Günther M, Ševčovič D. Pricing American Options with a Non-Constant Penalty Parameter. Journal of Risk and Financial Management. 2020; 13(6):124. https://doi.org/10.3390/jrfm13060124

Chicago/Turabian Style

Clevenhaus, Anna, Matthias Ehrhardt, Michael Günther, and Daniel Ševčovič. 2020. "Pricing American Options with a Non-Constant Penalty Parameter" Journal of Risk and Financial Management 13, no. 6: 124. https://doi.org/10.3390/jrfm13060124

Article Menu

Pricing American Options with a Non-Constant Penalty Parameter

Abstract

1. Introduction

2. Mathematical Modeling

2.1. Transformation

2.2. The Penalty Term

3. Discretization

4. Numerical Results

5. Conclusions

Author Contributions

Funding

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI