Modulus of continuity for solutions of non-local heat equations

We extend the method of modulus of continuity for solutions of parabolic equations--as used, for instance, to prove the Fundamental Gap Conjecture--to solutions of non-local heat equations on R^n and in dimension one with a non-local Neumann boundary condition. Specifically, we show that if a solution of a non-local heat equation has an initial modulus of continuity satisfying simple criteria, then this modulus of continuity is preserved at all subsequent times. In the process of trying to generalise our result in one dimension, we found a counterexample suggesting that a non-local analogue of the Payne-Weinberger inequality would depend on more than the diameter of a bounded (convex) domain.