Deep ReLU neural networks overcome the curse of dimensionality when approximating semilinear partial integro-differential equations

In this paper we consider PIDEs with gradient-independent Lipschitz continuous nonlinearities and prove that deep neural networks with ReLU activation function can approximate solutions of such semilinear PIDEs without curse of dimensionality in the sense that the required number of parameters in the deep neural networks increases at most polynomially in both the dimension $d$ of the corresponding PIDE and the reciprocal of the prescribed accuracy $\epsilon$.