The background info and comments for the MSE question Lie group heuristics for a raising operator for and the MO question Riemann zeta function at positive integers and an Appell sequence of poylnomials introduce an Appell sequence of polynomials containing the Euler-Mascheroni constant and the Riemann zeta function evaluated at the integers greater than one. The Appell sequence can be defined by its e.g.f.

The raising op for the Appell sequence

where is the digamma or Psi function, is associated with the infinitesimal generator

of a class of fractional integro-derivatives through a change of variables (). Exponentiation gives

As discussed in previous entries here on Appell sequences, every Appell sequence of polynomials has an umbral compositional inverse sequence of polynomials, which is also an Appell sequence whose raising op involves only a simple change in a sign in the formula for the inverse sequence’s raising op. The e.g.f. of the umbral inverse sequence in this case is

with the raising op

The associated raising op through the transformation is

which, when exponentiated, gives

where (cf. Pochhammer symbol of a differential and hypergeometric polynomials for definitions of the symbols).

These relations are simple to derive given the property defining the raising ops:

so

Then changing variables gives

implying

or

As a sanity check, confirm

The relation to can be seen from the derivative

implying, from the notes attached to the MSE and MO questions,

where

with . (Note the integrand is related to the q (Jackson) derivative, and the Pincherle derivative / commutator is .)

Additionally, using well-known formulas for the digamma function with ,

(The last expression is consistent with the derivative of the formulas for above.)

Explicit formulas for the first few Appell polynomials associated with are

,

where , the Euler-Mascheroni constant, and is the Riemann zeta function.

For

See the MO question for more notes on the polynomials and a derivation of the recursion relation. See the entry here on Appell polynomials and cumulants for a general derivation of a recursion using the raising operator that is satisfied by all Appell sequences.

The sequence of polynomials associated to are given by negating the Euler-Mascheroni constant and the zeta values in the since the raising op has the same form as that for the cycle index partition polynomials for the symmetric groups (cf. the MO question and OEIS A036039):

the first few are

,

They comprise the umbral compositional inverses of that is,

For example,

**Relation to the Mellin transform, convolutions, and umbral calculus:**

Given the Mellin transform of as and using the Mellin transform convolution theorem,

where, using the Euler reflection formula for factorials,

and, with the Heaviside step function, for ,

so

where now can be extended beyond those functions that have Mellin transforms. Moreover, the last convolution integral can be analytically continued to positive values of as a Hadamard finite part integral for fractional differintegration that formally acts on a function analytic at that has been dilated by and :

Summarizing, the Mellin convolution theorem leads to the fractional differintegral

with the umbral compositional rep

The generalized shift operator gives a diff op rep for the umbral composition as

with the definition

or

with

From the Chu-Vandermonde identity,

so

in agreement with the other derivations above.

Note that when the diff op is acting on a function analytic about the origin, the simpler rep applies

but acting on the more general gives a Newton series interpolation, regenerating the generalized Chu-Vandermonde identity:

**The infinigen from a Cauchy c****ontour integral rep:**

is satisfied for analytic on the interval of integration with no finite part required for . For other values of , whichever integral is convergent can be taken as the analytic continuation of the divergent one. Taking the derivative of the Cauchy contour formula gives

and collapsing the contour about the real axis gives

**Consistency with Mellin i****nterpolation of operators:**

Probably the mathematician to most ingeniously wield the power of the Mellin transform in some sense as an interpolation of the coefficients of exponential generating functions or Taylor series was Ramanujan with his Master Formula / Theorem. The iconic Euler integral for the gamma function for

provides the scaffolding for similarly interpolating operators. Let’s first show that a natural interpolation of integer powers of the derivative is the fractional integroderivative of fractional calculus by using the Mellin transform to interpolate the op coefficients of the op e.g.f. i.e., the shift op, for the integer powers of the derivative:

Then specifically acting on the power function for

The last summation converges with no restriction on . So, we see that the Mellin transform does indeed interpolate the coefficients of the e.g.f. generated by the binomial theorem expansion to to give an interpolation of the coefficients of the shift op to consistent with fractional calculus.

Next let’s use the same method to interpolate

where denotes the Laguerre polynomials, and, again, by definition. Following the algorithm above, we have from the generating function for the Laguerre polynomials the action of the generalized shift op

for , where is a convergence factor approaching unity necessary when . We can anticipate that this sum converges to an interpolation of the coefficients of the e.g.f. to in agreement with the interpretation above of . With these interpolation heuristics, the summation formula can be verified by using the inverse Mellin transform and the singularities of the factorials.

So, analytically continued from ,

To confirm the presumed Mellin transform evaluation from the interpolation heuristic applied either to the real coefficients of the e.g.f. in the integrand of the transform or the operator e.g.f. do the inverse Mellin transform. The modified MT pair is

and our presumed MT result gives

so, for , the inverse Mellin transform to evaluate is

Closing the first contour to the left for picks up the singularities at , and closing the second contour to the right for picks up a sign for the clockwise path transversal and the singularities at , giving

since , confirming that the presumed solution is indeed correct.

On a more general note, the algorithm can be applied to differential operators encompassing the Witt Lie algebra to obtain

For , the Heaviside step function in the integral is to be interpreted as setting the upper limit of the integral as the first zero from the origin of , i.e., the upper limit should be .

**Related stuff:**

For Appell sequences:

“Bernoulli Appells” entry below.

For differential ops:

“Diff. Ops, Special Polynomials, Binomial Transforms, Interpolation, and the Inverse Mellin Transform” entry below,

“Goin’ with the Flow: Log of the Derivative” entry below,

“Infinigens, the Pascal Triangle, and the Witt and Virasoro Algebra” entry below,

“A Generalized Dobinski Relation and the Confluent Hypergeometric Functions” entry below,

“Lagrange a la Lah: Part I” entry below.

For Ramanujan’s Master Theorem / Formula:

“Ramanajun’s Master Theorem and Duality of Symmetric Spaces” by Bertram