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:
Then changing variables gives
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,
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.
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,
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 ,
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
From the Chu-Vandermonde identity,
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 contour 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 interpolation 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 .
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: