Abstract

Parabolic partial differential equations (PDEs) are widely used in the mathematical modeling of natural phenomena and man made complex systems. In particular, parabolic PDEs are a fundamental tool to determine fair prices of financial derivatives in the financial industry. The PDEs appearing in financial engineering applications are often nonlinear and high dimensional since the dimension typically corresponds to the number of considered financial assets. A major issue is that most approximation methods for nonlinear PDEs in the literature suffer under the so-called curse of dimensionality in the sense that the computational effort to compute an approximation with a prescribed accuracy grows exponentially in the dimension of the PDE or in the reciprocal of the prescribed approximation accuracy and nearly all approximation methods have not been shown not to suffer under the curse of dimensionality. Recently, a new class of approximation schemes for semilinear parabolic PDEs, termed full history recursive multilevel Picard (MLP) algorithms, were introduced and it was proven that MLP algorithms do overcome the curse of dimensionality for semilinear heat equations. In this paper we extend those findings to a more general class of semilinear PDEs including as special cases semilinear Black-Scholes equations used for the pricing of financial derivatives with default risks. More specifically, we introduce an MLP algorithm for the approximation of solutions of semilinear Black-Scholes equations and prove that the computational effort of our method grows at most polynomially both in the dimension and the reciprocal of the prescribed approximation accuracy. This is, to the best of our knowledge, the first result showing that the approximation of solutions of semilinear Black-Scholes equations is a polynomially tractable approximation problem.

Highlights

  • Parabolic partial differential equations (PDEs) are key mathematical tools to model natural phenomena and man-made complex systems

  • Parabolic PDEs are used in the financial industry to model fair prices of financial derivatives

  • The use of PDEs for option pricing originated in the work of Black, Scholes, & Merton which suggested that the price of a financial derivative satisfies a linear parabolic PDE, nowadays known as Black-Scholes equation

Read more

Summary

Introduction

Parabolic partial differential equations (PDEs) are key mathematical tools to model natural phenomena and man-made complex systems. The conclusion of Theorem 1.1 above states that for every d ∈ N, ε ∈ (0, 1] there exists a natural number N ∈ N such that (i) it holds that VNd,,0N (0, ξd) approximates ud(0, ξd) in the strong L2-sense with accuracy ε and such that (ii) it holds that the computational effort to compute VNd,,0N (0, ξd) is essentially of the order d1+2(P+pq)ε−2 (cf (1.4) in Theorem 1.1 above for the precise formulation) This is exactly the computational complexity of standard Monte Carlo approximations of the solutions of the PDEs in (1.1) in the special case where the nonlinearity f vanishes (cf., e.g., Graham & Talay [51]).

Time-discrete Gronwall inequalities
A priori moment bounds for solutions of SDEs
Temporal regularity properties for solutions of SDEs
Strong error estimates for Euler-Maruyama approximations
On identically distributed random variables
On random evaluations of random fields
Brownian motions and right-continuous filtrations
On a distributional flow property for solutions of SDEs
Stochastic fixed point equations and MLP approximations
A priori bounds for solutions of stochastic fixed point equations
Properties of MLP approximations
Analysis of approximation errors of MLP approximations
Expectations of MLP approximations
Biases of MLP approximations
Estimates for the variances of MLP approximations
On a geometric time-discrete Gronwall inequality
Error estimates for MLP approximations
T qM n
Complexity analysis for MLP approximation algorithms
Viscosity solutions of semilinear PDEs and stochastic fixed point equations
MLP approximations in fixed space dimensions
MLP approximations in variable space dimensions
MLP approximations for PDE models
MLP approximations for semilinear heat equations
MLP approximations for semilinear Black-Scholes equations
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call