Deep ReLU neural networks overcome the curse of dimensionality in the numerical approximation of semilinear partial integro-differential equations

We prove that deep neural networks with ReLU activation function are capable of approximating solutions of semilinear partial integro-differential equations in the case of gradient-independent and Lipschitz-continuous nonlinearities, while the required number of parameters in the neural networks grows at most polynomially in both the dimension $ d\in\mathbb{N} $ and the reciprocal of the prescribed accuracy $ \epsilon $.