Lattice Boltzmann methods for solving partial differential equations of exotic option pricing

Abstract

This paper establishes a lattice Boltzmann method (LBM) with two amending functions for solving partial differential equations (PDEs) arising in Asian and lookback options pricing. The time evolution of stock prices can be regarded as the movement of randomizing particles in different directions, and the discrete scheme of LBM can be interpreted as the binomial models. With the Chapman-Enskog multi-scale expansion, the PDEs are recovered correctly from the continuous Boltzmann equation and the computational complexity is O(N), where N is the number of space nodes. Compared to the traditional LBM, the coefficients of equilibrium distribution and amending functions are taken as polynomials instead of constants. The stability of LBM is studied via numerical examples and numerical comparisons show that the LBM is as accurate as the existing numerical methods for pricing the exotic options and takes much less CPU time.

Key words： Exotic option pricing; lattice Boltzmann method; Chapman-Enskog multi-scale expansion; stability; computational complexity

References

1 Al-Zoubi A, Brenner G. Comparative study of thermal flows with different finite volume and lattice Boltzmann schemes. Int J Mod Phys C, 2004, 15: 307–319

2 Bhantnagar J, Gross E P, Krook M K. A model for collision processes in gas I: Small amplitude processes in charged and neutral one-componet systems.Phys Rev, 1954, 94: 511–525

3 Black F, Scholes M. The pricing of options and corporate liabilities. J Polit Econom, 1973, 81: 637–654

4 Boyle P, Draviam T. Pricing exotic options under regime switching. Insurance: Math Econom, 2007, 40: 267–282

5 Broadie M, Detemple J. Option pricing: valuation models and applications. Management Sci, 2004, 50: 1145–1177

6 Chen L, Ma C F. A lattice Boltzmann model with an amending function for simulating nonlinear partial differential equations. Chin Phys B, 2010, 19: 1–8

7 Dubois F, Leli`evre T. Efficient pricing of Asian options by the PDE approach. J Comput Finan, 2005, 8: 55–64

8 Geman H, Yor M. Bessel processes Asian options and perpetuities. Math Finan, 1993, 3: 349–375

9 Goldman M B, Sosin H B, Gatto M A. Path-dependent options buy at the low, sell at the high. J Finan, 1979, 34: 1111–1127

10 Guo Z L, Zheng C G, Shi B C. Non-equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method. Chin Phys, 2002, 11: 366–374

11 Ingersoll J. Theory of Financial Decision Making. Totowa/New Jersey: Roman & Littlefield, 1987

12 Lai H L, Ma C F. A higher order lattice BGK model for simulating some nonlinear partial differential equations. Sci China Ser G, 2009, 52: 1053–1061

13 Lai H L, Ma C F. Lattice Boltzmann method for the generalized Kuramoto-Sivashinsky equation. Phys A, 2009, 388: 1405–1412

14 Lai H L, Ma C F. Lattice Boltzmann model for generalized nonlinear wave equations. Phys Rev E, 2011, 84: 046708

15 15. Lai H L, Ma C F. A new lattice Boltzmann model for solving the coupled viscous Burgers’ equation. Phys A, 2014, 395: 445–457

16 16. Merton R. Theory of rational option pricing. Bell J Econom Management Sci, 1973, 4: 141–183

17 17. Rogers L C G, Shi Z. The value of an Asian option. J Appl Probab, 1995, 32: 1077–1088

18 18. Shreve S E. Stochastic Calculus for Finance: the Binomial Asset Pricing Model. New York: Springer, 2003

19 19. Succi S. The Lattice Boltzmann Equation for Fluid Dynamics and Beyond. Oxford: Oxford University Press, 2001

20 20. Sukop M C, Thorne D T. Lattice Boltzmann Modeling: An Introduction for Geoscientists and Engineers. Heidelberg/Berlin/New York: Springer, 2006

21 21. Vecer J. A new PDE approach for pricing arithmetic average Asian options. J Comput Finan, 2001, 4: 105–113

22 22. Wilmott P, Dewynne J, Howison S. Option Pricing: Mathematical Models and Computation. Oxford: Oxford Financial Press, 1997