Front-fixing FEMs for the pricing of American options based on a PML technique

In this paper, efficient numerical methods are developed for the pricing of American options governed by the Black–Scholes equation. The front-fixing technique is first employed to transform the free boundary of optimal exercise prices to some a priori known temporal line for a one-dimensional parabolic problem via the change of variables. The perfectly matched layer (PML) technique is then applied to the pricing problem for the effective truncation of the semi-infinite domain. Finite element methods using the respective continuous and discontinuous Galerkin discretization are proposed for the resulting truncated PML problems related to the options and Greeks. The free boundary is determined by Newton’s method coupled with the discrete truncated PML problem. Stability and nonnegativeness are established for the approximate solution to the truncated PML problem. Under some weak assumptions on the PML medium parameters, it is also proved that the solution of the truncated PML problem converges to that of the unbounded Black–Scholes equation in the computational domain and decays exponentially in the perfectly matched layer. Numerical experiments are conducted to test the performance of the proposed methods and to compare them with some existing methods.

[1]  G. Burton.  Sobolev Spaces , 2013 .

[2]  L. Rogers.  Monte Carlo valuation of American options , 2002 .

[3]  Michael V. Klibanov,et al.  Global uniqueness and Hölder stability for recovering a nonlinear source term in a parabolic equation , 2005 .

[4]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[5]  P. Glasserman,et al.  Monte Carlo methods for security pricing , 1997 .

[6]  Ravi P. Agarwal,et al.  The One-Dimensional Heat Equation , 2009 .

[7]  Jing Zhao,et al.  An Artificial Boundary Method for American Option Pricing under the CEV Model , 2008, SIAM J. Numer. Anal..

[8]  Felix E. Browder,et al.  The One-Dimensional Heat Equation: Preface , 1984 .

[9]  Alexander V. Tikhonravov,et al.  Estimates of initial conditions of parabolic equations and inequalities in infinite domains via lateral Cauchy data , 2007 .

[10]  Frédéric Nataf,et al.  Perfectly matched layers for the heat and advection-diffusion equations , 2010, J. Comput. Phys..

[11]  H. G. Landau,et al.  Heat conduction in a melting solid , 1950 .

[12]  Robert L. Sproull,et al.  The Conduction of Heat in Solids , 1962 .

[13]  J. Hull.  Fundamentals of Futures and Options Markets , 2001 .

[14]  Bernardo Cockburn.  Discontinuous Galerkin methods , 2003 .

[15]  D. Lamberton,et al.  Variational inequalities and the pricing of American options , 1990 .

[16]  Xiaonan Wu,et al.  A Fast Numerical Method for the Black-Scholes Equation of American Options , 2003, SIAM J. Numer. Anal..

[17]  Walter Allegretto,et al.  Finite Element Error Estimates for a Nonlocal Problem in American Option Valuation , 2001, SIAM J. Numer. Anal..

[18]  G. Andrews,et al.  Encyclopedia of Mathematics and its Applications , 1990 .

[19]  Joseph E. Pasciak,et al.  Analysis of a finite PML approximation to the three dimensional elastic wave scattering problem , 2010, Math. Comput..

[20]  Béatrice Rivière,et al.  Discontinuous Galerkin methods for solving elliptic and parabolic equations - theory and implementation , 2008, Frontiers in applied mathematics.

[21]  Haijun Wu,et al.  An Adaptive Finite Element Method with Perfectly Matched Absorbing Layers for the Wave Scattering by Periodic Structures , 2003, SIAM J. Numer. Anal..

[22]  Francis A. Longstaff,et al.  Valuing American Options by Simulation: A Simple Least-Squares Approach , 2001 .

[23]  Chi-Wang Shu,et al.  Discontinuous Galerkin Methods: Theory, Computation and Applications , 2011 .

[24]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..

[25]  Eduardo S. Schwartz,et al.  Finite Difference Methods and Jump Processes Arising in the Pricing of Contingent Claims: A Synthesis , 1977 .

[26]  RAUL KANGRO,et al.  Far Field Boundary Conditions for Black-Scholes Equations , 2000, SIAM J. Numer. Anal..

[27]  A. Ern,et al.  Mathematical Aspects of Discontinuous Galerkin Methods , 2011 .

[28]  Gang Bao,et al.  Convergence Analysis of the Perfectly Matched Layer Problemsfor Time-Harmonic Maxwell's Equations , 2005, SIAM J. Numer. Anal..

[29]  A. Friedman.  Partial Differential Equations of Parabolic Type , 1983 .

[30]  S. Ross,et al.  Option pricing: A simplified approach☆ , 1979 .

[31]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[32]  Lishang Jiang.  Mathematical Modeling and Methods of Option Pricing , 2005 .

[33]  Hongtao Yang,et al.  A Front-Fixing Finite Element Method for the Valuation of American Options , 2008, SIAM J. Sci. Comput..

[34]  Y. Kwok.  Mathematical models of financial derivatives , 2008 .

[35]  J. Bérenger.  Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves , 1996 .

[36]  Christoph Walker,et al.  Elliptic and Parabolic Equations , 2015 .

[37]  Kai Zhang,et al.  Studies on some perfectly matched layers for one-dimensional time-dependent systems , 2009, Adv. Comput. Math..

[38]  Kaushik I. Amin,et al.  CONVERGENCE OF AMERICAN OPTION VALUES FROM DISCRETE‐ TO CONTINUOUS‐TIME FINANCIAL MODELS1 , 1994 .