Appell Polynomials in Shifted Asymptotic Expansions: the Mills ratio, Hermite polynomials, and Stieltjes bounds

A theorem of Burić, Elezović and Vuk\v sić states that translating the argument in an asymptotic expansion, $f(x)\sim\sum(-1)^na_nx^{-n-1}$, to $f(x+t)$, replaces the constant coefficients $a_n$ by the Appell polynomials $R_n(t)$ generated by $(a_n)$.

Their motivating examples came from the gamma and polygamma functions, where Bernoulli polynomials occur.

We extend this construction beyond that setting and identify the Borel--Laplace representation $L_A(x)=\int_0^\infty e^{-xs}A(-s)\,\dd s$, whenever the integral exists.

For the Gaussian kernel $A(-s)=e^{-s^2/2}$, this representation gives the \emph{Mills--Hermite} expansion $M(x+t)\sim\sum(-1)^n\He_n(t)x^{-n-1}$, extending the classical scalar expansion.

We establish an explicit remainder estimate and derive finite-difference formulas that remain in the Hermite polynomial algebra.

As an application, we obtain a large-threshold expansion of the non-central Gaussian tail with explicit polynomial dependence on the mean.

If $A(-s)$ is itself the Laplace transform of a positive measure, then the coefficients form a Stieltjes moment sequence.

The associated Padé convergents give a systematic hierarchy of two-sided bounds in which successive lower and upper approximants incorporate the moments one at a time.