It is well known that the fundamental solution of
$${u_t}left( {n,t} right) = uleft( {n + 1,t} right) - 2uleft( {n,t} right) + uleft( {n - 1,t} right),n in mathbb{Z},$$
with
u(
n, 0) =
δ nm for every fixed
m ∈ Z is given by
u(
n,
t) =
e ?2t I n?m (2
t), where
I k (
t) is the Bessel function of imaginary argument. In other words, the heat semigroup of the discrete Laplacian is described by the formal series
W t f(
n) = Σ
m∈Z e ?2t I n?m (2
t)
f(
m). This formula allows us to analyze some operators associated with the discrete Laplacian using semigroup theory. In particular, we obtain the maximum principle for the discrete fractional Laplacian, weighted ?
p (Z)-boundedness of conjugate harmonic functions, Riesz transforms and square functions of Littlewood-Paley. We also show that the Riesz transforms essentially coincide with the so-called discrete Hilbert transform defined by D. Hilbert at the beginning of the twentieth century. We also see that these Riesz transforms are limits of the conjugate harmonic functions. The results rely on a careful use of several properties of Bessel functions.