## Hypergeometric Functions and Carlitz Differential Equations over Function Fields

### Anatoly N. Kochubei Institute of Mathematics, National Academy of Sciences of Ukraine, Tereshchenkivska 3, Kiev, 01601 Ukraine E-mail: kochubei@i.com.ua

Abstract
The paper is a survey of recent results in analysis of additive functions over function fields motivated by applications to various classes of special functions including Thakur's hypergeometric function. We consider basic notions and results of calculus, analytic theory of differential equations with Carlitz derivatives (including a counterpart of regular singularity), umbral calculus, holonomic modules over the Weyl-Carlitz ring. 2000 Mathematics Subject Classification: Primary 11S80, 12H25, 33E50. Secondary 05A40, 11G09, 16S32, 32C38.

1 INTRODUCTION

Let $K$  be the field of formal Laurent series $t={\sum }_{j=N}^{\infty }{\xi }_{j}{x}^{j}$  with coefficients ${\xi }_{j}$  from the Galois field $\mathbb{F}q$  , ${\xi }_{N}\ne 0$  if $t\ne 0$  , $q={p}^{\upsilon }$  , $\upsilon \in {\mathbf{Z}}_{+}$  , where $p$  is a prime number. It is well known that any non-discrete locally compact field of characteristic $p$  is isomorphic to such $K$  . The absolute value on $K$  is given by $|t|={q}^{-N}$  , $|0|=0$  . The ring of integers $O=\left\{t\in K:|t|\le 1\right\}$  is compact in the topology corresponding to the metric $dist\left(t,s\right)=|t-s|$  . The absolute value $|\cdot |$  can be extended in a unique way onto the completion $\overline{K}c$  of an algebraic closure of $K$  .
Analysis over $K$  and $\overline{K}c$  initiated by Carlitz [3and developed subsequently by Wagner, Goss, Thakur, the author, and others (see the bibliography in [12, 37) is very different from the classical calculus. The new features begin with an appropriate version of the factorial invented by Carlitz – since the usual factorial $i!$  , seen as an element of $K$  , vanishes for $i\ge p$  , Carlitz introduced the new one as
 $\begin{array}{c}{D}_{i}=\left[i\right]\left[i-1{\right]}^{q}\dots \left[1{\right]}^{{q}^{i-1}},\left[i\right]={x}^{{q}^{i}}-x\left(i\ge 1\right),{D}_{0}=1.\end{array}$ (1.1)
An important feature is the availability of many non-trivial $\mathbb{F}q$  -linear functions, that is such functions $f$  defined on $\mathbb{F}q$  -subspaces ${K}_{0}\subset K$  that $f\left({t}_{1}+{t}_{2}\right)=f\left({t}_{1}\right)+f\left({t}_{2}\right),f\left(\alpha t\right)=\alpha f\left(t\right),$  for any $t,{t}_{1},{t}_{2}\in {K}_{0}$  , $\alpha \in \mathbb{F}q$  . Such are, for example, polynomials and power series of the form $\sum {a}_{k}{t}^{{q}^{k}}$  , in particular, the Carlitz exponential
 $\begin{array}{c}{e}_{C}\left(t\right)={\sum }_{n=0}^{\infty }\frac{{t}^{{q}^{n}}}{{D}_{n}},|t|<1,\end{array}$ (1.2)
and its composition inverse, the Carlitz logarithm
 $\begin{array}{c}{log}_{C}\left(t\right)={\sum }_{n=0}^{\infty }\left(-1{\right)}^{n}\frac{{t}^{{q}^{n}}}{{L}_{n}},|t|<1,\end{array}$ (1.3)
where ${L}_{n}=\left[n\right]\left[n-1\right]\cdots \left[1\right]$  ( $n\ge 1$  ), ${L}_{0}=1$  . The notion of the Carlitz exponential obtained a wide generalization in the theory of Drinfeld modules (see [12, 37).
Among other special classes of $\mathbb{F}q$  -linear functions there are various polynomial systems (see below), an analog of the Bessel functions [4, 34, and Thakur's hypergeometric function [35, 36, 37.
The latter is defined as follows.
For $n\in {\mathbb{Z}}_{+}$  , $a\in \mathbb{Z}$  , denote
 $\begin{array}{c}\left(a{\right)}_{n}=\left\{\begin{array}{cc}{D}_{n+a-1}^{-\left(a-1\right)},& \text{if}a\ge 1\text{;}\\ {L}_{-a-n}^{-{q}^{n}},& \text{if}a\le 0,n\le -a\text{;}\\ 0,& \text{if}a\le 0,n>-a\text{}.\end{array}\end{array}$ (1.4)
Then, for ${a}_{i},{b}_{i}\in \mathbb{Z}$  , such that the series below makes sense, we set
 $\begin{array}{c}{}_{r}{F}_{s}\left({a}_{1},\dots ,{a}_{r};{b}_{1},\dots ,{b}_{s};z\right)={\sum }_{n=0}^{\infty }\frac{\left({a}_{1}{\right)}_{n}\cdots \left({a}_{r}{\right)}_{n}}{\left({b}_{1}{\right)}_{n}\cdots \left({b}_{s}{\right)}_{n}{D}_{n}}{z}^{{q}^{n}}.\end{array}$ (1.5)
Thakur [35, 36, 37has carried out a thorough investigation of the functions (1.5) and obtained analogs of many properties known for the classical situation. In particular, he found an analog of the hypergeometric differential equation. Its main ingredients are the difference operator $\left(\Delta u\right)\left(t\right)=u\left(xt\right)-xu\left(t\right)$  introduced by Carlitz [3, the nonlinear ( $\mathbb{F}q$  -linear) operator $d=\sqrt[q]{}\circ \Delta$  , and the $\mathbb{F}q$  -linear Frobenius operator $\tau u={u}^{q}$  . For example, the function $y{=}_{2}{F}_{1}\left(a,b;c;z\right)$  is a solution of the equation
 $\begin{array}{c}\left(\Delta -\left[-a\right]\right)\left(\Delta -\left[-b\right]\right)y=d\left(\Delta -\left[1-c\right]\right)y.\end{array}$ (1.6)
Here we touch only a part of Thakur's results (he considered also hypergeometric functions corresponding to other places of $\mathbb{F}q\left(x\right)$  , a version of (1.5) with parameters from $K$  and its extensions etc).
The Carlitz exponential ${e}_{C}$  satisfies a much simpler equation of the same kind:
 $\begin{array}{c}d{e}_{C}={e}_{C},\end{array}$ (1.7)
so that the operator $d$  may be seen as an analog of the derivative. The operator $\tau$  is an analog of the multiplication by $t$  , so that $\Delta$  is the counterpart of $t\frac{d}{dt}$  .
The same operators appear in the positive characteristic analogs of the canonical commutation relations of quantum mechanics [15, 16. In the analog of the Schrödinger representation we consider, on the Banach space ${C}_{0}\left(O,\overline{K}c\right)$  of continuous $\mathbb{F}q$  -linear functions on $O$  , with values from $\overline{K}c$  (with the supremum norm), the “creation and annihilation operators” ${a}^{+}=\tau -I,{a}^{-}=d$  ( $I$  is the identity operator). Then
 $\begin{array}{c}{a}^{-}{a}^{+}-{a}^{+}{a}^{-}=\left[1{\right]}^{1/q}I,\end{array}$ (1.8)
the operator ${a}^{+}{a}^{-}$  possesses the orthonormal (in the non-Archimedean sense [32) eigenbasis $\left\{{f}_{i}\right\}$  ,
 $\begin{array}{c}\left({a}^{+}{a}^{-}\right){f}_{i}=\left[i\right]{f}_{i},i=0,1,2,\dots ;\end{array}$ (1.9)
${a}^{+}$  and ${a}^{-}$  act upon the basis as follows:
 $\begin{array}{c}{a}^{+}{f}_{i-1}=\left[i\right]{f}_{i},{a}^{-}{f}_{i}={f}_{i-1},i\ge 1;{a}^{-}{f}_{0}=0.\end{array}$ (1.10)
Here $\left\{{f}_{i}\right\}$  is the sequence of normalized Carlitz polynomials
 $\begin{array}{c}{f}_{i}\left(s\right)={D}_{i}^{-1}{\prod }_{m\in \mathbb{F}q\left[x\right]degm (1.11)
which forms an orthonormal basis in ${C}_{0}\left(O,\overline{K}c\right)$  . The spectrum of the “number operator” ${a}^{+}{a}^{-}$  is the set of elements $\left[i\right]$  , so that even this notation (proposed by Carlitz in 1935) becomes parallel to the usual quantum mechanical situation.
An analog of the Bargmann-Fock representation is obtained if we consider the operators of almost the same form, ${\stackrel{~}{a}}^{+}=\tau ,{\stackrel{~}{a}}^{-}=d,$  but on the Banach space $H$  of power series $u\left(t\right)={\sum }_{n=0}^{\infty }{a}_{n}\frac{{t}^{{q}^{n}}}{{D}_{n}}$  with ${a}_{n}\in \overline{K}c$  , ${a}_{n}\to 0$  as $n\to \infty$  . These new operators satisfy the same relations (1.8)-(1.10), but this time instead of the Carlitz polynomials ${f}_{n}$  we get the eigenfunctions ${\stackrel{~}{f}}_{n}=\frac{{t}^{{q}^{n}}}{{D}_{n}}$  .
The above results motivated the author to begin to develop analysis and theory of differential equations for $\mathbb{F}q$  -linear functions over $K$  and $\overline{K}c$  , that is for the case which can be seen as a concentrated expression of features specific for the analysis in positive characteristic. This paper is a brief survey of some achievements in this direction. In particular, we consider the counterparts of the basic notions of calculus, analytic theory of differential equations (in the regular case and the case of regular singularity), their applications to some special functions, like the power function, logarithm and polylogarithms, Thakur's hypergeometric function etc.
An umbral calculus and a theory of holonomic modules are initiated for this case. Like in the classical situation (see [5), it is shown that some basic objects of the function field arithmetic generate holonomic modules.

2 Calculus

2.1. Higher Carlitz operators ${\Delta }^{\left(n\right)}$  are introduced recursively,
 $\begin{array}{c}\left({\Delta }^{\left(n\right)}u\right)\left(t\right)={\Delta }^{\left(n-1\right)}u\left(xt\right)-{x}^{{q}^{n-1}}{\Delta }^{\left(n-1\right)}u\left(t\right),n\ge 2.\end{array}$ (2.1)
For $n=1$  , the formula (2.1) coincides with the definition of $\Delta ={\Delta }^{\left(1\right)}$  , if we set ${\Delta }^{\left(0\right)}=I$  .
The first application of these operators is the reconstruction formula [16for the coefficients ${a}_{n}$  of a power series $u\in H$  . Note that the classical formula does not make sense here because it contains the expression ${u}^{\left(n\right)}\left(t\right)/n!$  where both the numerator and denominator vanish.
Theorem 2.1. If $u\in H$  , then ${a}_{n}={lim}_{t\to 0}\frac{{\Delta }^{\left(n\right)}u\left(t\right)}{{t}^{{q}^{n}}},n=0,1,2,\dots .$
For a continuous non-holomorphic $\mathbb{F}q$  -linear function $u$  the behaviour of the functions ${\mathfrak{D}}^{k}u\left(t\right)={t}^{-{q}^{k}}{\Delta }^{\left(k\right)}u\left(t\right),t\in O\\left\{0\right\},$  near the origin measures the smoothness of $u$  . We say that $u\in {C}_{0}^{k+1}\left(O,\overline{K}c\right)$  if ${\mathfrak{D}}^{k}u$  can be extended to a continuous function on $O$  . This includes the case ( $k=0$  ) of differentiable functions.
The next theorem proved in [16gives a characterization of the above smoothness in terms of coefficients of the Fourier-Carlitz expansion. It includes, as a particular case ( $k=0$  ), the characterization of differentiable $\mathbb{F}q$  -linear functions obtained by Wagner [41.
Theorem 2.2. A function $u={\sum }_{n=0}^{\infty }{c}_{n}{f}_{n}\in {C}_{0}\left(O,\overline{K}c\right)$  belongs to ${C}_{0}^{k+1}\left(O,\overline{K}c\right)$  if and only if ${q}^{n{q}^{k}}|{c}_{n}|\to 0\text{for}n\to \infty .$  In this case ${sup}_{t\in O}|{\mathfrak{D}}^{k}u\left(t\right)|={sup}_{n\ge k}{q}^{\left(n-k\right){q}^{k}}|{c}_{n}|.$
For a generalization to some classes of not necessarily $\mathbb{F}q$  -linear functions see [43.
Similarly [16, a function $u$  is analytic (that is, $u\left(t\right)=\sum {a}_{i}{t}^{{q}^{i}}$  ) on the ball $O$  if and only if ${q}^{\frac{{q}^{n}}{q-1}}|{c}_{n}|\to 0$  , as $n\to \infty$  . A more refined result by Yang [42, useful in many applications, which makes it possible to find an exact domain of analyticity, is as follows (again we consider only $\mathbb{F}q$  -linear functions while in [42a more general class is studied).
Theorem 2.3 (Yang). A function $u={\sum }_{n=0}^{\infty }{c}_{n}{f}_{n}\in {C}_{0}\left(O,\overline{K}c\right)$  is locally analytic if and only if
 $\begin{array}{c}\gamma ={liminf}_{n\to \infty }\left\{-{q}^{-n}{log}_{q}|{c}_{n}|\right\}>0,\end{array}$ (2.2)
and if (2.2) holds, then $u$  is analytic on any ball of the radius ${q}^{-l}$  , $l=max\left(0,\left[-\left(log\left(q-1\right)+log\gamma \right)/logq\right]+1\right).$
2.2. Viewing $d$  as a kind of a derivative, it is natural to introduce an antiderivative $S$  setting $Sf=u$  where $u$  is a solution of the equation $du=f$  , with the normalization $u\left(1\right)=0$  .
It is easy to find $Sf$  explicitly if $f$  is given by its Fourier-Carlitz expansion (see [16).
Next, we introduce a Volkenborn-type integral of a function $f\in {C}_{0}^{1}\left(O,\overline{K}c\right)$  (see [32for a similar integration theory over ${\mathbb{Z}}_{p}$  ) setting ${\int }_{O}f\left(t\right)dt\stackrel{\text{def}}{=}{lim}_{n\to \infty }\frac{Sf\left({x}^{n}\right)}{{x}^{n}}=\left(Sf{\right)}^{\prime }\left(0\right).$  The integral is a $\mathbb{F}q$  -linear continuous functional on ${C}_{0}^{1}\left(O,\overline{K}c\right)$  , ${\int }_{O}cf\left(t\right)dt={c}^{q}{\int }_{O}f\left(t\right)dt,c\in \overline{K}c,$  possessing the following “invariance” property (related, in contrast to the case of ${\mathbf{Z}}_{p}$  , to the multiplicative structure):
${\int }_{O}f\left(xt\right)dt=x{\int }_{O}f\left(t\right)dt-{f}^{q}\left(1\right).$  Calculating the integrals of some important functions we obtain new relations between them.
In addition to the Carlitz exponential ${e}_{C}$  and the Carlitz polynomials ${f}_{n}$  (see (1.2) and (1.11)), we mention the Carlitz module function
 $\begin{array}{c}{C}_{s}\left(z\right)={\sum }_{i=0}^{\infty }{f}_{i}\left(s\right){z}^{{q}^{i}},s\in O,|z|<1.\end{array}$ (2.3)
Note that if $s\in \mathbb{F}q\left[x\right]$  , then only the terms with $i\le degs$  are different from zero in (2.3).
We have ${\int }_{O}{t}^{{q}^{n}}dt=-\frac{1}{\left[n+1\right]},n=0,1,2,\dots ;$  ${\int }_{O}{f}_{n}\left(t\right)dt=\frac{\left(-1{\right)}^{n+1}}{{L}_{n+1}},n=0,1,2,\dots ;$  ${\int }_{O}{C}_{s}\left(z\right)ds={log}_{C}\left(z\right)-z,z\in K,|z|<1;$  ${\int }_{O}{e}_{C}\left(st\right)ds=t-{e}_{C}\left(t\right),t\in K,|t|<1.$  For the proofs see [16.

3 Differential Equations for $\mathbb{F}q$  -Linear Functions

3.1. Let us consider function field analogs of linear differential equations with holomorphic or polynomial coefficients. Note that in our situation the meaning of a polynomial coefficient is not a usual multiplication by a polynomial, but the action of a polynomial in the operator $\tau$  .
We begin with the regular case and consider an equation (actually, a system)
 $\begin{array}{c}dy\left(t\right)=P\left(\tau \right)y\left(t\right)+f\left(t\right)\end{array}$ (3.1)
where for each $z\in {\left(\overline{K}c\right)}^{m}$  , $t\in K$  ,
 $\begin{array}{c}P\left(\tau \right)z={\sum }_{k=0}^{\infty }{\pi }_{k}{z}^{{q}^{k}},f\left(t\right)={\sum }_{j=0}^{\infty }{\phi }_{j}\frac{{t}^{{q}^{j}}}{{D}_{j}},\end{array}$ (3.2)
${\pi }_{k}$  are $m×m$  matrices with elements from $\overline{K}c$  , ${\phi }_{j}\in {\left(\overline{K}c\right)}^{m}$  , and it is assumed that the series (3.2) have positive radii of convergence. The action of the operator $\tau$  upon a vector or a matrix is defined component-wise, so that ${z}^{{q}^{k}}=\left({z}_{1}^{{q}^{k}},\dots ,{z}_{m}^{{q}^{k}}\right)$  for $z=\left({z}_{1},\dots ,{z}_{m}\right)$  .
We seek a $\mathbb{F}q$  -linear solution of (3.1) on some neighbourhood of the origin, of the form
 $\begin{array}{c}y\left(t\right)={\sum }_{i=0}^{\infty }{y}_{i}\frac{{t}^{{q}^{i}}}{{D}_{i}},{y}_{i}\in {\left(\overline{K}c\right)}^{m},\end{array}$ (3.3)
where ${y}_{0}$  is a given element, so that the “initial” condition for our situation is
 $\begin{array}{c}{lim}_{t\to 0}{t}^{-1}y\left(t\right)={y}_{0}.\end{array}$ (3.4)
The next theorem, proved in [17, is the function field analog of the Cauchy theorem from the classical analytic theory of differential equations.
Theorem 3.1. For any ${y}_{0}\in {\left(\overline{K}c\right)}^{m}$  the equation (3.1) has a unique local solution of the form (3.3), which satisfies (3.4), with the series having a positive radius of convergence.
Thus, regular equations with Carlitz derivatives behave more or less as their classical counterparts. The situation is different for singular equations. Let us consider scalar equations of arbitrary order
 $\begin{array}{c}{\sum }_{j=0}^{m}{A}_{j}\left(\tau \right){d}^{j}u=f\end{array}$ (3.5)
where $f\left(t\right)={\sum }_{n=0}^{\infty }{\phi }_{n}\frac{{t}^{{q}^{n}}}{{D}_{n}}$  , ${A}_{j}\left(\tau \right)$  are power series having (as well as the one for $f$  ) positive radii of convergence.
We investigate formal solutions of (3.5), of the form
 $\begin{array}{c}u\left(t\right)={\sum }_{n=0}^{\infty }{u}_{n}\frac{{t}^{{q}^{n}}}{{D}_{n}},{u}_{n}\in \overline{K}c.\end{array}$ (3.6)
One can apply an operator series $A\left(\tau \right)={\sum }_{k=0}^{\infty }{\alpha }_{k}{\tau }^{k}$  (even without assuming its convergence) to a formal series (3.6), setting ${\tau }^{k}u\left(t\right)={\sum }_{n=0}^{\infty }{u}_{n}^{{q}^{k}}\left[n+1{\right]}^{{q}^{k-1}}\dots \left[n+k\right]\frac{{t}^{{q}^{n+k}}}{{D}_{n+k}},k\ge 1,$  and $A\left(\tau \right)u\left(t\right)={\sum }_{l=0}^{\infty }\frac{{t}^{{q}^{l}}}{{D}_{l}}{\sum }_{n+k=l}{\alpha }_{k}{u}_{n}^{{q}^{k}}\left[n+1{\right]}^{{q}^{k-1}}\dots \left[n+k\right]$  where the factor $\left[n+1{\right]}^{{q}^{k-1}}\dots \left[n+k\right]$  is omitted for $k=0$  . These formal manipulations are based on the identity $\tau \left(\frac{{t}^{{q}^{i-1}}}{{D}_{i-1}}\right)=\left[i\right]\frac{{t}^{{q}^{i}}}{{D}_{i}}.$  Using also the relation $d\left(\frac{{t}^{{q}^{i}}}{{D}_{i}}\right)=\frac{{t}^{{q}^{i-1}}}{{D}_{i-1}},$  now we can give a meaning to the notion of a formal solution of the equation (3.5).
Theorem 3.2. Let $u\left(t\right)$  be a formal solution (3.6) of the equation (3.5), where the series for ${A}_{j}\left(\tau \right)z$  , $z\in \overline{K}c$  , and $f\left(t\right)$  , have positive radii of convergence. Then the series (3.6) has a positive radius of convergence.
This result (proved in [17) is in a strong contrast to the classical theory. Note that in the $p$  -adic case a similar phenomenon takes place for equations satisfying certain strong conditions upon zeros of indicial polynomials [1, 6, 25, 33. In our case such a behavior is proved for any equation, which resembles the (much simpler) case [25of differential equations over a field of characteristics zero, whose residue field also has characteristic zero.
3.2. The equations (3.1) and (3.5) behave like linear equations, though they are actually only $\mathbb{F}q$  -linear. Theorem 3.1 can be extended [20to the case of strongly nonlinear equations (containing self-compositions $y\circ y\circ \cdots \circ y$  ).
On the other hand, it is natural to consider some equations of this kind in wider classes of $\mathbb{F}q$  -linear functions resembling meromorphic functions of a complex variables. The set ${\mathcal{ℛ}}_{K}$  of locally convergent $\mathbb{F}q$  -linear holomorphic functions forms a non-commutative ring with respect to the composition operation (the pointwise multiplication violates the $\mathbb{F}q$  -linearity). The non-commutativity of ${\mathcal{ℛ}}_{K}$  makes the algebraic structures related to Carlitz differential equations much more complicated compared to their classical counterparts. So far their understanding is only at its initial stage. It is known, however, that ${\mathcal{ℛ}}_{K}$  can be imbedded into a skew field of $\mathbb{F}q$  -linear “meromorphic” series containing terms like ${t}^{{q}^{-k}}$  (see [20). A deep investigation of bi-infinite series of this kind convergent on the whole of $\overline{K}c$  has been carried out by Poonen [24.
A specific class of equations with solutions meromorphic in the above sense is the class of scalar Riccati-type equations
 $\begin{array}{c}dy\left(t\right)=\lambda \left(y\circ y\right)\left(t\right)+\left(P\left(\tau \right)y\right)\left(t\right)+R\left(t\right)\end{array}$ (3.7)
where $\lambda \in \overline{K}c$  , $\left(P\left(\tau \right)y\right)\left(t\right)={\sum }_{k=1}^{\infty }{p}_{k}{y}^{{q}^{k}}\left(t\right),R\left(t\right)={\sum }_{k=0}^{\infty }{r}_{k}{t}^{{q}^{k}},$  ${p}_{k},{r}_{k}\in \overline{K}c$  (note that the right-hand side of (3.7) does not contain the linear term). The following theorem is proved in [20.
Theorem 3.3. If $0<|\lambda |\le {q}^{-1/{q}^{2}}$  , $|{p}_{k}|\le {q}^{-1/{q}^{2}}$  , $|{r}_{k}|\le {q}^{-1/{q}^{2}}$  for all $k$  , then the equation (3.7) possesses solutions of the form $y\left(t\right)=c{t}^{1/q}+{\sum }_{n=0}^{\infty }{a}_{n}{t}^{{q}^{n}},c,{a}_{n}\in \overline{K}c,c\ne 0,$  where the series converges on the open unit disk $|t|<1$  .

4 Regular Singularity

4.1. In analysis over $\mathbb{C}$  , a typical class of systems with regular singularity at the origin $\zeta =0$  over $\mathbb{C}$  consists of systems of the form
 $\begin{array}{c}\zeta {y}^{\prime }\left(\zeta \right)=\left(B+{\sum }_{k=1}^{\infty }{A}_{k}{\zeta }^{k}\right)y\left(\zeta \right)\end{array}$ (4.1)
where $B,{A}_{j}$  are constant matrices, and the series converges on a neighbourhood of the origin.
Such a system possesses a fundamental matrix solution of the form $W\left(\zeta \right){\zeta }^{C}$  where $W\left(\zeta \right)$  is holomorphic on a neighbourhood of zero, $C$  is a constant matrix, ${\zeta }^{C}=exp\left(Clog\zeta \right)$  is defined by the obvious power series. Under some additional assumptions regarding the eigenvalues of the matrix $B$  , one can take $C=B$  . For similar results over ${\mathbb{C}}_{p}$  see [9.
In order to investigate such a class of equations in the framework of $\mathbb{F}q$  -linear analysis over $K$  , one has to go beyond the class of locally analytic functions. Instead of power series expansions we can use the expansions in Carlitz polynomials on the compact ring $O\subset K$  . The property of local analyticity, if it takes place, can be recovered with the use of Theorem 2.3. Note that our approach would fail if we consider equations over $\overline{K}c$  instead of $K$  (our solutions may take their values from $\overline{K}c$  , but they are defined over subsets of $K$  ). In this sense our techniques are different from the ones developed for both the characteristic zero cases.
We begin with the simplest model scalar equation
 $\begin{array}{c}\tau du=\lambda u,\lambda \in \overline{K}c,\end{array}$ (4.2)
whose solution may be seen as a function field counterpart of the power function $t↦{t}^{\lambda }$  .
We look for a continuous $\mathbb{F}q$  -linear solution $u\left(t,\lambda \right)$  of the equation (4.2), with the “initial condition” $u\left(1,\lambda \right)=1$  , in the form
 $\begin{array}{c}u\left(t\right)={\sum }_{i=0}^{\infty }{c}_{i}{f}_{i}\left(t\right),t\in O,\end{array}$ (4.3)
where ${c}_{0}=1$  .
It is easy to see that the equation (4.2) has no continuous solutions if $|\lambda |\ge 1$  . If $|\lambda |<1$  , then the solution $u\left(t,\lambda \right)$  is unique, continuous on $O$  , and the coefficients from (4.3) have the form ${c}_{n}={\prod }_{j=0}^{n-1}\left(\lambda -\left[j\right]\right).$  The function $u\left(t,\lambda \right)$  is analytic on $O$  if and only if $\lambda =\left[j\right]$  for some $j\ge 0$  ; in this case $u\left(t,\lambda \right)=u\left(t,\left[j\right]\right)={t}^{{q}^{j}}$  . If $\lambda \ne \left[j\right]$  for any integer $j\ge 0$  , then $u\left(t,\lambda \right)$  is locally analytic on $O$  if and only if $\lambda =-x$  , and in that case $u\left(t,-x\right)=0$  for $|t|\le {q}^{-1}$  . The relation $u\left({t}^{{q}^{m}},\lambda \right)=u\left(t,{\lambda }^{{q}^{m}}+\left[m\right]\right),t\in O,$  holds for all $\lambda$  , $|\lambda |<1$  , and for all $m=0,1,2,\dots$  . For the proofs see [18.
Similarly, if in (4.1) $\lambda =\left({\lambda }_{ij}\right)$  is is a $m×m$  matrix with elements from $\overline{K}c$  , and we look for a matrix-valued solution of (4.1), then such a solution is given by the series (4.3) with the matrix coefficients ${c}_{i}=\left\{{\prod }_{j=0}^{i-1}\left(\lambda -\left[j\right]{I}_{m}\right)\right\}{c}_{0},i\ge 1$  ( ${I}_{m}$  is the unit matrix), if $|\lambda |\stackrel{\text{def}}{=}max|{\lambda }_{ij}|<1$  .
4.2. The analog, for our situation, of the system (4.1) is the system
 $\begin{array}{c}\tau du-P\left(\tau \right)u=0\end{array}$ (4.4)
where $P\left(\tau \right)$  is a matrix-valued analytic function, so that $P\left(\tau \right)z={\sum }_{k=0}^{\infty }{\pi }_{k}{z}^{{q}^{k}}$  . We assume that $|{\pi }_{k}|\le \gamma$  , $\gamma >0$  , for all $k$  , $|{\pi }_{0}|<1$  . Denote by $g\left(t\right)$  a solution of the equation $\tau dg={\pi }_{0}g$  .
Let ${\lambda }_{1},\dots ,{\lambda }_{m}\in \overline{K}c$  be the eigenvalues of the matrix ${\pi }_{0}$  .
Theorem 4.1. If
 $\begin{array}{c}{\lambda }_{i}-{\lambda }_{j}^{{q}^{k}}\ne \left[k\right],i,j=1,\dots ,m;k=1,2,\dots ,\end{array}$ (4.5)
then the system (4.4) has a matrix solution $u\left(t\right)=W\left(g\left(t\right)\right),W\left(s\right)={\sum }_{k=0}^{\infty }{w}_{k}{s}^{{q}^{k}},{w}_{0}={I}_{m},$  where the series for $W$  has a positive radius of convergence.
The paper [18contains, apart from the proof of Theorem 4.1, a discussion of some situations (the Euler type equations) where its conditions are violated, as well as of the meaning of the conditions (4.5). Here we only mention that in the scalar case $m=1$  the condition (4.5) is equivalent to the assumption ${\pi }_{0}\ne -x$  , so that it excludes the case where solutions of the equation $\tau dg={\pi }_{0}g$  has pathological properties.
4.3. For the above equation, continuous solutions were found as Fourier-Carlitz expansions
 $\begin{array}{c}u\left(t\right)={\sum }_{n=0}^{\infty }{c}_{n}{f}_{n}\left(t\right),\end{array}$ (4.6)
and we had to impose certain conditions upon coefficients of the equation, in order to guarantee the uniform convergence of the series on $O$  (which is equivalent to the fact that ${c}_{n}\to 0$  ).
However formally we could write the series (4.6) for the solutions without those conditions.
Thus, it is natural to ask whether the corresponding series (4.6) converge at some points $t\in O$  .
Note that (4.6) always makes sense for $t\in \mathbb{F}q\left[x\right]$  (for each such $t$  only a finite number of terms is different from zero). The question is whether the series converges on a wider set; if the answer is negative, such a formal solution is called strongly singular.
The available results regarding strong singularity of solutions of some equations are based on the following general fact [18.
Theorem 4.2. If $|{c}_{i}|\ge \rho >0$  for all $i\ge {i}_{0}$  (where ${i}_{0}$  is some natural number), then the function (4.6) is strongly singular.
It follows from Theorem 4.2 that non-trivial formal solutions of the equation (4.2) with $|\lambda |\ge 1$  are strongly singular. A more complicated example is provided by the equation
 $\begin{array}{c}\left(\Delta -\left[-a\right]\right)\left(\Delta -\left[-b\right]\right)u=d\Delta u,a,b\in \mathbb{Z},\end{array}$ (4.7)
for Thakur's hypergeometric function ${}_{2}{F}_{1}\left(a,b;1;t\right)$  .
A holomorphic solution of (4.7) is given by an appropriate specialization of (1.5). Classically (over $\mathbb{C}$  ), there exists the second solution with a logarithmic singularity. Here the situation is different. Looking for a solution of the form (4.6) we obtain a recursive relation
 $\begin{array}{c}\left({c}_{i+2}^{1/q}-{c}_{i+2}\right)+{c}_{i+1}^{1/q}\left[i+1{\right]}^{1/q}-{c}_{i+1}\left(\left[i\right]+\left[i+1\right]-\left[-a\right]-\left[-b\right]\right)-{c}_{i}\left(\left[i\right]-\left[-a\right]\right)\left(\left[i\right]-\left[-b\right]\right)=0,i=0,1,2,\dots .\end{array}$ (4.8)
Taking arbitrary initial coefficients ${c}_{0},{c}_{1}\in \overline{K}c$  we obtain a solution $u$  defined on $\mathbb{F}q\left[x\right]$  . On each step we have to solve the equation
 $\begin{array}{c}{z}^{1/q}-z=v.\end{array}$ (4.9)
If $|{c}_{i}|\le 1$  and $|{c}_{i+1}|\le 1$  , then in the equation for ${c}_{i+2}$  we have $|v|<1$  .
It can be shown [18that the equation (4.9) has a unique solution ${z}_{0}\in \overline{K}c$  , for which $|{z}_{0}|\le |v|$  , and $q-1$  other solutions $z$  , $|z|=1$  . It is natural to call a solution generic if, starting from a certain step of finding the coefficients ${c}_{n}$  , we always take the most frequent option corresponding to a solution of (4.9) with $|z|=1$  . Now Theorem 4.2 implies the following fact.
Theorem 4.3. A generic solution of the equation (4.7) is strongly singular.
Of course, in some special cases the recursion (4.8) can lead to more regular solutions, in particular, to the holomorphic solutions found by Thakur.

5 Polylogarithms and a Zeta Function

5.1. The Carlitz differential equations can be used for defining new special functions with interesting properties. Some examples are given in this section.
An analog of the function $-log\left(1-t\right)$  is defined via the equation
 $\begin{array}{c}\left(1-\tau \right)du\left(t\right)=t,t\in {K}_{\pi },\end{array}$ (5.1)
a counterpart of the classical equation $\left(1-t\right){u}^{\prime }\left(t\right)=1$  . The next results are taken from [21where the equation (5.1) is considered for an arbitrary finite place of $\mathbb{F}q\left(x\right)$  .
Let ${l}_{1}\left(t\right)$  be a $\mathbb{F}q$  -linear holomorphic solution of (5.1) with the zero initial condition (in the sense of (3.4)). Then it is easy to show that
 $\begin{array}{c}{l}_{1}\left(t\right)={\sum }_{n=1}^{\infty }\frac{{t}^{{q}^{n}}}{\left[n\right]},\end{array}$ (5.2)
and the series in (5.2) converges for $|t|\le {q}^{-1}$  .
Note that ${l}_{1}\left(t\right)$  is different from the well-known Carlitz logarithm ${log}_{C}$  (see (1.3)). Analogies motivating the introduction of special functions are not so unambiguous, and, for instance, from the composition ring viewpoint, ${log}_{C}$  is an analog of ${e}^{-t}$  , though in other respects it is a valuable analog of the logarithm. By the way, another possible analog of the logarithm is a continuous function $u\left(t\right)$  , $|t{|}_{\pi }\le 1$  , satisfying the equation $\Delta u\left(t\right)=t$  (an analog of $t{u}^{\prime }\left(t\right)=1$  ) and the condition $u\left(1\right)=0$  . In fact, $u={\mathcal{D}}_{1}$  , the first hyperdifferential operator; see [14.
Now we consider continuous non-holomorphic extensions of ${l}_{1}$  .
Theorem 5.1. The equation (5.1) has exactly $q$  continuous solutions on $O$  coinciding with (5.2) as $|t|\le {q}^{-1}$  . These solutions have the expansions in the Carlitz polynomials $u={\sum }_{i=0}^{\infty }{c}_{i}{f}_{i}$  where ${c}_{1}$  is an arbitrary solution of the equation ${c}_{1}^{q}-{c}_{1}+1=0$  , higher coefficients are found from the relation ${c}_{n}={\sum }_{j=0}^{\infty }{\left({c}_{n-1}\left[n-1\right]\right)}^{{q}^{j+1}},n\ge 2,$  and the coefficient ${c}_{0}$  is determined by the relation ${c}_{0}={\sum }_{i=1}^{\infty }\left(-1{\right)}^{i+1}\frac{{c}_{i}}{{L}_{i}},$
Below we denote by ${l}_{1}$  an arbitrary fixed “branch” of extensions of (5.2).
The polylogarithms ${l}_{n}\left(t\right)$  are defined recursively by the equations
 $\begin{array}{c}\Delta {l}_{n}={l}_{n-1},n\ge 2,\end{array}$ (5.3)
which agree with the classical ones $t{l}_{n}^{\prime }\left(t\right)={l}_{n-1}\left(t\right)$  . Analytic $\mathbb{F}q$  -linear solutions of (5.3), such that ${t}^{-1}{l}_{n}\left(t\right)\to 0$  as $t\to 0$  , are found easily by induction:
 $\begin{array}{c}{l}_{n}\left(t\right)={\sum }_{j=1}^{\infty }\frac{{t}^{{q}^{j}}}{\left[j{\right]}^{n}},|t|\le {q}^{-1}.\end{array}$ (5.4)
Theorem 5.2. For each $n\ge 2$  , there exists a unique continuous $\mathbb{F}q$  -linear solution of the equation (5.3) coinciding for $|t|\le {q}^{-1}$  with the polylogarithm (5.4). The solution is given by the Carlitz expansion ${l}_{n}={\sum }_{i=0}^{\infty }{c}_{i}^{\left(n\right)}{f}_{i}$  with $|{c}_{i}^{\left(n\right)}|\le {C}_{n}{q}^{-{q}^{i-1}},{C}_{n}>0,i\ge 1,$  ${c}_{0}^{\left(n\right)}={\sum }_{i=1}^{\infty }\left(-1{\right)}^{i+1}\frac{{c}_{i}^{\left(n\right)}}{{L}_{i}},$
5.2. Now that the above polylogarithms have been extended onto the disk $\left\{|t|\le 1\right\}$  , we can interpret their values at $t=1$  as “special values” of a kind of a zeta function. In order to define the latter, we introduce the operator ${\Delta }^{\left(\alpha \right)}$  , $\alpha \in O$  , a function field analog of the Hadamard fractional derivative ${\left(t\frac{d}{dt}\right)}^{\alpha }$  from real analysis (see [31).
Denote by ${\mathcal{D}}_{k}\left(t\right)$  , $k\ge 0$  , $t\in O$  , the sequence of hyperdifferentiations defined initially on monomials by the relations ${\mathcal{D}}_{0}\left({x}^{n}\right)={x}^{n}$  , ${\mathcal{D}}_{k}\left(1\right)=0$  , $k\ge 1$  , ${\mathcal{D}}_{k}\left({x}^{n}\right)=\left(\genfrac{}{}{0}{}{n}{k}\right){x}^{n-k},$  where it is assumed that $\left(\genfrac{}{}{0}{}{n}{k}\right)=0$  for $k>n$  . ${\mathcal{D}}_{k}$  is extended onto $\mathbb{F}q\left[x\right]$  by $\mathbb{F}q$  -linearity, and then onto $O$  by continuity [40. The sequence $\left\{{\mathcal{D}}_{k}\right\}$  is an orthonormal basis of the space of continuous $\mathbb{F}q$  -linear functions on $O$  [14, 7.
Let $\alpha \in O$  , $\alpha ={\sum }_{n=0}^{\infty }{\alpha }_{n}{x}^{n}$  , ${\alpha }_{n}\in \mathbb{F}q$  . Denote $\stackrel{^}{\alpha }={\sum }_{n=0}^{\infty }\left(-1{\right)}^{n}{\alpha }_{n}{x}^{n}$  . For an arbitrary continuous $\mathbb{F}q$  -linear function $u$  on $O$  we define its “fractional derivative” ${\Delta }^{\left(\alpha \right)}u$  at a point $t\in O$  by the formula $\left({\Delta }^{\left(\alpha \right)}u\right)\left(t\right)={\sum }_{k=0}^{\infty }\left(-1{\right)}^{k}{\mathcal{D}}_{k}\left(\stackrel{^}{\alpha }\right)u\left({x}^{k}t\right).$  The function $\alpha ↦\left({\Delta }^{\left(\alpha \right)}u\right)\left(t\right)$  is continuous and $\mathbb{F}q$  -linear.
Our understanding of ${\Delta }^{\left(\alpha \right)}$  as a kind of a fractional derivative is justified by the following properties:
${\Delta }^{\left({x}^{n}\right)}={\Delta }^{n},n=1,2,\dots ;$  ${\Delta }^{\left(\alpha \right)}\left({\Delta }^{\left(\beta \right)}u\right)\left(t\right)=\left({\Delta }^{\left(\alpha \beta \right)}u\right)\left(t\right),$  for any $\alpha ,\beta \in O$  .
5.3. We define $\zeta \left(t\right)$  , $t\in K$  , setting $\zeta \left(0\right)=0$  , $\zeta \left({x}^{-n}\right)={l}_{n}\left(1\right),n=1,2,\dots ,$  and $\zeta \left(t\right)=\left({\Delta }^{\left({\theta }_{0}+{\theta }_{1}x+\cdots \right)}{l}_{n}\right)\left(1\right),n=1,2,\dots ,$  if $t={x}^{-n}\left({\theta }_{0}+{\theta }_{1}x+\cdots \right)$  , ${\theta }_{j}\in \mathbb{F}q$  . $\zeta$  is a continuous $\mathbb{F}q$  -linear function on ${K}_{x}$  with values in $\overline{K}c$  .
In particular, we have $\zeta \left({x}^{m}\right)=\left({\Delta }^{m+1}{l}_{1}\right)\left(1\right),m=0,1,2,\dots .$  The above definition is of course inspired by the classical polylogarithm relation $\left(z\frac{d}{dz}\right){\sum }_{n=1}^{\infty }\frac{{z}^{n}}{{n}^{s}}={\sum }_{n=1}^{\infty }\frac{{z}^{n}}{{n}^{s-1}}.$  In contrast to Goss's zeta function defined on natural numbers and interpolated onto ${\mathbb{Z}}_{p}$  (see [12, 37), the above $\zeta$  is purely an object of the characteristic $p$  arithmetic.
Let us write some relations for special values of our $\zeta$  ; for the details see [21.
As we saw, ${\sum }_{j=1}^{\infty }\frac{{t}^{{q}^{j}}}{\left[j{\right]}^{n}}={\sum }_{i=0}^{\infty }\zeta \left({x}^{-n+i}\right){\mathcal{D}}_{i}\left(t\right),|t|\le {q}^{-1}.$  Next, let us consider the double sequence ${A}_{n,r}\in K$  , ${A}_{n,1}=\left(-1{\right)}^{n-1}{L}_{n-1}$  , ${A}_{n,r}=\left(-1{\right)}^{n+r}{L}_{n-1}{\sum }_{0<{i}_{1}<\dots <{i}_{r-1}  These elements appear as the coefficients of the expansion [40of a hyperdifferentiation ${\mathcal{D}}_{r}$  in the normalized Carlitz polynomials, as well as in the expression [13of the operators ${\Delta }^{\left(n\right)}$  from (2.1) via the iterations ${\Delta }^{r}$  . Here we have the identity $\zeta \left({x}^{-n}\right)={\sum }_{i=1}^{\infty }\left(-1{\right)}^{i+1}{L}_{i}^{-1}{\sum }_{r=1}^{i}{A}_{i,r}\zeta \left({x}^{r-n}\right)$  which may be seen as a distant relative of Riemann's functional equation for the classical zeta.
Finally, consider the coefficients ${c}_{i}$  of the Carlitz expansion of ${l}_{1}$  (see Theorem 5.1). They are expressed via zeta values:
${c}_{i}={\sum }_{r=1}^{i}{A}_{i,r}\zeta \left({x}^{r-1}\right).$  By Theorem 5.1, for $i\ge 2$  we have
 $\begin{array}{c}{c}_{i}={\sum }_{j=0}^{\infty }\left({z}_{i}{\right)}^{{q}^{j}},{z}_{i}={c}_{i-1}^{q}\left[i-1{\right]}^{q}.\end{array}$ (5.5)
The series in (5.5) may be seen as an analog of ${\sum }_{j}{j}^{-z}$  . This analogy becomes clearer if, for a fixed $z\in \overline{K}c$  , $|z|<1$  , we consider the set $S$  of all convergent power series ${\sum }_{n=1}^{\infty }{z}^{{q}^{{j}_{n}}}$  corresponding to sequences $\left\{{j}_{n}\right\}\subset \mathbb{N}$  . Let us introduce the multiplication $\odot$  in $S$  setting ${z}^{{q}^{i}}\odot {z}^{{q}^{j}}={z}^{{q}^{ij}}$  and extending the operation distributively (for a similar construction in the framework of $q$  -analysis in characteristic 0 see [23). Denoting by ${{\prod }_{p}}^{\odot }$  the product in S of elements indexed by prime numbers we obtain in a standard way the identity ${c}_{i}={{\prod }_{p}}^{\odot }{\sum }_{n=0}^{\infty }\left({z}_{i}{\right)}^{{q}^{{p}^{n}}}$  (the infinite product is understood as a limit of the partial products in the topology of $\overline{K}c$  ), an analog of the Euler product formula.

6 Umbral Calculus

6.1. Classical umbral calculus [29, 28is a set of algebraic tools for obtaining, in a unified way, a rich variety of results regarding structure and properties of various polynomial sequences. There exists a lot of generalizations extending umbral methods to other classes of functions. However there is a restriction common to the whole literature on umbral calculus – the underlying field must be of zero characteristic. An attempt to mimic the characteristic zero procedures in the positive characteristic case [10revealed a number of pathological properties of the resulting structures. More importantly, these structures were not connected with the existing analysis in positive characteristic based on a completely different algebraic foundation.
A version of umbral calculus inmplementing such a connection was developed by the author [19, and we summarize it in this section. Its basic notion is motivated by the following identity for the non-normalized Carlitz polynomials ${e}_{i}={D}_{i}{f}_{i}$  :
 $\begin{array}{c}{e}_{i}\left(st\right)={\sum }_{n=0}^{i}\left(\genfrac{}{}{0}{}{i}{n}{\right)}_{K}{e}_{n}\left(t\right)\left\{{e}_{i-n}\left(s\right){\right\}}^{{q}^{n}}\end{array}$ (6.1)
where the “ $K$  -binomial coefficients” $\left(\genfrac{}{}{0}{}{i}{n}{\right)}_{K}$  are defined as $\left(\genfrac{}{}{0}{}{i}{n}{\right)}_{K}=\frac{{D}_{i}}{{D}_{n}{D}_{i-n}^{{q}^{n}}}.$  Computing the absolute values of the Carlitz factorials directly from their definition (1.1), it is easy to show that $|\left(\genfrac{}{}{0}{}{i}{n}{\right)}_{K}|=1,0\le n\le i.$  In fact, $\left(\genfrac{}{}{0}{}{i}{n}{\right)}_{K}\in \mathbb{F}q\left(x\right)$  , and we can consider also other places of $\mathbb{F}q\left(x\right)$  , that is other non-equivalent absolute values. It can be proved [22that $\left(\genfrac{}{}{0}{}{i}{n}{\right)}_{K}$  belongs to the ring of integers for any finite place of $\mathbb{F}q\left(x\right)$  .
We see the relation (6.1) as a function field counterpart of the classical binomial identity [29, 28satisfied by many classical polynomials. Now, considering a sequence ${u}_{i}$  of $\mathbb{F}q$  -linear polynomials with coefficients from $\overline{K}c$  , we call it a sequence of $K$  -binomial type if $deg{u}_{i}={q}^{i}$  and for all $i=0,1,2,\dots$
 $\begin{array}{c}{u}_{i}\left(st\right)={\sum }_{n=0}^{i}\left(\genfrac{}{}{0}{}{i}{n}{\right)}_{K}{u}_{n}\left(t\right){\left\{{u}_{i-n}\left(s\right)\right\}}^{{q}^{n}},s,t\in K.\end{array}$ (6.2)
As in the conventional umbral calculus, the dual notion is that of a delta operator. However, in contrast to the classical situation, here the delta operators are only $\mathbb{F}q$  -linear, not linear.
Denote by ${\rho }_{\lambda }$  the operator of multiplicative shift, $\left({\rho }_{\lambda }u\right)\left(t\right)=u\left(\lambda t\right)$  . We call a linear operator $T$  , on the $\overline{K}c$  -vector space $\overline{K}c\left\{t\right\}$  of all $\mathbb{F}q$  -linear polynomials, invariant if it commutes with ${\rho }_{\lambda }$  for each $\lambda \in K$  .
A $\mathbb{F}q$  -linear operator $\delta ={\tau }^{-1}{\delta }_{0}$  , where ${\delta }_{0}$  is a linear invariant operator on $\overline{K}c\left\{t\right\}$  , is called a delta operator if ${\delta }_{0}\left(t\right)=0$  and ${\delta }_{0}\left(f\right)\ne 0$  for $degf>1$  . A sequence $\left\{{P}_{n}{\right\}}_{0}^{\infty }$  of $\mathbb{F}q$  -linear polynomials is called a basic sequence corresponding to a delta operator $\delta ={\tau }^{-1}{\delta }_{0}$  , if $deg{P}_{n}={q}^{n}$  , ${P}_{0}\left(1\right)=1$  , ${P}_{n}\left(1\right)=0$  for $n\ge 1$  ,
 $\begin{array}{c}\delta {P}_{0}=0,\delta {P}_{n}=\left[n{\right]}^{1/q}{P}_{n-1},n\ge 1,\end{array}$ (6.3)
or, equivalently,
 $\begin{array}{c}{\delta }_{0}{P}_{0}=0,{\delta }_{0}{P}_{n}=\left[n\right]{P}_{n-1}^{q},n\ge 1.\end{array}$ (6.4)
It is clear that $d={\tau }^{-1}\Delta$  is a delta operator. It follows from well-known identities for the Carlitz polynomials ${e}_{i}$  [11(see also (1.10)) that the sequence $\left\{{e}_{i}\right\}$  is basic with respect to the operator $d$  .
Theorem 6.1. For any delta operator $\delta ={\tau }^{-1}{\delta }_{0}$  , there exists a unique basic sequence $\left\{{P}_{n}\right\}$  , which is a sequence of $K$  -binomial type. Conversely, given a sequence $\left\{{P}_{n}\right\}$  of $K$  -binomial type, define the action of ${\delta }_{0}$  on ${P}_{n}$  by the relations (6.4), extend it onto $\overline{K}c\left\{t\right\}$  by linearity and set $\delta ={\tau }^{-1}{\delta }_{0}$  . Then $\delta$  is a delta operator, and $\left\{{P}_{n}\right\}$  is the corresponding basic sequence.
The analogs of the higher Carlitz difference operators (2.1) in the present general context are the operators ${\delta }_{0}^{\left(l\right)}={\tau }^{l}{\delta }^{l}$  . The identity
 $\begin{array}{c}{\delta }_{0}^{\left(l\right)}{P}_{j}=\frac{{D}_{j}}{{D}_{j-l}^{{q}^{l}}}{P}_{j-l}^{{q}^{l}}\end{array}$ (6.5)
holds for any $l\le j$  . If $f$  is a $\mathbb{F}q$  -linear polynomial, $degf\le {q}^{n}$  , then a generalized Taylor formula
 $\begin{array}{c}f\left(st\right)={\sum }_{l=0}^{n}\frac{\left({\delta }_{0}^{\left(l\right)}f\right)\left(s\right)}{{D}_{l}}{P}_{l}\left(t\right)\end{array}$ (6.6)
holds for any $s,t\in K$  . For the Carlitz polynomials ${e}_{i}$  , the formulas (6.5) and (6.6) are well known [11. It is important that, in contrast to the classical umbral calculus, the linear operators involved in (6.6) are not powers of a single linear operator.
Any linear invariant operator $T$  on $\overline{K}c\left\{t\right\}$  admits a representation
 $\begin{array}{c}T={\sum }_{l=0}^{\infty }{\sigma }_{l}{\delta }_{0}^{\left(l\right)},{\sigma }_{l}=\frac{\left(T{P}_{l}\right)\left(1\right)}{{D}_{l}}.\end{array}$ (6.7)
The infinite series in (6.7) becomes actually a finite sum if both sides of (6.7) are applied to any $\mathbb{F}q$  -linear polynomial. Conversely, any such series defines a linear invariant operator on $\overline{K}c\left\{t\right\}$  .
Let us consider the case where $\delta =d$  , so that ${\delta }_{0}^{\left(l\right)}={\Delta }^{\left(l\right)}$  . The next result leads to new delta operators and basic sequences.
Theorem 6.2. The operator $\theta ={\tau }^{-1}{\theta }_{0}$  , where ${\theta }_{0}={\sum }_{l=1}^{\infty }{\sigma }_{l}{\Delta }^{\left(l\right)},$  is a delta operator if and only if
 $\begin{array}{c}{S}_{n}\stackrel{\text{def}}{=}{\sum }_{l=1}^{n}\frac{{\sigma }_{l}}{{D}_{n-l}^{{q}^{l}}}\ne 0\text{for all}n=1,2,\dots .\end{array}$ (6.8)
Example 1. Let ${\sigma }_{l}=1$  for all $l\ge 1$  , that is
 $\begin{array}{c}{\theta }_{0}={\sum }_{l=1}^{\infty }{\Delta }^{\left(l\right)}.\end{array}$ (6.9)
Estimates of $|{D}_{n}|$  which follow directly from (1.1) show that $|{S}_{n}|={q}^{\frac{{q}^{n}-q}{q-1}}$  , so that (6.8) is satisfied. Comparing (6.9) with a classical formula from [29we may see the polynomials ${P}_{n}$  for this case as analogs of the Laguerre polynomials.
Example 2. Let ${\sigma }_{l}=\frac{\left(-1{\right)}^{l+1}}{{L}_{l}}$  . For this case it can be shown [19that ${S}_{n}={D}_{n}^{-1}$  , $n=1,2,\dots$  ; ${\theta }_{0}\left({t}^{{q}^{j}}\right)={t}^{{q}^{j}}$  for all $j\ge 1$  (of course, ${\theta }_{0}\left(t\right)=0$  ), and ${P}_{0}\left(t\right)=t$  , ${P}_{n}\left(t\right)={D}_{n}\left({t}^{{q}^{n}}-{t}^{{q}^{n-1}}\right)$  for $n\ge 1$  .
6.2. As in the $p$  -adic case [38, 39, 27, the umbral calculus can be used for constructing new orthonormal bases in ${C}_{0}\left(O,\overline{K}c\right)$  .
Let $\left\{{P}_{n}\right\}$  be the basic sequence corresponding to a delta operator $\delta ={\tau }^{-1}{\delta }_{0}$  ,
 $\begin{array}{c}{\delta }_{0}={\sum }_{l=1}^{\infty }{\sigma }_{l}{\Delta }^{\left(l\right)}.\end{array}$ (6.10)
The sequence ${Q}_{n}=\frac{{P}_{n}}{{D}_{n}}$  , $n=0,1,2,\dots$  , called the normalized basic sequence, satisfies the identity ${Q}_{i}\left(st\right)={\sum }_{n=0}^{i}{Q}_{n}\left(t\right){\left\{{Q}_{i-n}\left(s\right)\right\}}^{{q}^{n}},$  another form of the $K$  -binomial property. Though it resembles its classical counterpart, the presence of the Frobenius powers is a feature specific for the case of a positive characteristic.
Theorem 6.3. If $|{\sigma }_{1}|=1$  , $|{\sigma }_{l}|\le 1$  for $l\ge 2$  , then the sequence $\left\{{Q}_{n}{\right\}}_{0}^{\infty }$  is an orthonormal basis of the space ${C}_{0}\left(O,\overline{K}c\right)$  – for any $f\in {C}_{0}\left(O,\overline{K}c\right)$  there is a uniformly convergent expansion $f\left(t\right)={\sum }_{n=0}^{\infty }{\psi }_{n}{Q}_{n}\left(t\right),t\in O,$  where ${\psi }_{n}=\left({\delta }_{0}^{\left(n\right)}f\right)\left(1\right)$  , $|{\psi }_{n}|\to 0$  as $n\to \infty$  , $\parallel f\parallel ={sup}_{n\ge 0}|{\psi }_{n}|.$
By Theorem 6.3, the Laguerre-type polynomial sequence from Example 1 is an orthonormal basis of ${C}_{0}\left(O,\overline{K}c\right)$  . The sequence from Example 2 does not satisfy the conditions of Theorem 6.3.
Note that the conditions of Theorem 6.3 imply that ${S}_{n}\ne 0$  for all $n$  , so that the series (6.10) considered in Theorem 6.3 always correspond to delta operators.
In [19recursive formulas and generating functions for normalized basic sequences are also given.

7 The Weyl-Carlitz Ring and Holonomic Modules

7.1. The theory of holonomic modules over the Weyl algebra and more general algebras of differential or $q$  -difference operators is becoming increasingly important, both as a crucial part of the general theory of D-modules and in view of various applications (see, for example, [2, 5, 30).
It is clear from the above results that in the positive characteristic case a natural counterpart of the Weyl algebra is, for the case of a single variable, the ring ${\mathfrak{A}}_{1}$  generated by $\tau ,d$  , and scalars from $\overline{K}c$  , with the relations
 $\begin{array}{c}d\tau -\tau d=\left[1{\right]}^{1/q},\tau \lambda ={\lambda }^{q}\tau ,d\lambda ={\lambda }^{1/q}d\left(\lambda \in \overline{K}c\right).\end{array}$ (7.1)
The ring consists of finite sums
 $\begin{array}{c}a={\sum }_{i,j}{\lambda }_{ij}{\tau }^{i}{d}^{j},{\lambda }_{ij}\in \overline{K}c,\end{array}$ (7.2)
and the representation of an element in the form (7.2) is unique.
The next result proved in [17shows that basic algebraic properties of ${\mathfrak{A}}_{1}$  are similar to those of the Weyl algebra in characteristic 0 and quite different from the case of the algebra of usual differential operators over a field of positive characteristic [26.
Theorem 7.1. The ring ${\mathfrak{A}}_{1}$  is left and right Noetherian, without zero divisors; its centre coincides with $\mathbb{F}q$  . ${\mathfrak{A}}_{1}$  possesses no non-trivial two-sided ideals stable with respect to the mapping ${\sum }_{i,j}{\lambda }_{ij}{\tau }^{i}{d}^{j}↦{\sum }_{i,j}{\lambda }_{ij}^{q}{\tau }^{i}{d}^{j}.$
A generalization of ${\mathfrak{A}}_{1}$  to the case of several variables is not straightforward because the Carlitz derivatives ${d}_{s}$  and ${d}_{t}$  do not commute on a monomial $f\left(s,t\right)={s}^{{q}^{m}}{t}^{{q}^{n}}$  , if $m\ne n$  .
Moreover, if $m>n$  , then ${d}_{s}^{m}f$  is not a polynomial, nor even a holomorphic function in $t$  (since the action of $d$  is not linear and involves taking the $q$  -th root).
A reasonable generalization is inspired by Zeilberger's idea (see [5) to study holonomic properties of sequences of functions making a transform with respect to the discrete variables, which reduces the continuous-discrete case to the purely continuous one (simultaneously in all the variables). In our situation, if $\left\{{P}_{k}\left(s\right)\right\}$  is a sequence of $\mathbb{F}q$  -linear polynomials with $deg{P}_{k}\le {q}^{k}$  , we set
 $\begin{array}{c}f\left(s,t\right)={\sum }_{k=0}^{\infty }{P}_{k}\left(s\right){t}^{{q}^{k}},\end{array}$ (7.3)
and ${d}_{s}$  is well-defined. In the variable $t$  , we consider not ${d}_{t}$  but the linear operator ${\Delta }_{t}$  . The latter does not commute with ${d}_{s}$  either, but satisfies the commutation relations ${d}_{s}{\Delta }_{t}-{\Delta }_{t}{d}_{s}=\left[1{\right]}^{1/q}{d}_{s},{\Delta }_{t}\tau -\tau {\Delta }_{t}=\left[1\right]\tau ,$  so that the resulting ring ${\mathfrak{A}}_{2}$  resembles a universal enveloping algebra of a solvable Lie algebra.
More generally, denote by ${\mathcal{ℱ}}_{n+1}$  the set of all germs of functions of the form
 $\begin{array}{c}f\left(s,{t}_{1},\dots ,{t}_{n}\right)={\sum }_{{k}_{1}=0}^{\infty }\dots {\sum }_{{k}_{n}=0}^{\infty }{\sum }_{m=0}^{min\left({k}_{1},\dots ,{k}_{n}\right)}{a}_{m,{k}_{1},\dots ,{k}_{n}}{s}^{{q}^{m}}{t}_{1}^{{q}^{{k}_{1}}}\dots {t}_{n}^{{q}^{{k}_{n}}}\end{array}$ (7.4)
where ${a}_{m,{k}_{1},\dots ,{k}_{n}}\in \overline{K}c$  are such that all the series are convergent on some neighbourhoods of the origin. We do not exclude the case $n=0$  where ${\mathcal{ℱ}}_{1}$  will mean the set of all $\mathbb{F}q$  -linear power series ${\sum }_{m}{a}_{m}{s}^{{q}^{m}}$  convergent on a neighbourhood of the origin. ${\stackrel{^}{\mathcal{ℱ}}}_{n+1}$  will denote the set of all polynomials from ${\mathcal{ℱ}}_{n+1}$  , that is the series (7.4) in which only a finite number of coefficients is different from zero.
The ring ${\mathfrak{A}}_{n+1}$  is generated by the operators $\tau ,{d}_{s},{\Delta }_{{t}_{1}},\dots {\Delta }_{{t}_{n}}$  on ${\mathcal{ℱ}}_{n+1}$  , and the operators of multiplication by scalars from $\overline{K}c$  . To simplify the notation, we write ${\Delta }_{j}$  instead of ${\Delta }_{{t}_{j}}$  and identify a scalar $\lambda \in \overline{K}c$  with the operator of multiplication by $\lambda$  . The operators ${\Delta }_{j}$  are $\overline{K}c$  -linear, so that
 $\begin{array}{c}{\Delta }_{j}\lambda =\lambda {\Delta }_{j},\lambda \in \overline{K}c,\end{array}$ (7.5)
while the operators $\tau ,{d}_{s}$  satisfy the commutation relations (7.1). In the action of each operator ${d}_{s},{\Delta }_{j}$  (acting in a single variable), other variables are treated as scalars. The operator $\tau$  acts simultaneously on all the variables and coefficients. We have the relations involving ${\Delta }_{j}$  :
 $\begin{array}{c}{\Delta }_{j}\tau -\tau {\Delta }_{j}=\left[1\right]\tau ,{d}_{s}{\Delta }_{j}-{\Delta }_{j}{d}_{s}=\left[1{\right]}^{1/q}{d}_{s},j=1,\dots ,n.\end{array}$ (7.6)
Using the commutation relations (7.1), (7.5), and (7.6), we can write any $a\in {\mathfrak{A}}_{n+1}$  , in a unique way, as a finite sum
 $\begin{array}{c}a=\sum {c}_{l,\mu ,{i}_{1},\dots ,{i}_{n}}{\tau }^{l}{d}_{s}^{\mu }{\Delta }_{1}^{{i}_{1}}\dots {\Delta }_{n}^{{i}_{n}}.\end{array}$ (7.7)
Let us introduce a filtration in ${\mathfrak{A}}_{n+1}$  (an analog of the Bernstein filtration) denoting by ${\Gamma }_{\nu }$  , $\nu \in {\mathbb{Z}}_{+}$  , the $\overline{K}c$  -vector space of operators (7.7) with $max\left\{l+\mu +{i}_{1}+\cdots +{i}_{n}\right\}\le \nu$  where the maximum is taken over all the terms of (7.7). Then ${\mathfrak{A}}_{n+1}$  is a left and right Noetherian filtered ring.
In a standard way (see [8) we define filtered left modules over ${\mathfrak{A}}_{n+1}$  . All the basic notions regarding a filtered module $M$  (like those of the graded module $gr\left(M\right)$  , dimension $d\left(M\right)$  , multiplicity $m\left(M\right)$  , good filtration etc) are introduced just as their counterparts in the theory of modules over the Weyl algebra.
If we consider ${\mathfrak{A}}_{n+1}$  as a left module over itself, then
 $\begin{array}{c}d\left({\mathfrak{A}}_{n+1}\right)=n+2,m\left({\mathfrak{A}}_{n+1}\right)=1.\end{array}$ (7.8)
For any finitely generated left ${\mathfrak{A}}_{n+1}$  -module $M$  , we have $d\left(M\right)\le n+2$  . By (7.8), this bound cannot be improved in general. However, if $I$  is a non-zero left ideal in ${\mathfrak{A}}_{n+1}$  , then
 $\begin{array}{c}d\left({\mathfrak{A}}_{n+1}/I\right)\le n+1.\end{array}$ (7.9)
For the module ${\stackrel{^}{\mathcal{ℱ}}}_{n+1}$  of $\mathbb{F}q$  -linear polynomials (7.4), we have $d\left({\stackrel{^}{\mathcal{ℱ}}}_{n+1}\right)=n+1,m\left({\stackrel{^}{\mathcal{ℱ}}}_{n+1}\right)=n!$  The proofs of all these results, as well as the ones given in this section below, can be found in [22.
It is natural to call an ${\mathfrak{A}}_{n+1}$  -module $M$  holonomic if $d\left(M\right)=n+1$  . Thus, ${\stackrel{^}{\mathcal{ℱ}}}_{n+1}$  is an example of a holonomic module.
The next theorem demonstrates, already for the case of ${\mathfrak{A}}_{1}$  -modules, a sharp difference from the case of modules over the Weyl algebras. In particular, we see that an analog of the Bernstein inequality (see [8) does not hold here without some additional assumptions.
Theorem 7.2.
• (i) For any $k=1,2,\dots$  , there exists such a nontrivial ${\mathfrak{A}}_{1}$  -module $M$  that $dimM=k$  ( $dim$  means the dimension over $\overline{K}c$  ), that is $d\left(M\right)=0$  .
• (ii) Let $M$  be a finitely generated ${\mathfrak{A}}_{1}$  -module with a good filtration. Suppose that there exists a “vacuum vector” $v\in M$  , such that ${d}_{s}v=0$  and ${\tau }^{m}\left(v\right)\ne 0$  for all $m=0,1,2,\dots$  . Then $d\left(M\right)\ge 1$  .
7.3. Let us consider the case of holonomic submodules of the ${\mathfrak{A}}_{n+1}$  -module ${\mathcal{ℱ}}_{n+1}$  , consisting of $\mathbb{F}q$  -linear functions (7.4) polynomial in $s$  and holomorphic near the origin in ${t}_{1},\dots ,{t}_{n}$  .
Let $0\ne f\in {\mathcal{ℱ}}_{n+1}$  , ${I}_{f}=\left\{\phi \in {\mathfrak{A}}_{n+1}:\phi \left(f\right)=0\right\}.$  ${I}_{f}$  is a left ideal in ${\mathfrak{A}}_{n+1}$  . The left ${\mathfrak{A}}_{n+1}$  -module ${M}_{f}={\mathfrak{A}}_{n+1}/{I}_{f}$  is isomorphic to the submodule ${\mathfrak{A}}_{n+1}f\subset {\mathcal{ℱ}}_{n+1}$  – an element $\phi \left(f\right)\in {\mathfrak{A}}_{n+1}f$  corresponds to the class of $\phi \in {\mathfrak{A}}_{n+1}$  in ${M}_{f}$  . A natural good filtration in ${M}_{f}$  is induced from that in ${\mathfrak{A}}_{n+1}$  .
As we know (see (7.9)), if ${I}_{f}\ne \left\{0\right\}$  , then $d\left({M}_{f}\right)\le n+1$  . We call a function $f$  holonomic if the module ${M}_{f}$  is holonomic, that is $d\left({M}_{f}\right)=n+1$  . The condition ${I}_{f}\ne \left\{0\right\}$  means that $f$  is a solution of a non-trivial “differential equation” $\phi \left(f\right)=0$  , $\phi \in {\mathfrak{A}}_{n+1}$  . The case $n=0$  is quite simple.
Theorem 7.3. If a non-zero function $f\in {\mathcal{ℱ}}_{1}$  satisfies an equation $\phi \left(f\right)=0$  , $0\ne \phi \in {\mathfrak{A}}_{1}$  , then $f$  is holonomic.
In particular, any $\mathbb{F}q$  -linear polynomial of $s$  is holonomic, since it is annihilated by ${d}_{s}^{m}$  , with a sufficiently large $m$  .
If $n>0$  , the situation is more complicated. We call the module ${M}_{f}$  (and the corresponding function $f$  ) degenerate if $D\left({M}_{f}\right)  (by the Bernstein inequality, there is no degeneracy phenomena for modules over the complex Weyl algebra). The simplest example of a degenerate function (for $n=1$  ) is $f\left(s,{t}_{1}\right)=g\left(s{t}_{1}\right)\in {\mathcal{ℱ}}_{2}$  where the function $g$  belongs to ${\mathcal{ℱ}}_{1}$  and satisfies an equation $\phi \left(g\right)=0$  , $\phi \in {\mathfrak{A}}_{1}$  . It can be shown that $d\left({M}_{f}\right)=1$  .
In order to exclude the degenerate case, we introduce the notion of a non-sparse function.
A function $f\in {\mathcal{ℱ}}_{n+1}$  of the form (7.4) is called non-sparse if there exists such a sequence ${m}_{l}\to \infty$  that, for any $l$  , there exist sequences ${k}_{1}^{\left(i\right)},{k}_{2}^{\left(i\right)},\dots ,{k}_{n}^{\left(i\right)}\ge {m}_{l}$  (depending on $l$  ), such that ${k}_{\nu }^{\left(i\right)}\to \infty$  as $i\to \infty$  ( $\nu =1,\dots ,n$  ), and ${a}_{m,{k}_{1}^{\left(i\right)},\dots ,{k}_{n}^{\left(i\right)}}\ne 0$  .
Theorem 7.4. If a function $f$  is non-sparse, then $d\left({M}_{f}\right)\ge n+1$  . If, in addition, $f$  satisfies an equation $\phi \left(f\right)=0$  , $0\ne \phi \in {\mathfrak{A}}_{n+1}$  , then $f$  is holonomic.
7.4 We use Theorem 7.4 to prove that the functions (7.4) obtained via the sequence-to-function transform (7.3) or its multi-index generalizations, from some well-known sequences of polynomials over $K$  are holonomic. In all the cases below the non-sparseness is evident, and we have only to prove that the corresponding function satisfies a non-trivial Carlitz differential equation.
a) The Carlitz polynomials. The transform (7.3) of the sequence $\left\{{f}_{k}\right\}$  is the Carlitz module function ${C}_{s}\left(t\right)$  ; see (2.3). It is easy to check that ${d}_{s}{C}_{s}\left(t\right)={C}_{s}\left(t\right)$  . Therefore the Carlitz module function is holonomic, jointly in both its variables.
b) Thakur's hypergeometric polynomials. We consider the polynomial case of Thakur's hypergeometric function (1.5), that is
 $\begin{array}{c}{}_{l}{F}_{\lambda }\left(-{a}_{1},\dots ,-{a}_{l};-{b}_{1},\dots ,-{b}_{\lambda };z\right)={\sum }_{m}\frac{\left(-{a}_{1}{\right)}_{m}\dots \left(-{a}_{l}{\right)}_{m}}{\left(-{b}_{1}{\right)}_{m}\dots \left(-{b}_{\lambda }{\right)}_{m}{D}_{m}}{z}^{{q}^{m}}\end{array}$ (7.10)
where ${a}_{1},\dots ,{a}_{l},{b}_{1},\dots ,{b}_{\lambda }\in {\mathbb{Z}}_{+}$  . It is seen from (1.4) that the terms in (7.10), which make sense and do not vanish, are those with $m\le min\left({a}_{1},\dots ,{a}_{l},{b}_{1},\dots ,{b}_{\lambda }\right)$  . Let the function $f\in {\mathcal{ℱ}}_{l+\lambda +1}$  be given by
$f\left(s,{t}_{1},\dots ,{t}_{l},{u}_{1},\dots ,{u}_{\lambda }\right)={\sum }_{{k}_{1}=0}^{\infty }\dots {\sum }_{{k}_{l}=0}^{\infty }{\sum }_{{\nu }_{1}=0}^{\infty }\dots {{\sum }_{{\nu }_{\lambda }=0}^{\infty }}_{l}{F}_{\lambda }\left(-{k}_{1},\dots ,-{k}_{l};-{\nu }_{1},\dots ,-{\nu }_{\lambda };s\right){t}_{1}^{{q}^{{k}_{1}}}\dots {t}_{l}^{{q}^{{k}_{l}}}{u}_{1}^{{q}^{{\nu }_{1}}}\dots {u}_{\lambda }^{{q}^{{\nu }_{\lambda }}}.$
It is known ([37, Sect. 6.5) that
 $\begin{array}{c}{{d}_{s}}_{l}{F}_{\lambda }\left(-{k}_{1},\dots ,-{k}_{l};-{\nu }_{1},\dots ,-{\nu }_{\lambda };s\right){=}_{l}{F}_{\lambda }\left(-{k}_{1}+1,\dots ,-{k}_{l}+1;-{\nu }_{1}+1,\dots ,-{\nu }_{\lambda }+1;s\right)\end{array}$ (7.11)
if all the parameters ${k}_{1},\dots ,{k}_{l},{\nu }_{1},\dots ,{\nu }_{\lambda }$  are different from zero. If at least one of them is equal to zero, then the left-hand side of (7.11) equals zero. This property implies the identity ${d}_{s}f=f$  , the same as that for the Carlitz module function. Thus, $f$  is holonomic.
c). $K$  -binomial coefficients. It can be shown [22that the $K$  -binomial coefficients $\left(\genfrac{}{}{0}{}{k}{m}{\right)}_{K}$  (see Sect. 6) satisfy the Pascal-type identity
 $\begin{array}{c}\left(\genfrac{}{}{0}{}{k}{m}{\right)}_{K}=\left(\genfrac{}{}{0}{}{k-1}{m-1}{\right)}_{K}^{q}+\left(\genfrac{}{}{0}{}{k-1}{m}{\right)}_{K}^{q}{D}_{m}^{q-1}\end{array}$ (7.12)
where $0\le m\le k$  and it is assumed that $\left(\genfrac{}{}{0}{}{k}{-1}{\right)}_{K}=\left(\genfrac{}{}{0}{}{k-1}{k}{\right)}_{K}=0$  .
Consider a function $f\in {\mathcal{ℱ}}_{2}$  associated with the $K$  -binomial coefficients, that is
 $\begin{array}{c}f\left(s,t\right)={\sum }_{k=0}^{\infty }{\sum }_{m=0}^{k}\left(\genfrac{}{}{0}{}{k}{m}{\right)}_{K}{s}^{{q}^{m}}{t}^{{q}^{k}}.\end{array}$ (7.13)
The identity (7.12) implies the equation ${d}_{s}f\left(s,t\right)={\Delta }_{t}f\left(s,t\right)+\left[1{\right]}^{1/q}f\left(s,t\right)$  for the function (7.13). Therefore $f$  is holonomic.
Acknowledgements This work was supported in part by CRDF (Grant UM1-2567-OD-03), DFG (Grant 436 UKR 113/72), and the Ukrainian Foundation for Fundamental Research (Grant 01.07/027).
References

1. F. Baldassarri, Differential modules and singular points of $p$  -adic differential equations, Adv. Math. 44 (1982), 155–179.
2. Yu. Berest and A. Kasman, $\mathcal{D}$  -modules and Darboux transformations, Lett. Math. Phys. 43 (1998), 279–294.
3. L. Carlitz, On certain functions connected with polynomials in a Galois field, Duke Math. J. 1 (1935), 137–168.
4. L. Carlitz, Some special functions over $GF\left(q,x\right)$  , Duke Math. J. 27 (1960), 139–158.
5. P. Cartier, Démonstration “automatique” d'identités et fonctions hypergéometriques (d'après D. Zeilberger), Astérisque 206 (1992), 41–91.
6. D. N. Clark, A note on the $p$  -adic convergence of solutions of linear differential equations, Proc. Amer. Math. Soc. 17 (1966), 262–269.
7. K. Conrad, The digit principle, J. Number Theory 84 (2000), 230–237.
8. S. C. Coutinho, A Primer of Algebraic D-modules, Cambridge University Press, 1995.
9. B. Dwork, G. Gerotto, and F. J. Sullivan, An Introduction to $G$  -Functions, Princeton University Press, 1994.
10. L. Ferrari, An umbral calculus over infinite coefficient fields of positive characteristic, Comp. Math. Appl. 41 (2001), 1099–1108.
11. D. Goss, Fourier series, measures, and divided power series in the theory of function fields, K-Theory 1 (1989), 533–555.
12. D. Goss, Basic Structures of Function Field Arithmetic, Springer, Berlin, 1996.
13. S. Jeong, Continuous linear endomorphisms and difference equations over the completions of $\mathbb{F}q\left[T\right]$  , J. Number Theory 84 (2000), 276–291.
14. S. Jeong, Hyperdifferential operators and continuous functions on function fields, J. Number Theory 89 (2001), 165–178.
15. A. N. Kochubei, Harmonic oscillator in characteristic $p$  , Lett. Math. Phys. 45 (1998), 11–20.
16. A. N. Kochubei, $\mathbb{F}q$  -linear calculus over function fields, J. Number Theory 76 (1999), 281–300.
17. A. N. Kochubei, Differential equations for $\mathbb{F}q$  -linear functions, J. Number Theory 83 (2000), 137–154.
18. A. N. Kochubei, Differential equations for $\mathbb{F}q$  -linear functions II: Regular singularity, Finite Fields Appl. 9 (2003), 250–266.
19. A. N. Kochubei, Umbral calculus in positive characteristic, Adv. Appl. Math. 34 (2005), 175–191.
20. A. N. Kochubei, Strongly nonlinear differential equations with Carlitz derivatives over a function field, Ukrainian Math. J. (to appear), math.NT/ 0405542.
21. A. N. Kochubei, Polylogarithms and a zeta function for finite places of a function field, Contemporary Math. (to appear); math.NT/0405544.
22. A. N. Kochubei, Holonomic modules in positive characteristic, math.RA/0503398.
23. M. B. Nathanson, Additive number theory and the ring of quantum integers, math.NT/0204006.
24. B. Poonen, Fractional power series and pairings on Drinfeld modules, J. Amer. Math. Soc. 9 (1996), 783–812.
25. M. van der Put, Meromorphic differential equations over valued fields, Indag. Math. 42 (1980), 327–332.
26. M. van der Put, Differential equations in characteristic $p$  , Compositio Math. 97 (1995), 227–251.
27. A. M. Robert, A Course in $p$  -Adic Analysis, Springer, New York, 2000.
28. S. Roman, The Umbral Calculus, Academic Press, London, 1984.
29. G.-C. Rota, D. Kahaner and A. Odlyzko, On the foundations of combinatorial theory. VIII. Finite operator calculus, J. Math. Anal. Appl. 42 (1973), 684–760.
30. C. Sabbah, Systèmes holonomes d'équations aux q-differences. In: D-modules and Microlocal Geometry (M. Kashiwara et al., eds.), Walter de Gruyter, Berlin, 1993, pp. 125–147.
31. S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives: Theory and Applications, Gordon and Breach, New York, 1993.
32. W. Schikhof, Ultrametric Calculus, Cambridge University Press, 1984.
33. M. Setoyanagi, Note on Clark's theorem for $p$  -adic convergence, Proc. Amer. Math. Soc. 125 (1997), 717–721.
34. D. Sinnou and D. Laurent, Indépendence algebrique sur les $T$  -modules, Compositio Math. 122 (2000), 1–22.
35. D. S. Thakur, Hypergeometric functions for function fields, Finite Fields and Their Appl. 1 (1995), 219–231.
36. D. S. Thakur, Hypergeometric functions for function fields II, J. Ramanujan Math. Soc. 15 (2000), 43–52.
37. D. S. Thakur, Function Field Arithmetic, World Scientific, Singapore, 2004.
38. L. Van Hamme, Continuous operators, which commute with translations, on the space of continuous functions on ${\mathbb{Z}}_{p}$  . In: $p$  -Adic Functional Analysis (J. M. Bayod et al., eds.), Lect. Notes Pure Appl. Math. 137, Marcel Dekker, New York, 1992, pp. 75–88.
39. A. Verdoodt, Umbral calculus in non-Archimedean analysis. In: $p$  -Adic Functional Analysis (A. K. Katsaras et al., eds.), Lect. Notes Pure Appl. Math. 222, Marcel Dekker, New York, 2001, pp. 309–322.
40. J. F. Voloch, Differential operators and interpolation series in power series fields, J. Number Theory 71 (1998), 106–108.
41. C. G. Wagner, Linear operators in local fields of prime characteristic, J. Reine Angew. Math. 251 (1971), 153–160.
42. Z. Yang, Locally analytic functions over completions of ${\mathbf{F}}_{r}\left(U\right)$  , J. Number Theory 73 (1998), 451–458.
43. Z. Yang, ${C}^{n}$  -functions over completions of ${\mathbb{F}}_{r}\left[T\right]$  at finite places of ${\mathbb{F}}_{r}\left(T\right)$  , J. Number Theory 108 (2004), 346–374.