  Harmonic and polyanalytic functional calculi have been recently defined for
bounded commuting operators. Their definitions are based on the Cauchy formula
of slice hyperholomorphic functions and on the factorization of the Laplace
operator in terms of the Cauchy-Fueter operator
  $\mathcal{D}$ and of its conjugate $\overline{\mathcal{D}}$. Thanks to the
Fueter extension theorem when we apply the operator $\mathcal{D}$ to slice
hyperholomorphic functions we obtain harmonic functions and via the Cauchy
formula of slice hyperholomorphic functions we establish an integral
representation for harmonic functions. This integral formula is used to define
the harmonic functional calculus on the $S$-spectrum. Another possibility is to
apply the conjugate of the Cauchy-Fueter operator to slice hyperholomorphic
functions. In this case, with a similar procedure we obtain the class of
polyanalytic functions, their integral representation and the associated
polyanalytic functional calculus. The aim of this paper is to extend the
harmonic and the polyanalytic functional calculi to the case of unbounded
operators and to prove some of the most important properties. These two
functional calculi belong to so called fine structures on the $S$-spectrum in
the quaternionic setting. Fine structures on the $S$-spectrum associated with
Clifford algebras constitute a new research area that deeply connects different
research fields such as operator theory, harmonic analysis and hypercomplex
analysis.
