\section{Introduction}
In \cite{ace15} the authors presented a matrix representation
approach to all types of Appell polynomial sequences $\left\{p_k
(x)\right\}_{k\geq0}$ of one real variable which relies
essentially on the matrix $H$ defined by
\begin{equation} \label{H}
(H)_{ij}=\left\{
    \begin{array}{ll}
      i, & \quad  i=j+1 \\
      0, & \quad  \hbox{otherwise,}  \qquad i,j=0,1,\ldots,m.
    \end{array}
  \right.
\end{equation}
In fact,  considering that a sequence of real polynomials  of degree $k$
$\{p_k(x)\}_{k\geq 0}$ is called \textit{Appell polynomial sequence} if
\begin{equation}   \label{Apdiff1}
\displaystyle{\frac{d}{dx}}p_k(x) = k \, p_{k-1}(x), \quad
k=1,2,\dots,
\end{equation}
cf.  \cite{app80}, by introducing the vector
${\bf p}(x)=[p_0(x)\;\;p_1(x)\;\cdots\;p_m(x)]^T,$
from relation (\ref{Apdiff1}) we obtain the  first order differential equation
\begin{eqnarray} \label{ode}
\frac{d}{dx}{\bf p}(x) = H \, {\bf p}(x),
\end{eqnarray}
whose general solution is
\begin{eqnarray}   \label{genPasc}
{\bf p}(x) = e^{Hx} \, {\bf p}(0).
\end{eqnarray}
It is evident that the role of $H$ is twofold: it acts both as
\textit{derivation matrix} and  also as some type of \textit{creation matrix};
in fact, the different kinds of Appell polynomials are uniquely determined
from (\ref{genPasc})  by choosing the entries of the initial value vector
${\bf p}(0).$ In some sense, this matrix reveals the arithmetical
tape which glues the different types of Appell polynomials together.

The prototype of an Appell sequence are the monomials $p_k(x)=x^k.$ The
Bernoulli
polynomials as the most prominent representatives of Appell polynomials are in
Numerical Analysis almost so important as orthogonal polynomials
are \cite{cos99,cos06,milo13}. But the list of Appell polynomials includes,
among others, also the classical polynomials named after Frobenius-Euler,
Hermite, Laguerre, and Chebyshev.

Although an old subject, during the last two decades the
interest in Appell polynomials and their applications has
significantly increased.
As very few examples, recent applications of Appell polynomials exist
in fields like probability theory and statistics, cf.
\cite{ans09, sal11}, linear elasticity \cite{bock} or approximation of
3D-mappings in \cite{mafa10a}. Results in the
framework of noncommutative Clifford algebras and related to the
Gel'fand-Tsetlin branching approach gave evidence to Appell
polynomial sequences with shift, cf. \cite{pe13}, as
sequences of orthogonal polynomials in several variables, cf.
\cite{brasls10, lav12}. Operational approaches based on Appell
polynomials for generalizing Jacobi, Laguerre, Gould-Hopper,
and Chebyshev polynomials are used in the recent papers
\cite{cael14,cafama11a, cama11}. Employing methods of
representation theory, they are also tools for applications in
quantum physics \cite{wein95}.

Some authors were concerned with
finding new characterizations of Appell polynomials themselves through new
approaches. We mention, for instance,
the approach developed in \cite{yan09}, which makes
use of the generalized Pascal functional matrices and
the characterization proposed in  \cite{cos10}
which is based on a determinantal definition.
As previously quoted, the authors introduced in \cite{ace15} a matrix
approach in the real case which has already been applied in the context of
``image synthesis'' \cite{rom14}, and in connection with ``evolution equations''
 for the construction of a new algorithm of deflation \cite{yam15}.

In addition,  Appell polynomial sequences were subject to innumerable
generalizations, mostly depending from the type of applications where they could
be advantageously used, including multivariate
commutative or noncommutative settings.

This paper intends to confirm the advantageously use of the
matrix approach to real Appell polynomial sequences developed
in \cite{ace15} for the case of monogenic polynomials in
arbitrary dimensions and in a noncommutative hypercomplex
setting, i.e. for Clifford algebra valued polynomials in the
kernel of a generalized Cauchy-Riemann operator, cf.
\cite{desosu92}. In \cite{guema99} it has been shown that
monogenic functions are hypercomplex differentiable. Therefore
the classical definition in (\ref{Apdiff1}) by the derivative
property can be analogously used for defining monogenic Appell
sequences with respect to the hypercomplex derivative, cf.
\cite{bogue10, fama06, fama07, mafa07}.

The paper is organized as follows. Basic concepts and notations used in Hypercomplex
Function Theory are given in Section 2.  In Section 3  the
extension to the hypercomplex case of the unifying matrix approach developed in
\cite{ace15} is introduced. Since the
hypercomplex generalized Cauchy-Riemann operator is a linear
combination of two different types of first order differential
operators (a scalar one and a vector operator in the sense of
the underlying Clifford algebra) also two different nilpotent
matrices $H$ and $\tilde{H}$ are involved. Their relationship is essential
for guaranteeing the desired Appell property. The relation between
$H$ and $\tilde{H}$ reveals in a new and condensed matrix form
the transition from the real to the hypercomplex case.
Applying the results of Section 3 and the concept of a transfer
matrix, in Section 4 some of the main hypercomplex monogenic
counterparts of real Appell sequences are listed.

\section{Basic concepts and notation}

For an independent reading, in this section we repeat briefly some notions and
results useful in the sequel,  mainly following
\cite{cafama12,desosu92, GHSp08, ma04}.
\begin{definition}
Let $\{e_1,e_2,\dots,e_n\}$ be an orthonormal basis of the
Euclidean vector space $\mathbb{R}^{n}$ provided with a non-commutative
product according to the  multiplication rules
$e_i e_j +e_j e_i =-2\delta_{ij},\,i,j=1,2,\dots,n,$
where $\delta_{ij}$ is the Kronecker symbol. The associative
$2^n-$dimensional Clifford algebra ${\mathcal{C} \ell}_{0,n}$ over $\mathbb{R}$ is the set
of numbers $\alpha \in {\mathcal{C} \ell}_{0,n}$ of the form
$\sum_{A}\alpha_{A} e_{A},$
 with the basis $\{e_A: A\subseteq\; \{ 1,\dots,n\}\}$
 formed by $e_A=e_{h_1}e_{h_2}\dots e_{h_r}$, $1\leq h_1< \dots < h_r\leq n,  e_{\emptyset}=e_0=1\;$ and where
the components $\alpha_A$ are real numbers.
 The conjugate of $\alpha$ is defined
 by $\bar{\alpha}=\sum_A \alpha_A\bar{e}_A,$ with
$\bar{e}_A=\bar{e}_{h_r}\bar{e}_{h_{r-1}}\dots \bar{e}_{h_1};\;
\bar{e}_k=-e_k,
k=1,\ldots, n, \bar{e}_0=e_0=1.$
\end{definition}

In general, the vector space $\mathbb{R}^{n+1}$ is embedded
in ${\mathcal{C} \ell}_{0,n}$ by identifying the element $(x_0,x_1,\dots,x_n)\in
\mathbb{R}^{n+1}$ with an element of the real vector space $\mathcal{A}_n:=\hbox{span}_{\mathbb{R}}\{1,e_1,\ldots,e_n\}\subset {\mathcal{C} \ell}_{0,n}.$
For our purpose we only consider such elements $x \in
\mathcal{A}_n$  of the form $x=x_0+\sum_{k=1}^ne_kx_k=x_0+\underline{x}$, called
\emph{paravectors } (naturally, $x_0$ and $\underline{x}$ are called the
\emph{scalar part} and the \emph{vector part} of $x$,
respectively). Similarly to the complex case the conjugate
$\bar{x}$ and the norm $|x|$ of $x$ are given by
$\bar{x}=x_0-\underline{x}$ and
$|x|=(x\bar{x})^{1/2}=(\bar{x}x)^{1/2}=(x_0^2+x_1^2+\dots+x_n^2)^{1/2},$
respectively.

In Hypercomplex Function Theory $\, {\mathcal{C} \ell}_{0,n}-$valued functions are studied. They are
functions $f: \Omega \subseteq \mathbb{R}^{n+1}\cong \mathcal{A}_n \rightarrow {\mathcal{C} \ell}_{0,n}$ defined in an
open subset by $f(z)=\sum_A
f_A(z)e_A,$ where $f_A(z)$ are real valued functions. We will focus on a special
class of these functions analogous to complex holomorphic functions and
connected with them via the following concept.
\begin{definition}
A function $f$ is called left (right) monogenic in  $\Omega$ if
it is a solution of the differential equation
$\overline{\partial}f=0$ ($f\overline{\partial}=0$) where
$$\overline{\partial}:=\frac{1}{2}(\partial_0+\partial_{\underline{x}}), \;\;\hbox{with}\;\;
\partial_0:=\frac{\partial}{\partial {x_0}},\;\;\hbox{and}\;\; \partial_{\underline{x}}:=\sum_{k=1}^ne_k\frac{\partial}{\partial x_k},$$
generalizes the Cauchy-Riemann operator ($n=1, e_1\equiv i, z=x+iy$)
\[\frac{\partial}{\partial {\bar{z}}}=\frac{1}{2}\left(\frac{\partial}{\partial
{x}}+i \frac{\partial}{\partial {y}}\right).\]
The operator $\partial:=\frac{1}{2}(\partial_0-\partial_{\underline{x}})$ is called the
conjugate generalized Cauchy-Riemann operator or the hypercomplex differential
operator.
\end{definition}
Hereafter we only deal with left monogenic functions and we shall refer to them
simply as monogenic functions (right monogenic functions are treated
analogously). Notice that in the case of a paravector-valued function $f$ of the
variable $x \in \mathcal{A}_n,$ the hypercomplex partial differential equation
$\overline{\partial}f=0$  is, except for the
complex case ($n=1$), equivalent to an over-determined Cauchy-Riemann system of
$\frac{n(n+1)+2}{2}$ first order differential equations for the $(n+1)$
component functions of $f.$

We remark that hypercomplex differentiability as generalization of complex
differentiability has to be understood in the following way: a
function $f$ defined in an open domain $\Omega \subseteq \mathbb{R}^{n+1}$ is
hypercomplex differentiable supposed there exists in each point of $\Omega$ a
uniquely defined areolar derivative $f'.$ Then $f$ is real differentiable and
$f':=\partial f$. On the other hand,  $f$ is hypercomplex
differentiable in $\Omega$ if and only if $f$ is monogenic, cf.
\cite{guema99}. Consequently, if a hypercomplex function is monogenic then the
existence of the hypercomplex derivative is guaranteed and can be obtained as
result of the application of the conjugate generalized Cauchy-Riemann operator.

\section{Sequences of monogenic hypercomplex Appell polynomials}

Considering that for each $x \in \mathcal{A}_n$,
\[\overline{\partial}x^n=\frac{1}{2}(1-n), \qquad n \in \mathbb{N},\]
only in the complex case ($n=1$) the function $x^n$ belongs to the set of
monogenic functions. Obviously, this fact causes
problems for the consideration of monogenic polynomials in the
ordinary way and, moreover, for the whole understanding of a
suitable analog to power series in monogenic function theory,
cf.  \cite{ma90b}. Indeed, in some sense the problem of
embedding integer powers of $x$ in a theory of monogenic
functions was in the 90-ties of the last century for more than
a decade a driving force for modifying or extending the class
of monogenic functions. The result was, for example, the
consideration of a \textit{Modified Clifford Analysis} (H.
Leutwiler et al.) or the introduction of \textit{holomorphic
Cliffordian functions} (G. Laville et al.). The third way to overcome the
problem `inside' of the class of ordinary monogenic
functions, namely by generalizing Appell's concept of
power-like polynomials will be explained now (for further details,
see \cite{cafama12}).

Motivated by relation  (\ref{Apdiff1}) and the previously mentioned concept of
hypercomplex derivative, Appell sequences of homogeneous
monogenic
polynomials in the framework of Clifford Analysis have been
introduced in the following way, \cite{fama07, mafa07}:

\begin{definition}\label{generalizedAS}
A sequence of multivariate  homogeneous polynomials $\{\phi_k(x)\}_{k\geq 0}$ of
degree $k$ in the variable $x \in \mathcal{A}_n$  is called a
\it{generalized Appell sequence} with respect to  the
hypercomplex differential operator  $\partial$  if the following conditions
are satisfied:
\begin{description}
\item[(i)] $\overline{\partial}\phi_k(x)=0,$ i.e.,  $\phi_k(x)$ is monogenic
for each $\, k \ge 0;$
\item[(ii)] $\partial \phi_k(x)=k\phi_{k-1}(x),\quad    k=1,2,\dots.$
\end{description}
\end{definition}

 \begin{oss}
 Hereafter, if $x=x_0+ \underline{x} \in \mathcal{A}_n,$ we shall refer to $\phi_k(x)$ by using
equivalently the following notations: $\phi_k(x_0+\underline{x})$ or $\phi_k(x_0,
\underline{x}).$
 \end{oss}

Before generalizing the matrix approach to hypercomplex Appell
polynomials, it is important to observe that, setting $p_0(x) \equiv c_0 \neq
0,$  from  (\ref{Apdiff1}) the following explicit representation  for the
truncated sequence $\{p_k(x)\}_{k=0}^m$ of Appell polynomials in one real
variable occurs:
\begin{eqnarray*}
p_0(x) &=& c_0    \nonumber\\
p_1(x)&=& c_1+c_0 \, x   \nonumber \\
p_2(x)&=& c_2+2 \,c_1 \,x +c_0 \, x^2  \nonumber  \\
&\vdots& \nonumber\\
p_m(x)&=& c_m+ {m\choose 1}c_{m-1}\, x+{m\choose 2}c_{m-2}\, x^2+\cdots
+{m\choose m}c_{0}\, x^m,
\end{eqnarray*}
or, equivalently, in a more  compact form
\begin{equation}\label{compacAp}
p_k(x)= \sum_{j=0}^k \left(
                       \begin{array}{c}
                         k \\
                         j \\
                       \end{array}
                     \right) c_j \, x^{k-j},  \quad k=0,1,\dots,m, \qquad  c_0
\neq 0.
\end{equation}
In particular, in this notation
\begin{equation} \label{p0}
 {\bf p}(0) = [c_0  \; c_1 \; \ldots\; c_m]^T.
\end{equation}
The formal replacement in  (\ref{compacAp}) and (\ref{p0}) of the real
variable $x$ by $x_0,$ the scalar part of the paravector $x=x_0 + \underline{x} \in \mathcal{A}_n,$ as well
as that of $c_j$ by  $\tilde{c}_j \,  \underline{x}^j$ leads to  multivariate  homogeneous polynomials
 \begin{equation}  \label{phipar}
 \phi_k(x) =  \sum_{j=0}^k \left(
                       \begin{array}{c}
                         k \\
                         j \\
                       \end{array}
                     \right) \tilde{c}_j \,  \underline{x}^j \, x_0^{k-j},  \quad k=0,1,\dots,m
 \end{equation}
 and to the vector $ {\bm \phi} (0,\underline{x}) = [\tilde{c}_0 \,  \underline{x}^0  \;\; \tilde{c}_1 \,  \underline{x}^1 \; \;\ldots\; \;\tilde{c}_m \,  \underline{x}^m]^T.$
Consequently,  recalling that the entries of the lower triangular
 \textit{generalized Pascal matrix} $P(x_0) \equiv e^{H x_0}$ are given by
  \begin{eqnarray}   \label{Px}
 (P(x_0))_{ij} = \left\{
     \begin{array}{cl}
     {i \choose j} \, x_0^{i-j}, & \quad  \textrm{$i\geq j $} \\
       0, & \quad  \hbox{otherwise,}  \qquad i,j=0,1,\ldots,m,
     \end{array}
   \right.
 \end{eqnarray}
and denoting by  $ {{\bm \phi}}(x) =  [\phi_0(x) \;\;\phi_1(x) \;\;\cdots\;\;\phi_m(x)]^T,$
the corresponding matrix form of  (\ref{phipar}) is
\[
 {\bm \phi}(x) = e^{H x_0}   {\bm \phi} (0,\underline{x}),
 \]
which represents the counterpart in the hypercomplex framework of (\ref{genPasc}).
 Introducing the diagonal matrix  $D_{\tilde{c}}=\hbox{diag}[\tilde{c}_0\;\;\tilde{c}_1\;\cdots \;\tilde{c}_m],$ and the vector
 \begin{equation}  \label{xv}
{\bm \xi} (\underline{x})= [1\;\;\underline{x}\;\;\underline{x}^2\;\;\cdots\;\;\underline{x}^{m}]^T
 \end{equation}
 the previous relation becomes
\begin{equation}  \label{vetfi}
 {\bm \phi}(x) = e^{H x_0}  D_{\tilde{c}}  \, \, {\bm \xi}(\underline{x}).
 \end{equation}
Now, the entries of such vector are generalized Appell polynomials with respect to $\partial$ if  the two
properties in Definition~\ref{generalizedAS}
are verified. It is evident that this could lead to impose some constraints on the diagonal coefficients of $D_{\tilde{c}}.$ In order to deduce them,
first of all we need to rewrite in matrix form the action of $\partial_{\underline{x}}$ on the vector ${\bm \xi}(\underline{x}).$ By virtue to the fact that  \cite[p. 219]{desosu92}
\begin{equation*}
\partial_{\underline{x}}(\underline{x}^{k})=\left\{ \begin{array}{ll}
         - k \, \underline{x}^{k-1}, & \mbox{$k$ even}\\
         - (n+k-1) \, \underline{x}^{k-1},  \quad & \mbox{$k$ odd,}
        \end{array} \right.
\end{equation*}
introducing the matrix $\tilde H$ defined  by
\begin{equation} \label{Ht}
(\tilde{H})_{ij}=\left\{
    \begin{array}{ll}
      -(n+i-1), & \quad  i=j+1\wedge j\; \hbox{even} \\
      -i, & \quad  i=j+1\wedge j\; \hbox{odd} \\
      0, & \quad  \hbox{otherwise,} \qquad \qquad \quad i,j=0,1,\dots,m,
    \end{array}
  \right.
\end{equation}
we obtain
 \begin{equation}\label{tildeH}
    \partial_{\underline{x}}  \, {\bm \xi}({\underline{x}}) = \tilde{H}{\bm \xi}(\underline{x}),
 \end{equation}
which means that $ \tilde{H}$ is the derivation matrix playing the role of $\partial_{\underline{x}}.$

 \begin{oss} \label{complex}
In the complex case ($n=1$)  we obtain simply $\tilde{H}=-H.$ Denoting by
$w=x_0 + \underline{x},$ with $\underline{x}= ix_1,$ it is
easily  checked that  for ${\bm \xi}({w})= [1\;\;w\;\;w^2\;\cdots\;w^{m}]^T$
one has  $\overline{\partial} {\bm \xi}({w})= \mathbf{0}^T,$ the null vector,
and $\partial {\bm \xi}({w})= H  {\bm \xi}({w}),$ i.e., the entries of
${\bm \xi}({w})$ are hypercomplex Appell polynomials.
 \end{oss}

Therefore,  by applying  to both sides   in (\ref{vetfi})  the hypercomplex differential operator
 we obtain
 \begin{eqnarray} \label{parfi}
\bar{\partial}{\bm \phi}(x)&=&\frac{1}{2}
(\partial_0+\partial_{\underline{x}})(e^{H x_{0}}D_{\tilde{c}} \ {\bm \xi}(\underline{x}))   \nonumber
\\
&=&\frac{1}{2}\left[(\partial_0 e^{H x_{0}})D_{\tilde{c}}\ {\bm \xi}(\underline{x})+
e^{H x_{0}}D_{\tilde{c}}(\partial_{\underline{x}} {\bm \xi}(\underline{x}))\right]  \nonumber  \\
&=&\frac{1}{2}\left[(H e^{H x_{0}})D_{\tilde{c}} \ {\bm \xi}(\underline{x})+e^{H x_{0}}D_{\tilde{c}}\ (\tilde{H}{\bm \xi}(\underline{x}))\right]  \nonumber \\
& =&\frac{1}{2}e^{H x_{0}}\left[HD_{\tilde{c}}+D_{\tilde{c}}\ \tilde{H}\right] {\bm \xi}(\underline{x}).
\end{eqnarray}
\begin{oss}
Notice that due to the uniqueness theorem for monogenic functions (cf. \cite{ma04} or \cite[p.180]{GHSp08}) $D_{\tilde{c}}$ must be \textit{non-singular} and we can suppose that all $\tilde{c}_k \neq 0,\; k=0,1,\dots.$ Otherwise the restriction of a component $\phi_j(x_0,\underline{x})$  of ${\bm \phi}(x)$ with $\tilde{c}_j=0$ to the hyperplane $x_0=0$ would have the value $\phi_j(0,\underline{x})=0$ and, consequently, be identically zero. But this contradicts the property of belonging to an Appell sequence.
\end{oss}
We can now prove the following result.
\begin{teo}\label{mainteo}
Suppose $\tilde{c}_0 \neq 0$ is a given real number. The polynomials in (\ref{vetfi})
are monogenic if their remaining coefficients satisfy the conditions
\begin{equation}\label{condic}
\tilde{c}_{2k}=\tilde{c}_{2k-1} = \frac{(2k-1)!!(n-2)!!}{(n+2k-2)!!} \, \tilde{c}_0, \qquad  k=1,
2,\dots, \, n>1.
\end{equation}
\end{teo}
\begin{proof}
For ${\bm \phi}(x)$ being monogenic, from (\ref{parfi}) the diagonal matrix $D_{\tilde{c}}$  should be of such a form that
\begin{equation}\label{equalprod}
HD_{\tilde{c}}+D_{\tilde{c}}\ \tilde{H} = O,
\end{equation}
the null matrix, which means that $H$ and $-\tilde{H}$ become similar matrices. By using the definitions of $H$ and $D_{\tilde{c}},$
the product $HD_{\tilde{c}}$ is given by
\begin{equation}\label{prod1}
(HD_{\tilde{c}})_{ij}=\left\{
    \begin{array}{ll}
      (j+1)\tilde{c}_j, & \quad  i=j+1 \\
      0, & \quad  \hbox{otherwise.}
    \end{array}
  \right.
\end{equation}
Similarly, from (\ref{Ht}) the entries of $D_{\tilde{c}}\ \tilde{H}$ are
\begin{equation}\label{prod2}
(D_{\tilde{c}}\tilde{H})_{ij}=\left\{
    \begin{array}{ll}
      (n+j)\ \tilde{c}_{j+1}, & \quad  i=j+1\wedge j\; \hbox{even} \\
      (j+1)\ \tilde{c}_{j+1}, & \quad  i=j+1\wedge j\; \hbox{odd} \\
      0, & \quad  \hbox{otherwise.}
    \end{array}
  \right.
\end{equation}
Thus, the identity (\ref{equalprod}) is evident for  $i\neq j+1$. For ($i=j+1 \wedge j\;  \hbox{odd}$) and
($i=j+1 \wedge j\; \hbox{even}$), one has
\[(j+1)\tilde{c}_{j+1}=(j+1)\tilde{c}_j, \qquad (n+j)\tilde{c}_{j+1}=(j+1)\tilde{c}_{j},\]
respectively, or, equivalently,
\[
 \tilde{c}_{2k}=\tilde{c}_{2k-1}, \qquad \tilde{c}_{2k-1} = \frac{2k-1}{n+2k-2}  \tilde{c}_{2k-2}.
\]
From these relations the assertion follows.
\end{proof}

\begin{oss}
 From (\ref{condic}) it is evident that the coefficients $\tilde{c}_j$ actually
depends on $n,$ i.e.,  $\tilde{c}_j \equiv \tilde{c}_j (n),$ for each $j \ge 1.$
\end{oss}

\begin{oss}\label{rem}
The choice of the coefficients (\ref{condic}) with $\tilde{c}_0=1$ gives the
generalized Appell sequence
$\{\mathscr{P}^{(n)}_k (x)\}_{k\geq 0}$ introduced in \cite{fama07} by a monogenic generating
exponential function.
\end{oss}

 Under the hypothesis of Theorem~\ref{mainteo} the polynomials $\phi_k(x)$ are monogenic. It remains to check if
property (ii) in  Definition~\ref{generalizedAS} is satisfied. Applying $(\ref{equalprod})$ we get
\begin{eqnarray*}
{\partial}{\bm \phi}(x)&=&\frac{1}{2}
\left[(\partial_0-\partial_{\underline{x}})(e^{H x_{0}}D_{\tilde{c}}\ {\bm \xi}(\underline{x}))\right] \\
&=& \frac{1}{2}\left[H e^{H x_{0}}D_{\tilde{c}}\ {\bm \xi}(\underline{x}) - e^{H x_{0}}D_{\tilde{c}}
\ \tilde{H} {\bm \xi}(\underline{x})\right]\\
&=&\frac{1}{2} e^{H x_{0}} \left[H D_{\tilde{c}} - D_{\tilde{c}} \ \tilde{H} \right] {\bm \xi}(\underline{x})
=He^{H x_{0}} D_{\tilde{c}}\ {\bm \xi}(\underline{x})= H{\bm \phi}(x),
\end{eqnarray*}
which, actually, is the corresponding matrix representation of the property (ii) in
Definition~\ref{generalizedAS}.

\begin{oss}
Following the recent article \cite{cafama15} on special properties of hypercomplex Appell polynomials like, for instance, three-term recurrence relations, further extensions including also all orthogonal Appell polynomial sequences obtained by the Gel'fand-Tsetlin procedure, cf. \cite{brasls10, lav12}, can be constructed.

In order to get such extensions, the suitable formal replacement in  (\ref{compacAp}) and (\ref{p0}) should be $x$ by $x_0$ and $c_j$ by  $\tilde{c}_j \,  \underline{x}^j\, Q_s(\underline{x}),$  where $Q_s(\underline{x})$ is an arbitrary chosen monogenic polynomial of fixed degree $s=0,1,2,\ldots.$

One should  notice that $Q_s(\underline{x})$ as an arbitrary chosen polynomial monogenic with respect to the generalized Cauchy-Riemann operator is automatically a generalized constant since it does not depend from $ x_0$ and therefore belongs also to the kernel of the conjugated generalized Cauchy-Riemann operator. This fact implies that from the point of view of the hypercomplex derivative the polynomial $Q_s(\underline{x})$ behaves like an ordinary constant number, i.e., its hypercomplex derivative is constant zero. Moreover, one could come to the conclusion that the first polynomial of hypercomplex Appell sequences could be such an initial monogenic generalized constant $Q_s(\underline{x})$  and all polynomials of higher degree would have $Q_s(\underline{x})$ as common factor. This idea is realized in the paper \cite{pe13}.

The matrix form of these multivariate homogeneous polynomials is
\[
 {\bm \phi}(x) = Q_s(\underline{x}) \, \,e^{H x_0}  D_{\tilde{c}}  \, \, {\bm \xi}(\underline{x})
 \]
and the constraints to be imposed on the diagonal entries of $ D_{\tilde{c}}$ in order to get generalized Appell polynomials are
\[
\tilde{c}_{2k}=\tilde{c}_{2k-1} = \frac{(2k-1)!!(n+2s-2)!!}{(n+2k+2s-2)!!} \, \tilde{c}_0, \qquad  k=1,
2,\dots, \, n>1.
\]

Such constraints are achieved like in Theorem~\ref{mainteo}, but taking into account that the action of $\partial_{\underline{x}}$ on the vector ${\bm \xi}(\underline{x})$ is performed by the matrix $\tilde{H^s}\, Q_s(\underline{x})$, where
\begin{equation*}
(\tilde{H^s})_{ij}=\left\{
    \begin{array}{ll}
      -(n+i+2s-1), & \quad  i=j+1\wedge j\; \hbox{even} \\
      -i, & \quad  i=j+1\wedge j\; \hbox{odd} \\
      0, & \quad  \hbox{otherwise,} \qquad \qquad \quad i,j=0,1,\dots,m.
    \end{array}
  \right.
\end{equation*}
This matrix results from the fact that
\begin{equation*}
\partial_{\underline{x}}(\underline{x}^{k}Q_s(\underline{x}))=\left\{ \begin{array}{ll}
         - k \, \underline{x}^{k-1}\, Q_s(\underline{x}), & \mbox{$k$ even}\\
         - (n+k+2s-1) \, \underline{x}^{k-1}\, Q_s(\underline{x}),  \quad & \mbox{$k$ odd.}
        \end{array} \right.
\end{equation*}
\end{oss}

\section{Special families of monogenic Appell polynomials}
We recall that the real Appell polynomials, $p_k(x),$ may also
be characterized  in terms of their generating exponential function of the
form
\[
f(t)e^{tx}=\sum_{n=0}^{+\infty} p_n(x)\frac{t^n}{n!},
\]
 where
\begin{eqnarray}\label{formserie}
f(t)=\sum_{s=0}^{+\infty} c_s\frac{t^s}{s!},\qquad c_0\neq 0.
\end{eqnarray}
The creation matrix $H$ and formula (\ref{formserie}) are the essential tools to get a  \textit{transfer matrix}
 whose action is to change the vector of monomial powers ${\bm \xi}(x)$ into the
associated Appell vector ${\bf p}(x).$  Such matrix is $f(H)$ and it has been proved in
\cite[Theorem~3.2]{ace15} that
\[{\bf p}(x)=f(H)  {\bm \xi}(x).\]

The generalized Appell sequence  $\{\mathscr{P}^{(n)}_k (x)\}_{k\geq 0} $ referred in the Remark~\ref{rem} allows to define  a generalized exponential function in the hypercomplex setting by
\begin{equation}\label{gexp}
\hbox{Exp}_n(x)\equiv e^{x_0}F(\underline{x})=\sum_{k=0}^{+\infty}\frac{\mathscr{P}^{(n)}_k (x)}{k!}, \qquad x
\in \mathcal{A}_n,
\end{equation}
where \begin{eqnarray*}
F(\underline{x} )=\sum_{s=0}^{+\infty} \tilde{c}_s \frac{\underline{x}^s}{s!},
\end{eqnarray*}
with $\tilde{c}_0=1$ and, for each $s \ge 1,$  $\tilde{c}_s$ satisfying (\ref{condic}).

Due to the homogeneity of $\mathscr{P}^{(n)}_k (x),$ we have
\[\hbox{Exp}_n(tx)=\sum_{k=0}^{+\infty}\mathscr{P}^{(n)}_k (x)\frac{t^k}{k!},\]
showing that the generating function of the Appell sequence $\{\mathscr{P}^{(n)}_k
(x)\}_{k\geq 0} $ is of the form $G(x,t)=\hbox{Exp}_n(tx),\;x\in \mathcal{A}_n, t\in \mathbb{R}$.

\begin{oss}
In the complex case, $G(x,t)=e^{tx},\; x\in \mathbb{C}, \;t\in \mathbb{R}$  is
the generating function of the basic Appell sequence
$\{x^k\}_{k\geq 0}$. Thus, $\mathscr{P}^{(n)}_k (x)$ behave as monomial
functions in the sense of the complex powers
$w^k=(x_0+ix_1)^k,\;k \ge 0.$
\end{oss}

Using the hypercomplex exponential function (\ref{gexp}) we
define generalized Appell polynomials, in general
non-homogeneous, as follows:
\begin{definition}
The sequence $\{\phi_k(x)\}_{k\geq 0}$ whose generating
function is $G(x,t)=f(t)\hbox{Exp}_n(tx),\;x\in \mathcal{A}_n$ and
$f(t)$ a formal power series as in (\ref{formserie}) is called
generalized Appell sequence.
\end{definition}
Noting that the function $f(t)$ coincides with the one
appearing in the generating function of real Appell
polynomials, the transfer matrix for ${\bm \phi}(x), x\in \mathcal{A}_n,$ is also
$f(H)$. However, in this case, $f(H)$ transforms the vector
\[
{\bm \xi}(\mathscr{P}^{(n)}(x))=[\mathscr{P}_0^{(n)}(x)\;\;\mathscr{P}_1^{(n)}
(x)\cdots \mathscr{P}_m^{(n)}
(x)]^T
\]
in the generalized Appell vector, that is
$${\bm \phi}(x)=f(H) \, {\bm \xi}(\mathscr{P}^{(n)}(x)).$$

Different kinds of generalized Appell polynomials can be
derived by appropriate choice of $f(t).$  For
instance, the corresponding transfer matrices for
generalized Bernoulli, Frobenius-Euler, and monic Hermite
polynomials are (see \cite{ace15})
\[
\left(\sum_{k=0}^m\frac{H^k}{(k+1)!}\right)^{-1}, \qquad
(1-\lambda)(P-\lambda I)^{-1}, \qquad
e^{-H^2/4}=\sum_{k=0}^m\frac{(-H^2)^k}{2^{2k}k!},
\]
respectively. In particular, setting $\lambda=-1$ in the transfer
matrix of generalized Frobenius-Euler polynomials, we get generalized Euler
polynomials.

\begin{oss}
We notice that, when $\underline{x}=0$,
\[
{\bm \xi}(\mathscr{P}^{(n)}(x_0))\equiv {\bm \xi}(x_0)
\]
while the components of ${\bm \phi} (x_0)$ are the classical real Appell
polynomials entries of the vector ${\bf p} (x_0).$
\end{oss}