Converging Factors of a Class of Superfactorially Divergent Stieltjes Series

1. Introduction

Stieltjes series are fundamental tools in mathematical physics and continue to receive considerable attention. For instance, the, factorially divergent perturbation expansion of the energy eigenvalue of the $\mathcal{PT}$-symmetric Hamiltonian $\mathcal{H}\left(\lambda \right)=-{\displaystyle \frac{{\mathrm{d}}^2}{\mathrm{d}{x}^2}}+{\displaystyle \frac{1}{4}}{x}^2+\mathrm{i}\lambda x^3$ has been proved to be Stieltjes [1,2], as conjectured ten years earlier [3]. More recently, it was shown that even the character of the celebrated Bessel solution of Kepler’s equation [4] belongs to the Stieltjes family [5].

The present paper constitutes the follow up of a previous work on the study of convergence factors of Stieltjes series [6]. A well-established general convergence theory based on Padé approximants already exists for Stieltjes series. For example, if the moment sequence of a given Stieltjes series satisfies Carleman’s condition, then sequences of its diagonal or near-diagonal Padé approximants are guaranteed to converge to the Stieltjes function that generates the series itself (see for instance Baker and Graves-Morris [7]). This has consecrated Padé approximants as the main computational tool for the resummation of Stieltjes series. However, Padé approximants are subject to intrinsic limitation, especially when dealing with wildly divergent series (e.g., such that Calerman’s condition is not satisified). This has led to the development of new types of sequence transformations for summing divergent series. Among them, Levin-type transformations [8,9,10] proved to be particularly effective and powerful, in some cases outperforming Padé-based methods, such as Wynn’s epsilon algorithm [11]. Within the last four decades, an important literature has been produced, especially on nonlinear and nonregular sequence transformations, [9,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28].

About fifteen years ago, Ernst Joachim Weniger and I embarked on a challenging research project. We believed that a solid theoretical understanding of Levin-type sequence transformations – perhaps even a comprehensive convergence theory – could offer a valid alternative to Padé approximants for summing divergent Stieltjes series, particularly where Padé-based methods fail. Reference [28] demonstrated the remarkable computational effectiveness of Weniger’s $\delta$-transformation [9]. This specific Levin-type transformation not only successfully resummed the celebrated Euler series (a paradigm of factorial divergence) but did it with convergence rates greater than Padé. The proof in [28] relied on an inverse factorial expansion for Euler series converging factors, discovered a few years before [29]. Factorial series are often overlooked mathematical objects, largely unknown to non-specialists. Weniger’s merits include unearthing them while developing his $\delta$-transformation.

Ernst Joachim Weniger passed away on August 10th, 2022. Two years later, his contributions and legacy were celebrated in Ref. [6], where a constructive proof that converging factors of typical Stieltjes series can be expressed through inverse factorial series, was proposed. This proof was cast as an algorithm based on a first-order difference equation, which has been shown to be satisfied by the convergent factors of any Stieltjes series [6]. The present paper directly builds upon that tribute, providing a significant continuation of our work. Specifically, the algorithm proposed in [6] is here tested on a class of Stieltjes series with superfactorial moment growths. Weniger’s $\delta$-transformation has previously succeeded in resumming extremely divergent perturbation expansions, closely related to this class of Stieltjes series [30,31,32,33]. In some of these cases, Padé approximants often proved to be ineffective, as seen with the Rayleigh-Schrödinger perturbation series for the sextic anharmonic oscillator, or even failed in the more challenging octic case [34].

In Section 2, the main definitions and properties of Stieltjes series and Stieltjes functions are briefly reviewed, together with a resume of Ref. [6]. In Section 3, the class of superfactorially divergent series is presented, and the inverse factorial representation of the moment ratio sequence is analytically found. Section 4 illustrates some examples of application of our algorithm, to show its effectiveness and its ease of implementation. Finally, in Section 5 the converging factor representation problem is reformulated from a different perspective, by transforming the recurrence relation of [6] into a linear Cauchy’s problem, whose explicit solution is obtained using Faà di Bruno’s formula, together with the use of Bell’s polynomials [35] (Section 3.3).

For an improved readability, the most technical parts have been relegated to appendices.

2. Stieltjes Series Converging Factors, Factorial Series, and the Moment Problem

Consider a nondecreasing, real-valued function $\mu \left(t\right)$ defined for $t\in [0,\infty ]$, possessing infinitely many points of increase. This ensures that the associated measure, say $\mathrm{d}\mu$, is positive on $[0,\infty ]$. It will be assumed that all of its moments, defined as:

$${\displaystyle {\mu }_m={\int }_0^{\infty }{t}^m\mathrm{d}\mu ,m\ge 0,}$$

(1)

are finite. The formal power series

$${\displaystyle \sum _{m=0}^{\infty }{\displaystyle \frac{{(-)}^m}{z^{m+1}}}{\mu }_m,}$$

(2)

is called a Stieltjes series. Such series turns out to be asymptotic, in the sense of Poincaré, for $z\to \infty$, to the function $F\left(z\right)$ defined by

$${\displaystyle F\left(z\right)={\int }_0^{\infty }\frac{\mathrm{d}\mu }{z+t},|arg\left(z\right)|<\pi ,}$$

(3)

which is analytic in the complex plane cut along the negative real axis (i.e., $\mathbb{C}\setminus (-\infty ,0])$, and is called Stieltjes function.

The probably most known example of Stieltjes series is the Euler series [36], characterized by the moment sequence ${\{{\mu }_m=m!\}}_{m=0}^{\infty }$, and asymptotic to the Euler integral,

$${\displaystyle {\int }_0^{\infty }\frac{exp(-t)\mathrm{d}t}{z+t},|arg\left(z\right)|<\pi ,}$$

(4)

which has the form given in Eq. (3), with the measure $\mathrm{d}\mu =exp(-t)\mathrm{d}t$.

Given a sequence of moments ${\left\{{\mu }_m\right\}}_{m=0}^{\infty }$, is the corresponding Stieltjes function $F\left(z\right)$ uniquely determined? And, if so, is it possible to decode the asymptotic series into Eq. (2) to retrieve the correct value of F? The solution of such fundamental problem, which is known as the Stieltjes moment problem, depends on the growth rate of the moments. Carleman’s condition represents an important sufficient criterion to assess unicity to the moment problem. In particular, the Stieltjes moment problem is determinate if the series

$${\displaystyle \sum _{m=0}^{\infty }{\mu }_m^{-\frac{1}{2m}},}$$

(5)

diverges.

Any Stieltjes function $F\left(z\right)$ can be expressed as the sum of the nth-order partial sum of the associated asympotic series (2) and of a truncation error which has itself the form of a Stieltjes integral (see for example [9] (Theorem 13-1)). More precisely, we have

$$F\left(z\right)=\sum _{m=0}^n{\displaystyle \frac{{(-1)}^m}{z^{m+1}}}{\mu }_m+{\left(-{\displaystyle \frac{1}{z}}\right)}^{n+1}{\int }_0^{\infty }{\displaystyle \frac{t^{n+1}\mathrm{d}\mu }{z+t}},|arg\left(z\right)|<\pi ,$$

(6)

and the truncation error term can always be recast as follows:

$${\displaystyle {\displaystyle \frac{{(-)}^{n+1}}{z^{n+1}}}{\mu }_{n+1}{\phi }_{n+1}\left(z\right),}$$

(7)

where the function

$${\displaystyle {\phi }_m\left(z\right)=\frac{1}{{\mu }_m}{\int }_0^{\infty }{t}^m\frac{\mathrm{d}\mu }{t+z},m\in {\mathbb{N}}_0,|arg\left(z\right)|<\pi ,}$$

(8)

will be called the mth-order converging factor [37,38].1

From Eq. (7), it appears that, if reasonable estimates of the converging factor could be achieved without resorting to the numerical evaluation of the integral into Eq. (8), then the numerical evaluation of $F\left(z\right)$ via its asymptotics would be in principle possible. The search of techniques aimed at estimating convergence factors plays a role of pivotal importance to decode divergent asymptotic series. In [6] it has been shown that, given a Stieltjes series, the converging factor in Eq. (8) can always be represented as an inverse factorial series. Our proof was ultimately based on (i) the fact that ${\phi }_m\left(z\right)$ must satisfy the following first-order difference equation [6]:

$${\displaystyle {\phi }_{m+1}=\frac{{\mu }_m}{{\mu }_{m+1}}\left(1-z{\phi }_m\right),m\ge 0,}$$

(9)

and (ii) that inverse factorial series constitute natural tools for solving difference equations, similarly as inverse power series are used to solve differential equations. For reader’s convenience, the basic definitions and properties of factorial series will now be briefly recalled, although extensive reviews can be found, for instance, in [6,28]. More interested and motivated readers are encouraged to go through the Weniger paper [25], where a hystorical account of his re-discovery of factorial series can be found, together with a list of their most important computational features. For the scopes of the present paper, it is sufficient to limit ourselves to the following key points:

(i)

Let $f\left(x\right)$ be a complex function of a complex variable x. A factorial series representation of $f\left(x\right)$ is an expansion of the following type:

$$f\left(x\right)={\displaystyle \frac{a_0}{x}}+{\displaystyle \frac{1!}{x(x+1)}}{a}_1+{\displaystyle \frac{2!}{x(x+1)(x+2)}}{a}_2+\dots =\sum _{k=0}^{\infty }{\displaystyle \frac{k!}{{\left(x\right)}_{k+1}}}{a}_k,$$

(10)

where the symbol ${\left(x\right)}_{k+1}=x(x+1)\dots (x+k)$ denotes the Pochhammer symbol. In the following, it will be assumed $x>0$.

(ii)

Compared to inverse power series, factorial series often possess superior convergence properties. For example, consider the divergent sequence ${\{k!\}}_{k=0}^{\infty }$ and construct the following two infinite series:

$${\displaystyle \frac{1}{x}}-{\displaystyle \frac{1}{x^2}}+{\displaystyle \frac{2}{x^3}}-{\displaystyle \frac{6}{x^4}}+\dots =\sum _{k=0}^{\infty }{\displaystyle \frac{{(-1)}^{k}k!}{x^{k+1}}},$$

(11)

and

$${\displaystyle \frac{1}{x}}-{\displaystyle \frac{1}{{\left(x\right)}_2}}+{\displaystyle \frac{2}{{\left(x\right)}_3}}-{\displaystyle \frac{6}{{\left(x\right)}_4}}+\dots =\sum _{k=0}^{\infty }{\displaystyle \frac{{(-1)}^{k}k!}{{\left(x\right)}_{k+1}}}.$$

(12)

The power series diverges for all $\left|x\right|\in {\mathbb{R}}_{+}$, whereas the factorial series converges for all $x>1$. In other words, it may happen that a given function possessing a representation in terms of a divergent asymptotic series, also possesses a representation as a convergent factorial series.

More precisely, based on the following asymptotics:

$${\displaystyle \frac{k!}{{\left(x\right)}_{k+1}}}=\mathrm{O}\left(x^{-k}\right),k\to \infty ,$$

(13)

it is possible to prove that the factorial series into Eq. (10) converges if and only if the associated Dirichlet series, i.e.,

$$\sum _{n=1}^{\infty }{\displaystyle \frac{a_n}{n^x}},$$

(14)

converges too (see for example [39] (p. 262) or [40] (p. 167)).

(iii)

Factorial series admit useful integral representations. Starting from

$${\displaystyle \frac{k!}{{\left(x\right)}_{k+1}}}={\int }_0^1{t}^{x-1}{(1-t)}^k\mathrm{d}t,\mathrm{Re}\left(x\right)>0k\in {\mathbb{N}}_0,$$

(15)

and on inserting into Eq. (10), after interchanging integration and summation, the following integral representation is found [41] (Sec. I on p. 244):

$$f\left(x\right)={\int }_0^1{t}^{x-1}\Phi \left(t\right)\mathrm{d}t,\mathrm{Re}\left(x\right)>0,$$

(16)

where the function $\Phi \left(t\right)$ turns out to be

$$\Phi \left(t\right)=\sum _{k=0}^{\infty }{a}_k{(1-t)}^k.$$

(17)

Note that the function $f\left(x\right)$ can be thought of as a truncated Mellin transform of the function $\Phi \left(t\right)$, in terms of which the expanding coefficients into Eq. (10) take on the following form:

$$a_k={\displaystyle \frac{{(-)}^k}{k!}}{\Phi }^{\left(k\right)}\left(1\right),k\ge 1,$$

(18)

where the symbol ${\Phi }^{\left(k\right)}$ denotes the kth-order derivative of the function $\Phi \left(t\right)$. In other words, $\Phi \left(t\right)$ is nothing but the generating function of the sequence ${\{k!a_k\}}_{k=0}^{\infty }$.

The main results of the algorithm derived in Ref. [6] will now be briefly recalled. The idea consisted in searching the solution of the difference equation (9) in the form of the following inverse factorial series:

$${\displaystyle {\phi }_m\left(z\right)=\sum _{k=0}^{\infty }{\displaystyle \frac{k!}{{(m+\beta )}_{k+1}}}{c}_k,}$$

(19)

where $\beta$ is a real nonnegative parameter, and ${\left\{c_k\right\}}_{k=0}^{\infty }$ denotes a sequence which is independent of m. For the sake of simplicity, it will be assumed $\beta =0$ henceforth. What has been found in Ref. [6] is that the expanding coefficients $c_k$’s can be obtained in a very easy way if the moment ratio sequence ${\{{\mu }_m/{\mu }_{m+1}\}}_{m=0}^{\infty }$ admits itself an inverse factorial expansion, i.e.,

$${\displaystyle {\displaystyle \frac{{\mu }_m}{{\mu }_{m+1}}}=\sum _{k=0}^{\infty }\frac{k!}{{\left(m\right)}_{k+1}}{b}_k.}$$

(20)

In particular, the following recurrence relation holds [6]:

$$\left\{\begin{array}{c}{\displaystyle c_0=b_0,}\hfill \\ \\ c_k=c_{k-1}+b_k-z{d}_k,k\ge 1,\hfill \end{array}\right.$$

(21)

where the sequence ${\left\{d_k\right\}}_{k=1}^{\infty }$ is given by

$$\begin{array}{c}{\displaystyle d_k={\displaystyle \frac{1}{k!}}\sum _{\nu =0}^{k-1}\sum _{\mu =0}^{\nu }\nu !{(\mu +1)}_{k-1-\nu }{b}_{k-1-\nu }{c}_{\nu -\mu },}\hfill \end{array}$$

(22)

or, equivalently, by

$$\begin{array}{c}{\displaystyle d_k={\displaystyle \frac{1}{k!}}\sum _{\nu =0}^{k-1}\sum _{\lambda =0}^{\nu }\nu !{(\nu -\lambda +1)}_{k-1-\nu }{b}_{k-1-\nu }{c}_{\lambda }.}\hfill \end{array}$$

(23)

In Ref. [6], Eqs. (20) - (23) have been applied to re-derive and further generalize (to any $\beta >0$) the result found in [29] as far as the converging factor of the Euler series is concerned. In the next section, they will be employed to find the inverse factorial series representation of the converging factor of an important class of superfactorial divergent Stieltjes series, which can be thought as an important suitable generalization of the Euler series.

3. A Class of Superfactorially Divergent Stieltjes Series

The class of Stieltjes series under investigation are characterized by the following moment sequence:

$${\mu }_m=\Gamma (\nu m+1+q),m\ge 0,$$

(24)

where $\nu \in {\mathbb{Z}}_{+}$ and $q\in (-1,1)$. Euler series corresponds to the choice $(\nu ,q)=(1,0)$. Asymptotic series of Eq. (24) play a role of pivotal importance in theoretical as well as mathematical physics. They are strictly related to the theory of quantum anharmonic oscillators. In two seminal papers, Bender and Wu [42,43] shown that the Rayleigh-Schrödinger perturbative expansion of the energy levels of the Hamiltonian

$$\widehat{\mathcal{H}}\left(\epsilon \right)=-{\displaystyle \frac{{\mathrm{d}}^2}{\mathrm{d}{x}^2}}+x^2+\epsilon x^{2(\nu +1)},$$

(25)

with $\nu =1,2,3,\dots$, is asymptotically dominated by the Stieltjes series in Eq. (24), with $q=-1/2$. The summability of such asymptotic series has also been addressed in 1993 by Weniger et al. [33] through a series of important numerical experiments, aimed at comparing the retrieving performances of Levin’s type nonlinear transformations with those of Padé approximants, as far as the computation of the energy levels of quartic, sextic, and octic oscillators was concerned. In particular, they also considered the numerical resummation of Stieltjes series of the type into Eq. (24) for $q=-1/2$ and for some integer values of $\nu$. In particular, it was found that both Levin’s and Weniger’s transformations were able to decode correctly the series when $\nu =1,2,3$, whereas Padé approximants failed to achieve the task already for $\nu =3$, i.e., when the Calerman condition in Eq. (5) is not satisfied.

In the present section, the estimation of the converging factor of the Stieltjes series of Eq. (24) will be obtained on applying the algorithm of [6]. Our analysis will be carried out for any values of $\nu$ and q, while it will be set $\beta =0$ in all subsequent calculations, as said above. In order for the algorithm in Section 2 to be applied, it is mandatory to find the inverse factorial expantion of the moment ratio sequence for the class of Stieltjes series defined by Eq. (24), which can be recast as follows:

$$\begin{array}{c}{\displaystyle {\displaystyle \frac{{\mu }_m}{{\mu }_{m+1}}}={\displaystyle \frac{1}{(\nu m+q+1)(\nu m+q+2)\dots (\nu m+q+\nu )}}={\displaystyle \frac{{\nu }^{-\nu }}{{\displaystyle \prod _{\ell =1}^{\nu }\left(m+{\displaystyle \frac{q+\ell }{\nu }}\right)}}},\nu \in {\mathbb{Z}}_{+}.}\hfill \end{array}$$

(26)

On expanding the right side of Eq. (26) as a sum of partial fraction of the form

$${\displaystyle {\displaystyle \frac{{\mu }_m}{{\mu }_{m+1}}}=\sum _{\ell =1}^{\nu }{\displaystyle \frac{A_{\ell }}{m+{\displaystyle \frac{q+\ell }{\nu }}}},}$$

(27)

it is not difficult to prove that (see Appendix A)

$${\displaystyle A_{\ell }={\displaystyle \frac{{(-1)}^{\ell +1}}{\nu (\nu -\ell )!(\ell -1)!}}.\ell =1,2,\dots ,\nu ,}$$

(28)

The subsequent step consists in using Waring’s formula [41] (Eq. (3) on p. 77), namely

$${\displaystyle \frac{1}{z-w}}=\sum _{k=0}^{\infty }{\displaystyle \frac{{\left(w\right)}_k}{{\left(z\right)}_{k+1}}},\mathrm{Re}(z-w)>0,$$

(29)

which, together with Eqs. (27) and (28), after some algebra leads to Eq. (20) with the following expanding coefficient sequence (see Appendix B):

$$\begin{array}{c}{\displaystyle b_k={\displaystyle \frac{1}{\nu !k!}}\sum _{\ell =0}^{\nu -1}{(-1)}^{\ell }\left(\genfrac{}{}{0pt}{}{\nu -1}{\ell }\right){\left(-{\displaystyle \frac{\ell +q+1}{\nu }}\right)}_k,k\ge 0.}\hfill \end{array}$$

(30)

It should be noted how the mathematical structure of the right side of Eq. (30) is nothing but the $(\nu -1)$th-order forward difference of the kth-degree polynomial ${\mathcal{P}}_k\left(x\right)={\left(-{\displaystyle \frac{x}{\nu }}\right)}_k$, i.e.,

$$\begin{array}{c}{\displaystyle b_k={\displaystyle \frac{{(-1)}^{\nu -1}}{\nu !k!}}{\Delta }^{\nu -1}{\left\{{\mathcal{P}}_k\left(x\right)\right\}}_{x=q+1},k\ge 0,}\hfill \end{array}$$

(31)

where the forward difference operator is defined as $\Delta {\mathcal{P}}_k\left(x\right)={\mathcal{P}}_k(x+1)-{\mathcal{P}}_k\left(x\right)$. In particular, since ${\mathcal{P}}_k\left(x\right)$ is a kth-order polynomial, it follows at once that:

$$\begin{array}{c}{\displaystyle b_k=0,0\le k<\nu -1.}\hfill \end{array}$$

(32)

The above result can be easily extended to deal with the most general inverse factorial expansion of the moment ratio. In fact, it is not difficult to show that, in order to deal with $\beta \ne 0$, it is sufficient to recast Eq. (26) as follows:

$${\displaystyle {\displaystyle \frac{{\mu }_m}{{\mu }_{m+1}}}={\displaystyle \frac{{\nu }^{-\nu }}{{\displaystyle \prod _{\ell =1}^{\nu }\left[\left(m+\beta \right)+\left({\displaystyle \frac{q+\ell }{\nu }}-\beta \right)\right]}}},}$$

(33)

and then to proceed similarly as we did for the case $\beta =0$. The final result, which is given without proof, is the following:

$$\begin{array}{c}{\displaystyle {\displaystyle \frac{{\mu }_m}{{\mu }_{m+1}}}=\sum _{k=0}^{\infty }{\displaystyle \frac{k!}{{(m+\beta )}_{k+1}}}{b}_k,\beta >0,}\hfill \end{array}$$

(34)

with the expanding coefficient sequence ${\left\{b_k\right\}}_{k=0}^{\infty }$ being now given by

$$\begin{array}{c}{\displaystyle b_k={\displaystyle \frac{1}{\nu !k!}}\sum _{\ell =0}^{\nu -1}{(-1)}^{\ell }\left(\genfrac{}{}{0pt}{}{\nu -1}{\ell }\right){\left(\beta -{\displaystyle \frac{\ell +q+1}{\nu }}\right)}_k,k\ge 0.}\hfill \end{array}$$

(35)

Equations (34) and (35) represent the most general inverse factorial expansion of the moment ratio sequence into Eq. (26). However, as previously said, in the rest of the paper we shall limit ourselves to deal with the case $\beta =0$.

4. Numerical Results

In the present section, a few examples concerning the numerical estimation of the converging factors of some of the divergent superfactorial Stieltjes series analyzed in the previous section, will be presented. Our attention will be directed to a series of important numerical experiments carried out in 1993 by Weniger et al. [33]. In particular, our computational targets are two Stieltjes integrals, namely2

$$\begin{array}{c}{\displaystyle {\mathcal{J}}_3={\int }_0^{\infty }{\displaystyle \frac{t^{-1/2}exp(-t)}{1+{\displaystyle \frac{64}{45{\pi }^2}}{t}^2}}\mathrm{d}t,}\hfill \end{array}$$

(36)

and

$$\begin{array}{c}{\displaystyle {\mathcal{J}}_4={\int }_0^{\infty }{\displaystyle \frac{t^{-1/2}exp(-t)}{1+{\displaystyle \frac{225{\left[\Gamma (2/3)\right]}^9}{112{\pi }^6}}{t}^3}}\mathrm{d}t.}\hfill \end{array}$$

(37)

It is not difficult to show that both integrals into Eqs. (36) and (37) can be recast in the form of Eq. (3). For the first integral, it is sufficiente to change the integration variable from t to $\sqrt{t}$ and to let $z={\displaystyle \frac{45{\pi }^2}{64}}$, so to have

$$\begin{array}{c}{\displaystyle {\mathcal{J}}_3=z{\int }_0^{\infty }{\displaystyle \frac{\mathrm{d}{\mu }_3}{z+t}},}\hfill \end{array}$$

(38)

where the measure $\mathrm{d}{\mu }_3$ is given by

$$\begin{array}{c}{\displaystyle \mathrm{d}{\mu }_3={\displaystyle \frac{1}{2}}{t}^{-3/4}exp(-\sqrt{t})\mathrm{d}t,}\hfill \end{array}$$

(39)

and gives at once

$$\begin{array}{c}{\displaystyle {\mu }_m={\int }_0^{\infty }{t}^m\mathrm{d}{\mu }_3=\Gamma \left(2m+{\displaystyle \frac{1}{2}}\right),m\in {\mathbb{N}}_0,}\hfill \end{array}$$

(40)

corresponding to the pair $(\nu ,q)=\left(2,-{\displaystyle \frac{1}{2}}\right)$ in the model given into Eq. (24). More importantly, the mth-order converging factor ${\phi }_m\left(z\right)$ can be analytically evaluated for any $m>0$ in terms of hypergeometric functions,3

$$\begin{array}{c}{\displaystyle {\phi }_m\left(z\right)={\displaystyle \frac{\Gamma \left(2m-{\displaystyle \frac{3}{2}}\right)1F_2\left(1;{\displaystyle \frac{5}{4}}-m,{\displaystyle \frac{7}{4}}-m;-{\displaystyle \frac{z}{4}}\right)+\pi z^{m-\frac{3}{4}}cos\left[\pi \left(m+{\displaystyle \frac{1}{4}}\right)+\sqrt{z}\right]}{\Gamma \left(2m+{\displaystyle \frac{1}{2}}\right)}}.}\hfill \end{array}$$

(41)

Similar results hold also for the Stieltjes integral ${\mathcal{J}}_4$, which can be recast as

$$\begin{array}{c}{\displaystyle {\mathcal{J}}_4=z{\int }_0^{\infty }{\displaystyle \frac{\mathrm{d}{\mu }_4}{z+t}},}\hfill \end{array}$$

(42)

where now $z={\displaystyle \frac{112{\pi }^6}{225{\left[\Gamma (2/3)\right]}^9}}$ and

$$\begin{array}{c}{\displaystyle \mathrm{d}{\mu }_4={\displaystyle \frac{1}{3}}{t}^{-5/6}exp\left(-t^{1/3}\right)\mathrm{d}t,}\hfill \end{array}$$

(43)

so that

$$\begin{array}{c}{\displaystyle {\mu }_m={\int }_0^{\infty }{t}^m\mathrm{d}{\mu }_4=\Gamma \left(3m+{\displaystyle \frac{1}{2}}\right),m\in {\mathbb{N}}_0,}\hfill \end{array}$$

(44)

corresponding to $(\nu ,q)=\left(3,-{\displaystyle \frac{1}{2}}\right)$ in the model given into Eq. (24). Similarly as happened for ${\mathcal{J}}_3$, also the converging factor of ${\mathcal{J}}_4$ can be found exactly, although the resulting expression is quite complicated and annoying.

Our algorithm will now be tested starting from ${\mathcal{J}}_3$, by focusing the attention on the relative error obtained when the mth-order converging factor ${\phi }_m\left(z\right)$ into Eq. (41) is going to be evaluated through the following truncated inverse factorial series:

$${\displaystyle {\phi }_m^{\left(N\right)}\left(z\right)=\sum _{k=0}^N{\displaystyle \frac{k!}{{\left(m\right)}_{k+1}}}{c}_k.}$$

(45)

Figure 1 shows the behaviour of the relative error, defined as

$${\displaystyle {ϵ}_m^{\left(N\right)}\left(z\right)=\left|1-{\displaystyle \frac{{\phi }_m^{\left(N\right)}\left(z\right)}{{\phi }_m\left(z\right)}}\right|,}$$

(46)

related to the converging factor defined by Eq. (41), as far as the integral ${\mathcal{J}}_3$ is concerned, for $m=5$ (dashed curve), $m=10$ (dotted curve), and $m=15$ (solid curve).

It can be appreciated how the inverse factorial series representation of ${\phi }_m$ appears to converge as $N\to \infty$, with a convergence rate becoming larger and larger on increasing the order m, as it might be expected from Eq. (3) and from the general considerations about inverse factorial series recalled at the beginning of Section 2.

Even more interesting are the numerical results obtained as far as the estimate of the converging factor of the Stieltjes integral ${\mathcal{J}}_4$ is concerned, which are shown in Figure 2.

In particular, it should be noted how, for both cases, the performances of our algorithm (e.g., for what concerns the order of magnitude of the relative error as a function of the factorial series truncation order) are comparable, despite the fact that choosing $\nu =2$ or $\nu =3$ determines or not the validity of Calerman’s condition, respectively. As it was put into evidence in Ref. [33], such a circumstance greatly influenced the numerical performances of Padé approximants in retrieving the correct values of ${\mathcal{J}}_3$ as well as ${\mathcal{J}}_4$, while did not affect the retrieving action of Levin’s, as well as Weniger’s transformations [33]. These preliminary results corroborate our feeling about the robustness of Levin-type transformations in resumming superfactorially divergent Stieltjes series. The final part of the present paper is aimed at exploring further properties of the key role played by of the factorial expansion of converging factors, and to give some general guidelines for future estimations of the related convergence rates.

5. On Integral Representation of the Converging Factors for $Q=0$

In the present section, some of the results found in Section 3 will now be re-derived in a completely different way, on employing the integral representation of factorial series of Eq. (16). We are convinced that what is going to be presented could reveal of a certain importance in future, in order for the convergence features of Levin-type nonlinear transformations in resumming superfactorially divergent Stieltjes series to be explored, similarly to what has been found for the Euler series [28]. There, the starting point of our convergence analysis was just the truncated Mellin transform representation of the Euler series converging factor, which was obtained from the results published in [29]. The same results will now be rederived within a few steps, on using the recurrence relation in Eq. (9) for the converging factor ${\phi }_m\left(z\right)$, together with Eq. (24) written for $(\nu ,q)=(1,0)$, i.e.,

$${\displaystyle (m+1){\phi }_{m+1}=1-z{\phi }_m,m\ge 0.}$$

(47)

so that, on taking Eq. (16) into account, the following integral equation for the generating function $\Phi \left(t\right)$ is obtained:

$${\displaystyle {\int }_0^1(m+1)t^m\Phi \left(t\right)\mathrm{d}t+z{\int }_0^1{t}^{m-1}\Phi \left(t\right)\mathrm{d}t=1.}$$

(48)

On integrating by parts the first integral, we obtain

$${\displaystyle {\int }_0^1(m+1)t^m\Phi \left(t\right)\mathrm{d}t=\Phi \left(1\right)-{\int }_0^1{t}^{m+1}{\Phi }^{\prime }\left(t\right)\mathrm{d}t,}$$

(49)

with the superscript denoting derivation with respect to t. Finally, on substituting from Eq. (49) into Eq. (48), our integral equation is transformed into the following first-order Cauchy problem for $\Phi \left(t\right)$:

$$\left\{\begin{array}{c}{\displaystyle t^2{\Phi }^{\prime }\left(t\right)-z\Phi \left(t\right)=0,}\hfill \\ \\ \Phi \left(1\right)=1,\hfill \end{array}\right.$$

(50)

whose solution is

$${\displaystyle \Phi \left(t\right)=exp\left[-z(t^{-1}-1)\right],}$$

(51)

and that leads to [28] (Eq. (5.17)).4

It is then natural to ask whether similar results could be found also for the pairs $(\nu >1,q=0)$, following the same strategy. For the sake of simplicity, only the cases $\nu =2$ and $\nu =3$ will now be detailed. A conjecture, which we thought to be valid for $\nu >3$, will be proposed as an open problem at the end of the paper.

Consider now the case $\nu =2$. From Eq. (24) we have ${\mu }_m=\Gamma (2m+1)$ which, once substituted into Eq. (9), gives

$$4(m+1)\left(m+{\displaystyle \frac{1}{2}}\right){\phi }_{m+1}+z{\phi }_m=1,$$

(52)

that can also be recast as follows:

$$\begin{array}{c}{\displaystyle 4{(m+1)}^2{\phi }_{m+1}-2(m+1){\phi }_{m+1}+z{\phi }_m=1.}\hfill \end{array}$$

(53)

Similarly as done for the Euler series, to solve the difference equation the integral representation in Eq. (16) will now be used, together with the two relationships

$$\begin{array}{c}{\displaystyle m{\phi }_m=\Phi \left(1\right)-{\int }_0^1{t}^m{\Phi }^{\prime }\mathrm{d}t,}\hfill \\ \\ {\displaystyle m^2{\phi }_m=m\Phi \left(1\right)-{\Phi }^{\prime }\left(1\right)+{\int }_0^1{t}^m{\left(t{\Phi }^{\prime }\right)}^{\prime }\mathrm{d}t,}\hfill \\ \end{array}$$

(54)

that can be proved again via partial integration. Substitution from Eq. (54) into Eq. (53) gives the following integral equation for the function $F\left(t\right)$:

$$\begin{array}{c}{\displaystyle (4m+2)\Phi \left(1\right)-4{\Phi }^{\prime }\left(1\right)+{\int }_0^1{t}^{m-1}\left[4t^2{\left(t{\Phi }^{\prime }\right)}^{\prime }+2t^2{\Phi }^{\prime }+z\Phi \right]\mathrm{d}t=1,}\hfill \end{array}$$

(55)

which leads to the Cauchy problem

$$\left\{\begin{array}{c}{\displaystyle t^2{\left(t{\Phi }^{\prime }\right)}^{\prime }+\frac{1}{2}{t}^2{\Phi }^{\prime }+\frac{z}{4}\Phi =0,}\hfill \\ \\ {\displaystyle \Phi \left(1\right)=0,{\Phi }^{\prime }\left(1\right)=-\frac{1}{4}.}\hfill \end{array}\right.$$

(56)

On letting $\Phi \left(t\right)=G(t^{-1/2}-1)$, straigthforward algebra transforms Eq. (56) into

$$\left\{\begin{array}{c}{\displaystyle G^{\prime \prime }+zG=0,}\hfill \\ \\ {\displaystyle G\left(0\right)=0,G^{\prime }\left(0\right)=\frac{1}{2},}\hfill \end{array}\right.$$

(57)

whose solution is $G\left(\xi \right)=sin\left(\sqrt{z}\xi \right)/2\sqrt{z}$. Accordingly, we have

$$\begin{array}{c}{\displaystyle \Phi \left(t\right)=\frac{1}{2\sqrt{z}}sin\left[\sqrt{z}\left(\frac{1}{\sqrt{t}}-1\right)\right].}\hfill \end{array}$$

(58)

It is not difficult to check how, after substituting from Eq. (58) into Eq. (18), the whole sequence so generated coincides with the sequence ${\left\{c_k\right\}}_{m=0}^{\infty }$ obtained through the recursive algorithm into Eqs. (21) and (23), for any z.

The case $\nu =3$ is particularly intriguing, due to the fact that when ${\mu }_m=\Gamma (3m+1)$ Calerman’s condition is not satisfied and Padé approximants are no longer able to decode the associated Stieltjes series. As far as Eq. (9) is concerned, we have

$$27(m+1)\left(m+{\displaystyle \frac{2}{3}}\right)\left(m+{\displaystyle \frac{1}{3}}\right){\phi }_{m+1}+z{\phi }_m=1,$$

(59)

which can be recast as follows:

$$\begin{array}{c}{\displaystyle 27{(m+1)}^3{\phi }_{m+1}-27{(m+1)}^2{\phi }_{m+1}+6(m+1){\phi }_{m+1}+z{\phi }_m=1.}\hfill \end{array}$$

(60)

On recalling Eqs. (54) and on taking into account that

$$\begin{array}{c}{\displaystyle m^3{\phi }_m=m^2\Phi \left(1\right)-(m-1){\Phi }^{\prime }\left(1\right)+{\Phi }^{\prime \prime }\left(1\right)-{\int }_0^1{t}^m{\left[t{\left(t{\Phi }^{\prime }\right)}^{\prime }\right]}^{\prime }\mathrm{d}t,}\hfill \end{array}$$

(61)

after some algebra it is possible to show that the function $\Phi \left(t\right)$ must satisfy the following Cauchy problem:

$$\left\{\begin{array}{c}{\displaystyle t^2{\left[t{\left(t{\Phi }^{\prime }\right)}^{\prime }\right]}^{\prime }+t^2{\left(t{\Phi }^{\prime }\right)}^{\prime }+\frac{2}{9}{t}^2{\Phi }^{\prime }-\frac{z}{27}\Phi =0,}\hfill \\ \\ {\displaystyle \Phi \left(1\right)=0,{\Phi }^{\prime }\left(1\right)=0,{\Phi }^{\prime \prime }\left(1\right)={\displaystyle \frac{1}{27}},}\hfill \end{array}\right.$$

(62)

that, on letting $\Phi \left(t\right)=G(t^{-1/3}-1)$, transforms into

$$\left\{\begin{array}{c}{\displaystyle G^{\prime \prime \prime }+zG=0,}\hfill \\ \\ {\displaystyle G\left(0\right)=0,G^{\prime }\left(0\right)=0,G^{\prime \prime }\left(0\right)=\frac{1}{3},}\hfill \end{array}\right.$$

(63)

whose explicit solution is given by

$$\begin{array}{c}{\displaystyle G\left(\xi \right)=\frac{z^{-2/3}}{9}exp(-z^{1/3}\xi )-\frac{2}{9}{z}^{-2/3}exp\left(\frac{1}{2}{z}^{1/3}\xi \right)cos\left(\frac{\sqrt{3}}{2}{z}^{1/3}\xi +\frac{\pi }{3}\right).}\hfill \end{array}$$

(64)

The above results would suggest that a factorial expansion of the converging factor could be costructed also for $\nu >3$. In particular, we conjecture that the generating function $\Phi \left(t\right)$ can alwasy be recast as follows:

$$\begin{array}{c}{\displaystyle \Phi \left(t\right)=G\left(t^{-1/\nu }-1\right),}\hfill \end{array}$$

(65)

with the function $G\left(\xi \right)$ being the solution of the $\nu$th-order Cauchy problem

$$\left\{\begin{array}{c}{\displaystyle \frac{{\mathrm{d}}^{\nu }G}{\mathrm{d}{\xi }^{\nu }}+zG=0,}\hfill \\ \\ {\displaystyle G\left(0\right)=G^{\prime }\left(0\right)=\dots =G^{(\nu -2)}\left(0\right)=0,}\hfill \\ \\ {\displaystyle G^{(\nu -1)}\left(0\right)=\frac{1}{\nu }.}\hfill \end{array}\right.$$

(66)

In particular, the extraction of the sequence $\left\{c_k\right\}$ could be done directly through Eq. (18), without explicitly solving Eq. (66). To this end, it is sufficient to note that, thanks to Eq. (65), the sequence ${\left\{{\Phi }^{\left(k\right)}\left(1\right)\right\}}_{k=0}^{\infty }$ can be generated through

$$\left\{\begin{array}{c}{\displaystyle \Phi \left(1\right)=G\left(0\right),}\hfill \\ \\ {\displaystyle {\Phi }^{\prime }\left(1\right)=-\frac{1}{\nu }{G}^{\prime }\left(0\right),}\hfill \\ \\ {\displaystyle {\Phi }^{\prime \prime }\left(1\right)=\frac{1}{{\nu }^2}{G}^{\prime \prime }\left(0\right)+\frac{1+\nu }{{\nu }^2}{G}^{\prime }\left(0\right),}\hfill \\ \\ {\displaystyle {\Phi }^{\prime \prime \prime }\left(1\right)=-\frac{1}{{\nu }^3}{G}^{\prime \prime \prime }\left(0\right)-3\frac{1+\nu }{{\nu }^3}{G}^{\prime \prime }\left(0\right)-\frac{(1+\nu )(1+2\nu )}{{\nu }^3}{G}^{\prime }\left(0\right),}\hfill \\ \\ \dots \hfill \end{array}\right.$$

(67)

and so on, where the sequence and ${\left\{G^{\left(k\right)}\left(0\right)\right\}}_{k=0}^{\infty }$ should be achievable via recurrenced as follows:

$$\begin{array}{c}{\displaystyle G^{\left(k\right)}\left(0\right)=-z{G}^{(k-\nu )}\left(0\right),k\ge \nu ,}\hfill \end{array}$$

(68)

starting from the initial values given in Eq. (66). More importantly, Eq. (67) can be given in closed-form terms on using Faà di Bruno’s formula, as shown in Appendix C, i.e.,

$$\begin{array}{c}{\displaystyle {\Phi }^{\left(k\right)}\left(1\right)=\sum _{\ell =0}^k{G}^{(\ell )}\left(0\right)B_{k,\ell }\left[-{\left({\displaystyle \frac{1}{\nu }}\right)}_1,{\left({\displaystyle \frac{1}{\nu }}\right)}_2,\dots {(-1)}^{k-\ell +1}{\left({\displaystyle \frac{1}{\nu }}\right)}_{k-\ell +1}\right],}\hfill \end{array}$$

(69)

where the symbol $B_{k,\ell }(\dots )$ denotes the partial Bell polynomial [35] (Section 3.3).

6. Conclusions

In the present paper, we explored the converging factors for a class of Stieltjes series exhibiting superfactorial divergence using the algorithm developed in [6]. We successfully demonstrated the numerical effectiveness of the algorithm and its ease of implementation. Our algorithm proved to be robust in estimating converging factors even for series whose moments do not satisfy Carleman’s condition, where Padé approximants largely fail.

We also introduced novel theoretical perspectives on the converging factor representation problem. By reducing the recurrence relation satisfied by the converging factor to a linear Cauchy problem, we provided explicit solutions using the Faà di Bruno formula and Bell’s polynomials. This could help to develop new strategies for the theoretical exploration of Levin-type sequence transformations convergence features. In particular, the integral representations of the converging factor obtained for the class of Stieltjes series investigated here are expected to be instrumental in estimating the convergence rates of Weniger’s $\delta$-transformation for these class of divergent series, similar to what has been done for the Euler series [28].