“…In recent years, there has been great interest in using nonlocal integro-differential equations (IDEs) as a means to describe physical systems, due to their natural ability to describe physical phenomena at small scales and their reduced regularity requirements which lead to greater flexibility [54,8,64,29,30,63,27,25,49,42,41,13,22,18,37,24,3,17,52,14,31]. In particular, nonlocal problems with Neumann-type boundary constraints have received particular attention [15,16,20,34,47,7,21,23,26,55,51,1,46,63] due to their prevalence in describing problems related to: interfaces [2], free boundaries, and multiscale/multiphysics coupling problems [38,53,62,5,6]. Unlike classical PDE models, in the nonlocal IDEs the boundary conditions must be defined on a region with non-zero volume outside the surface [16,21,55], in contrast to more traditional engineering scenarios where boundary conditio...…”