Abstract:Machine learning for partial differential equations (PDEs) is a hot topic. In this paper we introduce and analyse a Deep BSDE scheme for nonlinear integro-PDEs with unbounded nonlocal operators -problems arising in e.g. stochastic control and games involving infinite activity jump-processes. The scheme is based on a stochastic forward-backward SDE representation of the solution of the PDE and (i) approximation of small jumps by a Gaussian process, (ii) simulation of the forward part, and (iii) a neural net regression for the backward part. Unlike grid-based schemes, it does not suffer from the curse of dimensionality and is therefore suitable for high dimensional problems. The scheme is designed to be convergent even in the infinite activity/unbounded nonlocal operator case. A full convergence analysis is given and constitutes the main part of the paper.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
The paper aims to solve the numerical solution problem of a class of nonlinear integro - partial differential equations (integro - PDEs) with unbounded non - local operators. Specifically, the paper proposes a method based on deep backward stochastic differential equations (Deep BSDE) to deal with such equations. These equations are widely used in stochastic control, game theory, and financial models involving infinite - activity jump processes.
### Specific problem description
Consider the nonlinear integro - partial differential equation in the following form:
\[
\begin{cases}
\frac{\partial u}{\partial t}(t,x)+\tilde{L}[u](t,x)+f(t,x,u,\sigma\nabla u,B[u]) = 0, &\text{for }(t,x)\in[0,T)\times\mathbb{R}^q\\
u(T,x)=g(x), &\text{for }x\in\mathbb{R}^q
\end{cases}
\]
where $\tilde{L}$ is a degenerate elliptic integro - differential operator, defined as:
\[
\tilde{L}u=b\cdot D_xu+\frac{1}{2}\text{Tr}(\sigma\sigma^T D_x^2u)+J[u]
\]
and $J$ and $B$ are possibly unbounded fractional/non - local operators, defined as:
\[
J[u](t,x)=\int_{E}\left[u(t,x +\beta(x,e))-u(t,x)-D_xu(t,x)\beta(x,e)\right]\nu(de)
\]
\[
B[u](t,x)=\int_{E}\left[u(t,x+\beta(x,e))-u(t,x)\right]\gamma(e)\nu(de)
\]
Here, $E=\mathbb{R}^q\setminus\{0\}$, $\nu$ is a Levy measure, and $\beta$ and $\gamma$ are appropriate functions.
### Main challenges
1. **High - dimensional problems**: Traditional grid methods (such as the finite difference method, the semi - Lagrangian method, and the finite element method) become infeasible for high - dimensional problems due to the curse of dimensionality.
2. **Unbounded non - local operators**: The non - local operators in the equations may be unbounded, which increases the difficulty of numerical solution.
3. **Infinite - activity jump processes**: The Poisson noise driving the equations may have infinite activity (i.e., there are an infinite number of jumps in each time interval), which requires special methods to handle.
### Solutions
The paper proposes a method based on deep backward stochastic differential equations (Deep BSDE) to solve the above problems through the following steps:
1. **Gaussian approximation of small jumps**: Deal with the infinite - activity jump processes by approximating small jumps as Gaussian processes.
2. **Simulation of the forward part**: Simulate the forward SDE part using an Euler - Maruyama - type discretization method.
3. **Neural network regression of the backward part**: Calculate the backward part through deep neural network regression with local - time optimization.
### Convergence analysis
The paper also provides a complete convergence analysis of this method, including weak and strong convergence rates. These results not only verify the effectiveness of the method but also provide theoretical support for practical applications.
In summary, the paper successfully solves the high - dimensional numerical solution problem of nonlinear integro - partial differential equations with unbounded non - local operators by introducing a new Deep BSDE method.