Harmonic functions for a class of integro-differential operators

We consider the operator $\text{sL}$ defined on $C^2({\text{bR}}^d)$ functions by \sL f(x)&=&{1/2}\sum_{i,j=1}^d a_{ij}(x)\frac{\partial^2f(x)}{\partial x_i\partial x_j}+\sum_{i=1}^d b_i(x)\frac{\partial f(x)}{\partial x_i} &+&\int_{\bR^d\backslash{0}}[f(x+h)-f(x)-1_{(|h|\leq1)}h\cdot \grad f(x)]n(x,h)dh. Under the assumption that the local part of the operator is uniformly elliptic and with suitable conditions on $n(x,h)$, we establish a Harnack inequality for functions that are nonnegative in ${\text{bR}}^d$ and harmonic in a domain. We also show that the Harnack inequality can fail without suitable conditions on $n(x,h)$. A regularity theorem for those nonnegative harmonic functions is also proved