We consider a discretization of the Laplace operator on $\mathbb Z^2$, <https://en.wikipedia.org/wiki/Discrete_Laplace_operator>


Then, it is natural to consider its fundamental solution $u$, i.e. $|u(x)|\leq C \ln|x|,\Delta u = 1$ at $(0,0)$ and $\Delta u=0$ elsewhere. I am sure that somewhere it is proven that the discrete derivatives of $u$ behave exactly as we expect them : $|u'(x)|\leq \frac{1}{|x|}, |u''(x)|\leq c\frac{1}{|x|^2}$.
But I can not find anything like that in the Internet. Could you provide me with a reference, please? Any information about $c$?



Fundamental solution of Discrete Laplace in the plane
For nearest neighbor Laplacian, you can find a short self-contained proof of the formula $$u(x)=\frac1{2\pi}\log |x|+c+O\left(\frac1{|x|^2}\right)$$ at [these lecture notes](https://wiki.helsinki.fi/download/attachments/155233241/LectureNotes20_04.pdf?version=2&modificationDate=1430322536394&api=v2), pages 6-7. It is a simple exercise in Laplace's method to work out further coefficients of the asymptotic expansion, if necessary. The proof there is not the shortest possible; the easiest way would be to compare the double integral formula defining $G\_0$ on page 6 with its counterpart for continuous Laplacian. But then you do not get the beautiful formula (4.4) along the way.


"Random Walk: a modern introduction" by Lawler and Limic does the same in much greater generality (not necessarily nearest neighbour), using local CLT with error bounds.


Historically, the first sources usually cited for the above formula are


Stohr, A., [50]: "Uber einige lineare partielle Differenzengleichungen mit konstanten Koefizienten," I, II, and Ill, Mathematische Nachrichten 3 (1950),


or even


W. H. McCrea and F. J. W. Whipple, Random paths in two and three dimensions, Proc. Roy. Soc. Edinburgh 60 (1940),


although the latter apparently proves a slightly weaker result (I have no access to it at the moment to check).


