Given consider the second order differential operator in $(0,\infty)$ $L_\alpha f \equiv \left(x^2\, \frac{d^2}{dx^2} + (2\alpha+3)x\, \frac{d}{dx} + x^2 + (\alpha+1)^2 \right)(f),$ which appears in the theory of Bessel functions. The purpose of this paper is to develop the corresponding harmonic analysis taking as the analogue to the classical Laplacian. Namely we study the boundedness properties of the heat and Poisson semigroups. These boundedness properties allow us to obtain some convergence results that can be used to solve the Cauchy problem for the corresponding heat and Poisson equations.