Next Article in Journal
Sequential Ignorability and Dismissible Treatment Components to Identify Mediation Effects
Previous Article in Journal
Invariants for Second Type Almost Geodesic Mappings of Symmetric Affine Connection Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Factorial Series Representation of Stieltjes Series Converging Factors

Departimento di Ingegneria Civile, Informatica e delle Tecnologie Aeronautiche, Università degli Studi “Roma Tre” 1, 00146 Rome, Italy
Mathematics 2024, 12(15), 2330; https://doi.org/10.3390/math12152330
Submission received: 1 June 2024 / Revised: 19 July 2024 / Accepted: 21 July 2024 / Published: 25 July 2024
(This article belongs to the Section Mathematical Physics)

Abstract

:
The practical usefulness of Levin-type nonlinear sequence transformations as numerical tools for the summation of divergent series or for the convergence acceleration of slowly converging series is nowadays beyond dispute. The Weniger transformation, in particular, is able to accomplish spectacular results when used to overcome resummation problems, often outperforming better-known resummation techniques, like, for instance, Padé approximants. However, our theoretical understanding of Levin-type transformations is still far from being satisfactory and is particularly bad as far as the decoding of factorially divergent series is concerned. The Stieltjes series represent a class of power series of fundamental interest in mathematical physics. In the present paper, it is shown how the converging factor of any order of typical Stieltjes series can be expressed as an inverse factorial series, whose terms are analytically retrieved through a simple recursive algorithm. A few examples of applications are presented, in order to show the effectiveness and implementation ease of the algorithm itself. We believe that further investigations of the asymptotic forms of the remainder terms, encoded within the converging factors, could eventually lead toward a more general theory of the asymptotic behavior of the Weniger transformation when it is applied to Stieltjes series in high transformation orders. It is a rather ambitious project, which should be worthy of being pursued in the future.

1. Prelude

The present paper is a fitting tribute to the memory of Ernst Joachim Weniger, mentor and friend, who passed away on 10 August 2022, sadly too soon. I had the good fortune to meet Ernst Joachim in 2009 at a conference organized at the Centre International de Rencontres Mathématiques in Marseille, France. Since then, I have had the pleasure of working with him for over a decade. In 2015, we published our first (and sole) joint paper, in which a convergence theory for the resummation of the Euler series via Levin-type nonlinear sequence transformations, was presented. In particular, we proved that his transformation, named the Weniger transformation, was not only able to resum the Euler series but also that, in accomplishing such a task, it turned out to be “exponentially faster” than Padé approximants. The Euler series is an example of the Stieltjes series, for which a well-established general convergence theory based on Padé approximants is well known to be available. After our 2015 joint paper, Ernst Joachim and I thought that trying to conceive a convergence theory also for Levin-type transformations, when applied to Stieltjes series, could have been an ambitious scientific project, but one worthy of being pursued.
What is contained in the present work should have been my second joint work with Ernst Joachim. In writing it, I tried to respect to the best of my ability all that Ernst Joachim had and would suggest to me, also thanks to the help of still unpublished material and several discussions that took place during his visits to my former Engineering Department at the University “Roma Tre.”

2. Introduction

Power series are among the most important mathematical objects, not only in classical analysis, but also in the mathematical treatment of several scientific and engineering problems. Most power series have a finite radius of convergence. If the argument of a power series is close to the expansion point, convergence is usually (very) good, but once the argument approaches the boundary of the circle of convergence, convergence can become arbitrarily bad. Outside the circle of convergence, power series diverge [1].
Especially in theoretical physics, wildly divergent power series occur abundantly which cannot be replaced easily by alternative expressions. The Stirling series for the logarithm of the gamma function ([2] Equation (5.11.1)), as well as the Poincaré seminal work on asymptotic series [3], are classic examples of useful divergent series. In quantum physics, divergent series have become indispensable. Dyson [4] first argued that perturbation expansions in quantum electrodynamics must diverge factorially. Bender and Wu [5,6,7] showed that factorially divergent perturbation expansions occur also in nonrelativistic quantum mechanics. Later, it was found that factorially divergent perturbation expansions are actually the rule in quantum physics rather than the exception (see for example the articles by ([8] Table 1) and [9], or the articles reprinted in the book by [10]). In view of the conceptual and technical problems caused by divergent series both in applied mathematics and theoretical physics, it should not be surprising that there has been extensive research on summation techniques. Recent reviews including numerous references can be found in [11,12,13,14].
The most important summation techniques are the Borel summation [15], which replaces a divergent perturbation expansion by a Laplace-type integral, and Padé approximants [16], which transform the partial sums of a (formal) power series to rational functions. There is endless research on these summation techniques. Borel summation is discussed, for example, in Costin [17], Shawyer and Watson [18], Sternin and Shatalov [19], while the probably most complete treatment of Padé approximants can still be found in the monograph by Baker and Graves-Morris [20]. In recent years, a considerable deal of work has been invested on the so-called sequence transformations, which turned out to be very useful numerical tools for the summation of divergent series. The extended recent literature, in particular on nonlinear and nonregular sequence transformations, is available in [12,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35].
Levin-type transformations [26,36,37], in particular, constitute an effective and powerful class of nonlinear sequence transformations, which differ substantially from other, better-known sequence transformations, such as, for example, the celebrated Wynn’s epsilon algorithm [38]. So many successful applications of Levin-type transformations have been achieved and described in the literature (see, for instance, the references in ([11] Appendix B), [13,39], or in the articles in [12]), that they are also treated in NIST Handbook ([2] Chapter 3.9(v) Levin’s and Weniger’s Transformations). The practical usefulness of sequence transformations as computational tools for the acceleration of convergence as well as the summation of divergent series is now established beyond doubt. However, from a purely theoretical perspective, our understanding of the convergence properties of Levin-type sequence transformations is still far from being satisfactory.
In the case of Padé, the situation is much better. For example, if the input data are the partial sums of a Stieltjes series, it can be proved rigorously that certain subsequences of the Padé table converge to the value of the corresponding Stieltjes function. More precisely, if a factorially divergent power series is a Stieltjes series and if it also satisfies the so-called Carleman condition, such a series would be Padé summable. In [13], a rigorous convergence study of the summation of the Euler series, namely E ( z ) n = 0 ( 1 ) n n ! z n , which is the paradigm of the factorially divergent power series (for a compact description, see ([13] Section 2)), was carried out by comparing the performances of Padé approximants with those of a certain class of nonlinear sequence transformations.
Readers are encouraged to go through Section 4 of [13], where an extensive summary of the most relevant properties of Levin-type transformations can be found. Levin-type transformations always admit explicit analytical expressions, which might be used as a natural starting point for building up a convergence theory. However, there is a main conceptual drawback: the technical problems occurring in any convergence studies, substantially differ among the various Levin-type transformations. For example, in [13], the convergence analysis of the summation of the Euler series turns out to be, in the case of the Weniger transformation, much easier than in the case of other Levin-type transformations. The principal reason stems in the role played by the so-called converging factor [40] (Converging factors must not be mixed up with convergence factors for summable series, which are, for instance, discussed in an article and a book by [41,42]), whose theory has been extensively discussed in articles [43,44,45,46,47,48] and a book ([49] Sections 21–26) by Dingle.
The Euler series is a Stieltjes series, and its partial sums are Padé summable to the so-called Euler integral. In the present paper, it is proved that, when the series under investigation is of the Stieltjes type, an inverse factorial representation of its converging factor can always be analytically retrieved via an algorithm based on the following points:
(i)
the converging factor of a typical Stieltjes series satisfies a first-order finite difference equation.
(ii)
there is a tight connection between inverse factorial series and recurrence relationships, the former being the natural mathematical tools for representing the solution of linear finite difference equations, similarly as solutions of differential equations are naturally expressed in terms of power series. A detailed discussion of the most important features of factorial series plus additional references can be found, for instance, in [50] and a more condensed one in ([13] Appendix B).
This paper is structured as follows: in Section 3, after the main definitions and properties of the Stieltjes series and Stieltjes functions have been briefly resumed, the recurrence rule for the converging factors of the Stieljes series is derived and formally solved through inverse factorial series expansions. In Section 4, a somewhat didactical derivation of Weniger’s transformation is proposed, in order to enlighten its connection with the factorial series representation of converging factors. Section 5 represents the core of this paper: the algorithm to generate the factorial series representation of a typical Stieltjes series converging factor is proposed and its derivation detailed. Section 6 is devoted to illustrate several examples of the application of the results established in Section 5. Finally, a few conclusive considerations and plans for future works are outlined in Section 7.
In order to improve this paper’s readability, the most technical parts of the paper, together with those containing the main definitions and notations, have been confined within several appendices, with the first of them (A) being a shorter version of ([13] Appendix B).

3. Converging Factors of Stieltjes Series

Consider an increasing real-valued function μ ( t ) defined for t [ 0 , ] , with infinitely many points of increase. The measure d μ is then positive on [ 0 , ] . Assume all positive quantities,
μ m = 0 t m d μ , m 0 ,
which will be called moments, to be finite. Then, the formal power series
m = 0 ( ) m z m + 1 μ m ,
is called a Stieltjes series. Such a series turns out to be asymptotic, for z , to the function F ( z ) defined by
F ( z ) = 0 d μ z + t , | arg ( z ) | < π ,
which is called the Stieltjes function. In other words,
F ( z ) m = 0 ( ) m z m + 1 μ m , z .
Whether such a series converges or diverges depends on the behavior of the moment sequence { μ m } m = 0 as m . Any Stieltjes function can always be expressed as the nth-order partial sum of the associated Stieltjes series (4), plus a truncation error term which has itself the form of a Stieltjes integral (see for example ([26] Theorem 13-1)), i.e.,
F ( z ) = m = 0 n ( 1 ) m z m + 1 μ m + 1 z n + 2 0 t n + 1 d μ z + t , | arg ( z ) | < π .
Such a truncation error term can always be recast as the product of the first term neglected, namely ( 1 / z ) n + 2 μ n + 1 , and a converging factor [40,49], say φ n + 1 ( z ) , which is defined as follows:
φ m ( z ) = 1 μ m 0 t m d μ t + z , m N 0 , | arg ( z ) | < π ,
in such a way that F ( z ) takes on the following form:
F ( z ) = m = 0 n ( 1 ) m z m + 1 μ m + ( 1 ) n + 1 z n + 2 μ n + 1 φ n + 1 ( z ) .
From Equation (7), it appears that the value of the Stieltjes function F ( z ) could be retrieved, in principle, starting from the knowledge of a finite number of terms of the associated Stieltjes series, provided that the corresponding converging factor is to be estimated in some way. In particular, if a sufficient approximation to the converging factor could be obtained by avoiding the explicit evaluation of the Stieltjes integral into Equation (6), then, at least in principle, all problems related to the numerical evaluation of the Stieltjes function F ( z ) via its asymptotic power series (4) would definitely be solved, even if this power series converges prohibitively slowly or not at all.
For example, the analysis in [13] was ultimately made possible thanks to the following explicit expansion obtained for the nth-order Euler series converging factor ([51] Equation (52)),
φ n ( z ) = k = 0 L k ( 1 ) ( 1 / z ) z k ! ( n ) k + 1 ,
where L k ( 1 ) ( · ) denotes a generalized Laguerre polynomial ([2] Table 18.3.1) and where the symbol ( n ) k + 1 = n ( n + 1 ) ( n + k ) denotes the so-called Pochhammer symbol. The expansion in Equation (8) is just an example of inverse factorial series in the discrete index n N 0 , which greatly facilitates the application of finite difference operators, on which Levin-type nonlinear transformations are ultimately based. For reader’s convenience, the basic definitions and properties of factorial series are briefly recalled in Appendix A. A more extensive review can be found, as previously said, in ([13] Appendix B).
It is a remarkably simple task to prove that the converging factors φ m ( z ) satisfy a first-order linear recurrence rule. To this end, it is sufficient to start from the integral representation of φ m + 1 given in Equation (6), i.e.,
φ m + 1 ( z ) = 1 μ m + 1 0 t m + 1 d μ t + z ,
which can be recast as follows:
φ m + 1 ( z ) = 1 μ m + 1 0 t m t t + z d μ = 1 μ m + 1 0 t m 1 z t + z d μ ,
so that, after taking Equations (1) and (6) into account, we have
φ m + 1 = μ m μ m + 1 1 z φ m , m 0 ,
Accordingly, the converging factor sequence { φ m } m = 0 must depend only on z and on the moment ratio sequence μ m μ m + 1 m = 0 , a fact which can hardly be obtained in terms of the sole integral representation into Equation (9). In particular, the limiting value φ = lim m φ m can be extracted directly from Equation (11), which yields
φ = ρ 0 1 + z ρ 0 ,
where
ρ 0 = lim m μ m μ m + 1 ,
is directly related to the series convergence radius.

4. Inverse Factorial Series and Weniger Transformation

That Stieltjes series converging factors have to satisfy a recurrence relationship could represent the key to understanding why the Weniger transformation turns out to be particularly fit for summing series of the Stieltjes type.
Given the nth-order converging factor φ n ( z ) , suppose there exists a sequence, say { c k ( z ) } k = 0 , whose terms are independent of n, such that the following representation holds:
φ n ( z ) = k = 0 k ! ( n ) k + 1 c k ( z ) .
Then, on substituting from Equation (14) into Equation (7), we would have
F ( z ) = m = 0 n a m ( z ) + a m + 1 ( z ) k = 0 k ! ( n + 1 ) k + 1 c k ( z ) ,
where, for brevity, the quantity a m ( z ) = ( ) m μ m / z m + 1 denotes the mth-order term of the Stieltjes series (4).
Suppose for a moment that the factorial series in Equation (14) consists only of a finite number of terms, say N. Then, it is easy to show that the correct value of F ( z ) can be retrieved within a finite number of steps and, more importantly, without explicitly knowing the sequence { c k ( z ) } k = 0 , according to the following steps:
  • Equation (15) is first recast in the form
    k = 0 N k ! ( n + 1 ) k + 1 c k ( z ) = 1 a n + 1 ( z ) F ( z ) s n ( z ) a n + 1 ( z ) ,
    where s n ( z ) = m = 0 n a m ( z ) denotes the nth-order partial sum of the series.
  • Multiply both sides of Equation (16) by ( n + 1 ) N + 1 , which gives
    ( n + 1 ) N + 1 k = 0 N k ! ( n + 1 ) k + 1 c k ( z ) = ( n + 1 ) N + 1 a n + 1 ( z ) F ( z ) ( n + 1 ) N + 1 s n ( z ) a n + 1 ( z ) .
  • The left side of Equation (17) is an Nth-order polynomial with respect to the variable n. This, in turn, implies that it can be annihilated by applying N + 1 times the first forward difference operator with respect to n, defined as
    Δ f n = f n + 1 f n .
    Accordingly, we have
    Δ N + 1 ( n + 1 ) N + 1 k = 0 N k ! ( n + 1 ) k + 1 c k ( z ) = 0 ,
    where the iterated Δ k operator turns out to be
    Δ k f n = ( 1 ) k j = 0 k ( ) j k j f n + j .
  • Finally, substitution from Equation (19) into Equation (17) eventually gives
    F ( z ) = Δ N + 1 ( n + 1 ) N + 1 s n ( z ) a n + 1 ( z ) Δ N + 1 ( n + 1 ) N + 1 a n + 1 ( z ) ,
    which proves the initial thesis.
It must be remarked again that Equation (21) represents an exact result only if the factorial expansion into Equation (14) reduces to a finite inverse factorial sum. If not, the quantity defined in the right side of Equation (21) can be interpreted as the Nth-order term of an infinite sequence, say { δ m ( n ) } m = 1 , which is built up from the partial sum sequence according to
δ m ( n ) = Δ m + 1 ( n + 1 ) m + 1 s n Δ s n Δ m + 1 ( n + 1 ) m + 1 Δ s n , m 1 , n 0 ,
where Δ s n = a n + 1 , and where, for simplicity, any dependence on the variable z has to tacitly be intended. Equation (22) defines the so-called Weniger transformation of the partial sum sequence { s m } m = 0 [26].
For the readers’ convenience, the explicit expansion of the right side of Equation (22) is given,
δ m ( n ) = j = 0 m + 1 ( 1 ) j m + 1 j ( 1 + n + j ) m + 1 s j Δ s j + 1 j = 0 m + 1 ( 1 ) j m + 1 j ( 1 + n + j ) m + 1 1 Δ s j + 1 ,
where Equation (20) has been taken into account.
What would be highly desirable is that, as m , the δ -sequence defined in Equation (22) could converge, in some sense, to the correct value F ( z ) , i.e.,
lim m δ m ( n ) = F ( z ) ,
regardless the value of the parameter n.
Presently, a general convergence theory does not exist that could validate, for a typical Stieltjes series, Equation (24). However, it seems reasonable that such a theory, assuming it does exist, should be based on the representation of the convergence factor given in (14). In the next section, an algorithm aimed at building up such a representation will be conceived.

5. On the Factorial Series Expansion of Stieltjes Series Converging Factor

Our idea consists in letting the solution of the difference Equation (11) in the form of an inverse factorial series, as follows:
φ m ( z ) = φ + k = 0 k ! ( m + β ) k + 1 c k ( z ) ,
where φ is given in Equation (12), β is a real nonnegative parameter, and { c k ( z ) } k = 0 denotes an infinite sequence independent of m (Clearly, it is expected the sequence { c k } k = 0 to depend also on β . However, for simplicity such a dependence has not been made explicit). For the moment, it will be assumed β = 0 .
Suppose now that also the sequence { μ m / μ m + 1 } m = 0 admits itself an inverse factorial expansion of the form
μ m μ m + 1 = ρ 0 + k = 0 k ! ( m ) k + 1 b k ,
with { b k } k = 0 also being a numerical sequence independent of m. As we shall see in a moment, Equation (26) can be taken for granted under rather general conditions. For example, consider again the Euler series, which is the asymptotic series (4) with F ( z ) = E ( z ) , where
E ( z ) = 0 exp ( t ) t + z d t ,
defines the so-called Euler integral. Euler series is a Stieltjes series, with μ ( t ) = 1 exp ( t ) and μ m = m ! , so that
μ m μ m + 1 = 1 m + 1 = 1 m 1 m + 1 m + 1 = 1 m 1 m ( m + 1 ) .
Equation (28) is of the form (26) with ρ 0 = 0 , b 0 = 1 , b 1 = 1 , and b k = 0 for k > 1 . It is worth noting that the factorial series expansion into Equation (26) can always be obtained, in principle, on assuming the moment ratio sequence to admit an inverse power series representation with respect to the index m,
μ m μ m + 1 = k = 0 ρ k m k = ρ 0 + k = 1 ρ k m k = ρ 0 + k = 0 ρ k + 1 m k + 1 , m 1 .
The last series expansion can always be transformed into an inverse factorial series, for example, by using the following conversion formula [50]:
μ m μ m + 1 = ρ 0 + k = 0 ρ k + 1 m k + 1 = ρ 0 + k = 0 k ! ( m ) k + 1 b k ,
where
b k = ( 1 ) k k ! μ = 0 k ( 1 ) μ s ( k , μ ) ρ μ + 1 ,
and the symbol s ( k , μ ) denotes the Stirling number of the first order (Mathematica Notation for the Stirling number: s S 1 ). In order to solve Equation (11), the following formula will also play a key role:
φ m + 1 = φ + c 0 m + k = 1 k ! ( m ) k + 1 ( c k c k 1 ) ,
whose derivation is detailed in Appendix B. Another formula which will be revealed to be of pivotal importance is that concerning the product of two factorial series. In particular, it is not difficult to prove that
k = 0 k ! ( m ) k + 1 b k k = 0 k ! ( m ) k + 1 c k = k = 1 k ! ( m ) k + 1 d k ,
with the sequence { d k } k = 1 being defined as follows:
d k = 1 k ! ν = 0 k 1 μ = 0 ν ν ! ( μ + 1 ) k 1 ν b k 1 ν c ν μ ,
or, equivalently,
d k = 1 k ! ν = 0 k 1 λ = 0 ν ν ! ( ν λ + 1 ) k 1 ν b k 1 ν c λ .
Now, all we have to do is to substitute from Equations (30), (25), and (32) into Equation (11). Nontrivial mathematical steps, detailed in Appendix C, lead to
c 1 = 0 , c 0 = b 0 ( 1 + z ρ 0 ) 2 , c k = c k 1 1 + z ρ 0 + b k ( 1 + z ρ 0 ) 2 z 1 + z ρ 0 d k , k 1 ,
which, together with Equations (26) and (34), represents the main result of the present paper.
We have proved that the nth-order truncation error of a Stieltjes series can be, under quite general hypotheses, expressed as the product of the first term of the series neglected by the nth-order partial sum, times an inverse factorial series whose single terms are defined through a triangular recurrence relation involving, apart from z, only the asymptotic expansion of the moment ratio sequence { μ m / μ m + 1 } for m .
It should be noted how, in the case of divergent series with ρ 0 = 0 , Equation (36) takes on the following form:
c 1 = 0 , c 0 = b 0 , c k = c k 1 + b k z d k , k > 1 .

6. A Few Examples of Applications

6.1. Revisiting Euler Series

In Ref. [51], the factorial representation (25) of the converging factor of the Euler series was found for the particular case of β = 0 , being
φ m ( z ) = k = 0 k ! ( m ) k + 1 z L k ( 1 ) ( z ) ,
where the symbol L k ( α ) ( · ) denotes the associate Laguerre polynomial [52]. Equation (38) was obtained by using the general approach proposed in [53]. This led to a nontrivial derivation involving the use of Bernoulli as well as exponential polynomials [51]. In the present section, Equation (38) will now be re-derived by using the results previously obtained and generalized to arbitrary positive values of β . To this end, the first, preliminary step is to find the factorial expansion of the moment ratio in terms of m + β . This can be achieved, for instance, by recasting the moment ratio as follows:
μ m μ m + 1 = 1 m + 1 = 1 m + β ( β 1 ) ,
and then by using the fundamental expansion ([54] Equation (3) on p. 77): (In the books by (Milne-Thomson [55] p. 291) and by (Paris and Kaminski [56] p. 177), this expansion is called Waring’s formula, but in Nielsen’s book ([54] Footnote 2 on p. 77), this formula is attributed to Stirling [57], and in the book by (Nörlund [58] p. 199), it is attributed to Nicole (compare also (Tweedie [59] p. 30))).
1 z w = n = 0 ( w ) n ( z ) n + 1 , Re ( z w ) > 0 ,
which gives at once
μ m μ m + 1 = k = 0 k ! ( m + β ) k + 1 ( β 1 ) k k ! .
On comparing Equation (41) with Equation (26) and on taking into account that ρ 0 = 0 , we have
b k = ( β 1 ) k k ! , k 0 , β 0 .
It is worth showing an alternative approach, based on the asymptotic expansion in Equation (29). In particular, the expanding coefficients { ρ k } k = 0 can be found by writing Equation (39) as follows:
μ m μ m + 1 = 1 m + β ( β 1 ) = 1 m + β 1 1 β 1 m + β ,
and by formally expanding the second factor as a geometric series, i.e.,
μ m μ m + 1 = k = 0 ( β 1 ) k ( m + β ) k + 1 .
On comparing Equation (44) with Equation (29), we have ρ 0 = 0 and ρ k = ( β 1 ) k 1 , for k 1 , so that Equation (31) leads to
b k = ( 1 ) k k ! μ = 0 k ( 1 ) μ s ( k , μ ) ( β 1 ) μ = ( β 1 ) k k ! ,
which proves again Equation (42).
Once the factorial expansion of the moment ratio has been obtained, it is sufficient to substitute from Equation (42) into Equation (36) to obtain, with a small amount of elementary algebra,
c 0 = 1 , c 1 = β z = L 1 ( β 1 ) ( z ) , c 2 = 1 2 β 2 + β + z 2 2 ( β + 1 ) z = L 2 ( β 1 ) ( z ) ,
and so on, for all k > 2 . Accordingly, Equation (38) becomes
φ m ( z ) = k = 0 k ! ( m + β ) k + 1 z L k ( β 1 ) ( z ) , β 0 .

6.2. Error Function

The following Stieltjes series:
m = 0 ( 1 ) m z m + 1 Γ m + 1 2 ,
where Γ ( · ) denotes the Gamma function [60], is asymptotic to the complementary error function erfc(·) [60], being
m = 0 ( 1 ) m z m + 1 Γ m + 1 2 π exp ( z ) erfc ( z ) z , z .
In order to find the factorial expansion of the corresponding converging factor, it is sufficient to recast the moment ratio as follows:
μ m μ m + 1 = Γ m + 1 2 Γ m + 3 2 = 1 m + 1 2 = 1 m 1 2 + 1 ,
and to note that it is formally identical to Equation (39), once the transformation m m 1 / 2 is applied. Accordingly, it is not difficult to prove that all steps leading to Equation (47) for the Euler series can be repeated for the asymptotic series (48), which yields
φ m ( z ) = k = 0 k ! m + β 1 2 k + 1 z L k ( β 1 ) ( z ) , m > 1 .

6.3. The Logarithm Function

Another important Stieltjes function is the natural logarithm. Consider the function L ( z ) defined as follows:
L ( z ) = 0 1 d t z + t = log 1 + 1 z .
The Stieltjes series of the form (4) asymptotic to L ( z ) , is characterized by the following moment sequence:
μ m = 1 m + 1 , m 0 ,
so that
μ m μ m + 1 = m + 2 m + 1 = 1 + 1 m + 1 .
Similarly, as for the error function, the factorial expansion of the converging factor of the logarithm can also be obtained in a surprisingly simple way. In fact, it is sufficient to compare Equation (54) with Equation (30), written with m + β in place of m and, on taking Equation (41) into account, to have
μ m μ m + 1 = 1 + k = 0 k ! ( m + β ) k + 1 ( β 1 ) k k ! ,
from which it follows that
ρ 0 = 1 b k = ( β 1 ) k k ! .
Then, on substituting from Equation (56) into Equation (36) after simple algebra, we obtain
c 0 = 1 ( 1 + z ) 2 , c 1 = β + ( β 2 ) z ( 1 + z ) 3 , c 2 = β ( 1 + β ) + 2 ( β 2 ) ( β + 1 ) z + ( β 2 ) ( β 1 ) z 2 2 ( 1 + z ) 4 ,
and so on, for all k > 2 . Remarkably, it is also possible to give Equation (57) the following closed-form expression:
c k ( z ) = ( 1 ) k ( 1 + z ) 2 P k ( 2 β k , β 1 ) z 1 z + 1 , k 0 ,
where the symbol P n ( α , β ) ( · ) denotes the Jacobi polynomial [60].

6.4. Lerch’s Trascendental Function

The logarithm function is a particular case of a more general Stieltjes function, namely the Lerch trascendental [60], whose asymptotic Stieltjes series (4) is characterized by the following moment sequence:
μ m = 1 ( m + α ) s , m 0 ,
where, for simplicity, it will be assumed α > 0 and s > 1 . The moment ratio then turns out to be
μ m μ m + 1 = 1 + 1 m + α s = 1 + 1 m + β + ( α β ) s , m 0 ,
whose inverse factorial expansion can be obtained by using Equations (29)–(31). To this end, consider first the case β = 0 . Then, the asymptotic inverse power expansion with respect to m gives
ρ 0 = 1 , ρ 1 = s , ρ 2 = 1 2 s ( s 1 ) s α , ρ 3 = 1 6 s ( s 1 ) ( s 2 ) s ( s 1 ) α + s α 2 , ρ 4 = 1 24 s ( s 1 ) ( s 2 ) ( s 3 ) 1 2 s ( s 1 ) ( s 2 ) α + 3 2 s ( s 1 ) α 2 s α 3 ,
Accordingly, on using Equation (31) again, the expanding coefficients b k of the factorial expansion of μ m μ m + 1 turn out to be
b 0 = s , b 1 = 1 2 s ( s 1 2 α ) , b 2 = 1 12 s ( s 2 1 6 s α + 6 α 2 ) , b 3 = 1 144 s ( s ( s ( 6 + s ) 1 ) 12 s ( s + 3 ) α + 36 ( s + 1 ) α 2 24 α 3 6 ) ,
and our recursive algorithm provides the following analytical expressions of the first terms of the sequence { c k } k = 0 :
c 0 = s ( 1 + z ) 2 , c 1 = ( s ( ( s + 1 ) ( z 1 ) + 2 ( 1 + z ) α ) 2 ( 1 + z ) 3 , c 2 = s ( ( 1 + s ) ( 5 + s + s ( z 4 ) z z ( 8 + z ) ) + 6 ( s ( z 1 ) 2 ) ( 1 + z ) α + 6 ( 1 + z ) 2 α 2 ) 12 ( 1 + z ) 4 ,
Note that, in order to recover the more general factorial expansion of the converging factor (25) valid for a typical β > 0 , it is sufficient, according to Equation (60), to replace α by α β into Equation (63). Differently from Equation (57), we were not able to give Equation (63) a closed form expression.

7. Conclusions

Although the practical usefulness of divergent series has always been more important than concerns about mathematical rigor, especially in physics, a theoretical understanding of the mathematical machinery behind Levin-type sequence transformations would be highly desirable. An important achievement in this direction has been established in [13], where it was rigorously proved that a particular type of Levin-type transformation, the Weniger transformation [26], is not only able to resum the Euler series, but also that it turns out to be “exponentially faster” than Padé approximants in accomplishing such a task. The Euler series is the paradigm of factorial divergence. Several other series expansions occurring in mathematical physics are factorially or even hyperfactorially diverging. The unquestionable superiority of the Weniger transformation over Padé in resumming the Euler series proved in [13] was not limited to such a specific case, as it was shown in the past [61,62,63,64,65].
Levin-type transformations in general, and Weniger’s transformation in particular, are able to dramatically accelerate the convergence of slowly convergent series as well as to resum wildly divergent series, some of them even not resummable via Padé. The retrieving capabilities of Weniger’s transformation are strictly related to the possibility of representing the series converging factor of any finite order as a suitable inverse factorial series.
In the present paper, which is the natural continuation of Ref. [13], it has been proved that the nth-order converging factor of a Stieltjes series can be expressed through an inverse factorial series, whose terms are uniquely determined, starting from the knowledge of the series moment ratio, through a linear recurrence relation. Our constructive proof has the form of an algorithm, developed starting from the first-order difference equation which has to be satisfied by the convergent factors of any Stieltjes series. A few examples of applications of our algorithm have also been presented, in order to show its effectiveness and implementation ease.
The class of Stieltjes series is of fundamental interest in mathematical physics. Recently, it was rigorously shown that the factorially divergent perturbation expansion for the energy eigenvalue of the PT -symmetric Hamiltonian H ( λ ) = p 2 + 1 4 x 2 + i λ x 3 is a Stieltjes series [66,67], thus proving a conjecture formulated ten years before by Bender and Weniger [68]. A much more recent result is the proof that the character of the celebrated Bessel solution of Kepler’s equation [69] is also Stieltjes [70]. Weniger’s transformation turned out to be successful, in the past, for resumming extremely violently divergent perturbation expansions [71,72,73,74], where Padé approximants were not powerful enough to achieve anything substantial, as in the case of Rayleigh-Schrödinger perturbation series for the sextic anharmonic oscillator, or even not able to sum, as for the more challenging octic case [75]. There is no doubt that a better theoretical understanding of the summation of hyper-factorially divergent expansions via Weniger’s transformation would also be highly desirable.
Ernst Joachim Weniger and I were strongly convinced that the results obtained in [13] could have been a preliminary step toward the development of a general convergence theory of the summation of Stieltjes series with Levin-type transformations in place of Padé approximants. What has been shown here is a small leap in this direction. The present work could also lead to a hypothetical general framework for the analysis of the convergence rate of Weniger’s transformation. A possible path in this direction (I must thank one of the anonymous reviewer for such an interesting suggestion), i.e., to obtain an estimate for the error of a delta transformation, could be to look at the first term of the inverse factorial series not eliminated at the m-th step into the δ -sequence of Equation (22), and to try approximating the remainder term of the transformation in some compact asymptotic form. This is not an easy task at all but worthy of being pursued in the future.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

I wish to thank one of the anonymous reviewers for his/her criticisms and suggestions. I wish to thank Turi Maria Spinozzi for his useful comments and help.

Conflicts of Interest

The author declares no conflicts of interest.

Appendix A. Basic Properties of Factorial Series

Let Ω ( z ) be a function that vanishes as z + . A factorial series for Ω ( z ) is an expansion in terms of inverse Pochhammer symbols of the following kind:
Ω ( z ) = b 0 z + b 1 1 ! z ( z + 1 ) + b 2 2 ! z ( z + 1 ) ( z + 2 ) + = ν = 0 b ν ν ! ( z ) ν + 1 .
In general, Ω ( z ) will have simple poles at z = m with m N 0 .
Factorial series had been employed already in Stirling’s classic book Methodus Differentialis [57,76], which appeared in 1730. In the nineteenth and the early twentieth century, the theory of factorial series was fully developed by a variety of authors. Fairly complete surveys of the older literature as well as thorough treatments of their properties can be found in books on finite differences by Meschkowski [77], Milne-Thomson [55], Nielsen [54], and Nörlund [58], Nörlund [78], Nörlund [79]. Factorial series are also discussed in the books by Knopp [80] and Nielsen [81] on infinite series.
Compared to inverse power series, factorial series possess vastly superior convergence properties. If we use (see for example (Olver et al. [2] Equation (5.11.12)))
Γ ( z + a ) / Γ ( z + b ) = z a b 1 + O ( 1 / z ) , z ,
we obtain the asymptotic estimate n ! / ( z ) n + 1 = O ( n z ) as n . Thus, the factorial series (A1) converges with the possible exception of the points z = m with m N 0 if and only if the associated Dirichlet series Ω ¯ ( z ) = n = 1 b n / n z converges (see for example ([80] p. 262) or ([82] p. 167)). Accordingly, a factorial series converges for sufficiently large Re ( z ) even if the reduced series coefficients b n in (A1) grow like a fixed power n α with Re ( α ) > 0 as n . Thus, a factorial series can converge even if its effective coefficients b n n ! diverge quite wildly as n .
Factorial series can also be expressed in terms of the beta function B ( x , y ) = Γ ( x ) Γ ( y ) / Γ ( x + y ) ([2] Equation (5.12.1)). Because of B ( z , n + 1 ) = n ! / ( z ) n + 1 , the factorial series in Equation (A1) can be interpreted as an expansion in terms of special beta functions (see for instance (Milne-Thomson [55] p. 288) or (Paris and Kaminski [56] [p. 175])):
Ω ( z ) = n = 0 b n B ( z , n + 1 ) .
The beta function possesses the integral representation (see for example ([2] Equation (5.12.1)))
B ( x , y ) = 0 1 t x 1 ( 1 t ) y 1 d t , Re ( x ) , Re ( y ) > 0 ,
which implies
B ( z , n + 1 ) = n ! ( z ) n + 1 = 0 1 t z 1 ( 1 t ) n d t , Re ( z ) > 0 n N 0 .
If we insert this into Euqation (A3) and interchange integration and summation, we obtain the following very important integral representation ([54] Section I on p. 244):
Ω ( z ) = 0 1 t z 1 φ Ω ( t ) d t , Re ( z ) > 0 ,
φ Ω ( t ) = n = 0 b n ( 1 t ) n .
Often, the properties of Ω ( z ) can be studied more easily via this integral representation than via the defining factorial series (A1) (see for example ([54] Ch. XVII)).

Appendix B. Proof of Equation (32)

First of all, φ m + 1 is recast as follows:
φ m + 1 = φ m + Δ φ m .
Then, on applying the forward difference operator Δ to both sides of Equation (25), we have
Δ φ m = Δ k = 0 k ! ( m ) k + 1 c k = k = 0 k ! c k Δ 1 ( m ) k + 1 ,
where (see ([26] Equation (13.2-12)))
Δ 1 ( m ) k + 1 = k + 1 ( m ) k + 2 ,
so that
Δ φ m = k = 1 k ! ( m ) k + 1 c k 1 .
Then, on taking again Equation (A7) into account, from Equations (25) and (A10), Equation (32) follows after simple rearranging.

Appendix C. Proof of Equation (36)

We start by recasting Equation (11) as follows:
φ m + 1 = μ m μ m + 1 z μ m μ m + 1 φ m ,
which, on taking Equations (25), (30), and (32) into account, can be recast as follows:
φ + k = 0 k ! ( m ) k + 1 Δ c k 1 = ρ 0 + k = 0 k ! ( m ) k + 1 b k z ρ 0 + k = 0 k ! ( m ) k + 1 b k φ + k = 0 k ! ( m ) k + 1 c k ,
where, for simplicity, it has been assumed formally c 1 = 0 . Equation (A12) can be rearranged, on taking Equation (12) into account, as follows:
k = 0 k ! ( m ) k + 1 Δ c k 1 = 1 1 + z ρ 0 k = 0 k ! ( m ) k + 1 b k z ρ 0 k = 0 k ! ( m ) k + 1 c k z k = 0 k ! ( m ) k + 1 b k k = 0 k ! ( m ) k + 1 c k ,
Finally, on using Equation (33) into Equation (A13), after simple algebra, we obtain
c 0 = b 0 1 + z ρ 0 z ρ 0 c 0 , Δ c k 1 = b k 1 + z ρ 0 z ρ 0 c k z d k , k 1 ,
from which Equation (36) follows.

References

  1. Hardy, G.H. Divergent Series; Clarendon Press: Oxford, UK, 1949. [Google Scholar]
  2. Olver, F.W.J.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. (Eds.) NIST Handbook of Mathematical Functions; Cambridge U. P.: Cambridge, UK, 2010. [Google Scholar]
  3. Poincaré, H. Sur les intégrales irrégulières des équations linéaires. Acta Math. 1886, 8, 295–344. [Google Scholar] [CrossRef]
  4. Dyson, D.J. Divergence of perturbation theory in quantum electrodynamics. Phys. Rev. 1952, 85, 32–33. [Google Scholar] [CrossRef]
  5. Bender, C.M.; Wu, T.T. Anharmonic oscillator. Phys. Rev. 1969, 184, 1231–1260. [Google Scholar] [CrossRef]
  6. Bender, C.M.; Wu, T.T. Large-order behavior of perturbation theory. Phys. Rev. Lett. 1971, 27, 461–465. [Google Scholar] [CrossRef]
  7. Bender, C.M.; Wu, T.T. Anharmonic oscillator. II. A study in perturbation theory in large order. Phys. Rev. D 1973, 7, 1620–1636. [Google Scholar] [CrossRef]
  8. Fischer, J. On the role of power expansions in quantum field theory. Int. J. Mod. Phys. A 1997, 12, 3625–3663. [Google Scholar] [CrossRef]
  9. Suslov, I.M. Divergent perturbation series. J. Exp. Theor. Phys. (JETP) 2005, 100, 1188–1234. [Google Scholar] [CrossRef]
  10. Le Guillou, J.C.; Zinn-Justin, J. (Eds.) Large-Order Behaviour of Perturbation Theory; North-Holland: Amsterdam, The Netherlands, 1990. [Google Scholar]
  11. Weniger, E.J. On the analyticity of Laguerre series. J. Phys. A 2008, 41, 425207-1–425207-43. [Google Scholar] [CrossRef]
  12. Brezinski, C.; Redivo Zaglia, M.; Weniger, E.J. Special Issue: Approximation and extrapolation of convergent and divergent sequences and series (CIRM, Luminy–France, 2009). Appl. Numer. Math. 2010, 60, 1183–1464. [Google Scholar] [CrossRef]
  13. Borghi, R.; Weniger, E.J. Convergence analysis of the summation of the factorially divergent Euler series by Padé approximants and the delta transformation. Appl. Numer. Math. 2015, 94, 149–178. [Google Scholar] [CrossRef]
  14. Borghi, R. Computational optics through sequence transformations. In Progress in Optics; Visser, T.D., Ed.; Academic Press: Amsterdam, The Netherlands, 2016; Volume 61, Chapter 1; pp. 1–70. [Google Scholar]
  15. Borel, E. Mémoires sur les séries divergentes. Ann. Sci. Éc. Norm. Sup. 1899, 16, 9–136. [Google Scholar] [CrossRef]
  16. Padé, H. Sur la représentation approachée d’une fonction par des fractions rationelles. Ann. Sci. Éc. Norm. Sup. 1892, 9, 3–93. [Google Scholar] [CrossRef]
  17. Costin, O. Asymptotics and Borel Summability; Chapman & Hall/CRC: Boca Raton, FL, USA, 2009. [Google Scholar]
  18. Shawyer, B.; Watson, B. Borel’s Method of Summability; Oxford U. P.: Oxford, UK, 1994. [Google Scholar]
  19. Sternin, B.Y.; Shatalov, V.E. Borel-Laplace Transform and Asymptotic Theory; CRC Press: Boca Raton, FL, USA, 1996. [Google Scholar]
  20. Baker, G.A., Jr.; Graves-Morris, P. Padé Approximants, 2nd ed.; Cambridge U. P.: Cambridge, UK, 1996. [Google Scholar]
  21. Brezinski, C. Accélération de la Convergence en Analyse Numérique; Springer: Berlin, Germany, 1977. [Google Scholar]
  22. Brezinski, C. Algorithmes d’Accélération de la Convergence—Étude Numérique; Éditions Technip: Paris, France, 1978. [Google Scholar]
  23. Brezinski, C.; Redivo Zaglia, M. Extrapolation Methods; North-Holland: Amsterdam, The Netherlands, 1991. [Google Scholar]
  24. Sidi, A. Practical Extrapolation Methods; Cambridge U. P.: Cambridge, UK, 2003. [Google Scholar]
  25. Wimp, J. Sequence Transformations and Their Applications; Academic Press: New York, NY, USA, 1981. [Google Scholar]
  26. Weniger, E.J. Nonlinear sequence transformations for the acceleration of convergence and the summation of divergent series. Comput. Phys. Rep. 1989, 10, 189–371. [Google Scholar] [CrossRef]
  27. Homeier, H.H. Analytical and numerical studies of the convergence behavior of the j transformation. J. Comput. Appl. Math. 1996, 69, 81–112. [Google Scholar] [CrossRef]
  28. Weniger, E.J. An introduction to the topics presented at the conference “Approximation and extrapolation of convergent and divergent sequences and series” CIRM Luminy: 28 September–2 October 2009. Appl. Numer. Math. 2010, 60, 1184–1187. [Google Scholar] [CrossRef]
  29. Aksenov, S.V.; Savageau, M.A.; Jentschura, U.D.; Becher, J.; Soff, G.; Mohr, P.J. Application of the combined nonlinear-condensation transformation to problems in statistical analysis and theoretical physics. Comput. Phys. Commun. 2003, 150, 1–20. [Google Scholar] [CrossRef]
  30. Bornemann, F.; Laurie, D.; Wagon, S.; Waldvogel, J. The SIAM 100-Digit Challenge: A Study in High-Accuracy Numerical Computing; Society of Industrial and Applied Mathematics: Philadelphia, PA, USA, 2004. [Google Scholar]
  31. Gil, A.; Segura, J.; Temme, N.M. Numerical Methods for Special Functions; SIAM: Philadelphia, PA, USA, 2007. [Google Scholar]
  32. Caliceti, E.; Meyer-Hermann, M.; Ribeca, P.; Surzhykov, A.; Jentschura, U.D. From useful algorithms for slowly convergent series to physical predictions based on divergent perturbative expansions. Phys. Rep. 2007, 446, 1–96. [Google Scholar] [CrossRef]
  33. Temme, N.M. Numerical aspects of special functions. Acta Numer. 2007, 16, 379–478. [Google Scholar] [CrossRef]
  34. Gil, A.; Segura, J.; Temme, N.M. Basic methods for computing special functions. In Recent Advances in Computational and Applied Mathematics; Simos, T.E., Ed.; Springer: Dordrecht, The Netherlands, 2011; pp. 67–121. [Google Scholar]
  35. Trefethen, L.N. Approximation Theory and Approximation Practice; SIAM: Philadelphia, PA, USA, 2013. [Google Scholar]
  36. Levin, D. Development of non-linear transformations for improving convergence of sequences. Int. J. Comput. Math. B 1973, 3, 371–388. [Google Scholar] [CrossRef]
  37. Weniger, E.J. Mathematical properties of a new Levin-type sequence transformation introduced by Čížek, Zamastil, and Skála. I. Algebraic theory. J. Math. Phys. 2004, 45, 1209–1246. [Google Scholar] [CrossRef]
  38. Wynn, P. On a device for computing the em(Sn) transformation. Math. Tables Aids Comput. 1956, 10, 91–96. [Google Scholar] [CrossRef]
  39. Chang, X.K.; He, Y.; Hu, X.B.; Sun, J.Q.; Weniger, E.J. New generalizations of Wynn’s epsilon and rho algorithm by solving finite difference equations in the transformation order. Numer. Algor. 2020, 83, 593–627. [Google Scholar] [CrossRef]
  40. Airey, J.R. The “converging factor” in asymptotic series and the calculation of Bessel, Laguerre and other functions. Philos. Mag. 1937, 24, 521–552. [Google Scholar] [CrossRef]
  41. Moore, C.N. On the introduction of convergence factors into summable series and summable integrals. Trans. Am. Math. Soc. 1907, 8, 299–330. [Google Scholar] [CrossRef]
  42. Moore, C.N. Summable Series and Convergence Factors; American Mathematical Society: Providence, RI, USA, 2008. [Google Scholar]
  43. Dingle, R.B. Asymptotic expansions and converging factors I. General theory and basic converging factors. Proc. Royal Soc. Lond. A 1958, 244, 456–475. [Google Scholar]
  44. Dingle, R.B. Asymptotic expansions and converging factors II. Error, Dawson, Fresnel, exponential, sine and cosine, and similar integrals. Proc. Royal Soc. Lond. A 1958, 244, 476–483. [Google Scholar]
  45. Dingle, R.B. Asymptotic expansions and converging factors III. Gamma, psi and polygamma functions, and Fermi-Dirac and Bose-Einstein integrals. Proc. Royal Soc. Lond. A 1958, 244, 484–490. [Google Scholar]
  46. Dingle, R.B. Asymptotic expansions and converging factors IV. Confluent hypergeometric, parabolic cylinder, modified Bessel, and ordinary Bessel functions. Proc. Royal Soc. Lond. A 1959, 249, 270–283. [Google Scholar]
  47. Dingle, R.B. Asymptotic expansions and converging factors V. Lommel, Struve, modified Struve, Anger and Weber functions, and integrals of ordinary and modified Bessel functions. Proc. Royal Soc. Lond. A 1959, 249, 284–292. [Google Scholar]
  48. Dingle, R.B. Asymptotic expansions and converging factors. VI. Application to physical prediction. Proc. Royal Soc. Lond. A 1959, 249, 293–295. [Google Scholar]
  49. Dingle, R.B. Asymptotic Expansions: Their Derivation and Interpretation; Academic Press: London, UK, 1973. [Google Scholar]
  50. Weniger, E.J. Summation of divergent power series by means of factorial series. Appl. Numer. Math. 2010, 60, 1429–1441. [Google Scholar] [CrossRef]
  51. Borghi, R. Asymptotic and factorial expansions of Euler series truncation errors via exponential polynomials. Appl. Numer. Math. 2010, 60, 1242–1250. [Google Scholar] [CrossRef]
  52. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions; National Bureau of Standards: Washington, DC, USA, 1972. [Google Scholar]
  53. Weniger, E.J. Asymptotic approximations to truncation errors of series representations for special functions. In Algorithms for Approximation; Iske, A., Levesley, J., Eds.; Springer: Berlin, Germany, 2007; pp. 331–348. [Google Scholar]
  54. Nielsen, N. Die Gammafunktion; Chelsea: New York, NY, USA, 1965. [Google Scholar]
  55. Milne-Thomson, L.M. The Calculus of Finite Differences; Chelsea: New York, NY, USA, 1981. [Google Scholar]
  56. Paris, R.B.; Kaminski, D. Asymptotics and Mellin-Barnes Integrals; Cambridge U.P.: Cambridge, UK, 2001. [Google Scholar]
  57. Stirling, J. Methodus Differentialis sive Tractatus de Summatione et Interpolatione Serierum Infinitarum; Strahan: London, UK, 1730. [Google Scholar]
  58. Nörlund, N.E. Vorlesungen über Differenzenrechnung; Chelsea: New York, NY, USA, 1954. [Google Scholar]
  59. Tweedie, C. Nicole’s contributions to the foundations of the calculus of finite differences. Proc. Edinb. Math. Soc. 1917, 36, 22–39. [Google Scholar] [CrossRef]
  60. Olver, F.W.J.; Daalhuis, A.B.O.; Lozier, D.W.; Schneider, B.I.; Boisvert, R.F.; Clark, C.W.; Miller, B.R.; Saunders, B.V.; Cohl, H.S.; McClain, M.A. (Eds.) NIST Digital Library of Mathematical Functions; Release 1.0.26 of 15 March 2020. Available online: http://dlmf.nist.gov/ (accessed on 1 June 2024).
  61. Weniger, E.J. Computation of the Whittaker function of the second kind by summing its divergent asymptotic series with the help of nonlinear sequence transformations. Comput. Phys. 1996, 10, 496–503. [Google Scholar] [CrossRef]
  62. Jentschura, U.D.; Mohr, P.J.; Soff, G.; Weniger, E.J. Convergence acceleration via combined nonlinear-condensation transformations. Comput. Phys. Commun. 1999, 116, 28–54. [Google Scholar] [CrossRef]
  63. Jentschura, U.D.; Mohr, P.J.; Soff, G. Calculation of the electron self-energy for low nuclear charge. Phys. Rev. Lett. 1999, 82, 53–56. [Google Scholar] [CrossRef]
  64. Jentschura, U.D.; Becher, J.; Weniger, E.J.; Soff, G. Resummation of QED perturbation series by sequence transformations and the prediction of perturbative coefficients. Phys. Rev. Lett. 2000, 85, 2446–2449. [Google Scholar] [CrossRef] [PubMed]
  65. Borghi, R.; Santarsiero, M. Summing Lax series for nonparaxial beam propagation. Opt. Lett. 2003, 28, 774–776. [Google Scholar] [CrossRef] [PubMed]
  66. Grecchi, V.; Maioli, M.; Martinez, A. Padé summability of the cubic oscillator. J. Phys. A 2009, 42, 425208-1–425208-17. [Google Scholar] [CrossRef]
  67. Grecchi, V.; Martinez, A. The spectrum of the cubic oscillator. Commun. Math. Phys. 2013, 319, 479–500. [Google Scholar] [CrossRef]
  68. Bender, C.M.; Weniger, E.J. Numerical evidence that the perturbation expansion for a non-Hermitian PT -symmetric Hamiltonian is Stieltjes. J. Math. Phys. 2001, 42, 2167–2183. [Google Scholar] [CrossRef]
  69. Colwell, P. Solving Kepler’s Equation Over Three Centuries; Willmann-Bell: Richmond, VA, USA, 1993. [Google Scholar]
  70. Borghi, R. On the Bessel Solution of Kepler’s Equation. Mathematics 2024, 12, 154. [Google Scholar] [CrossRef]
  71. Weniger, E.J. A convergent renormalized strong coupling perturbation expansion for the ground state energy of the quartic, sextic, and octic anharmonic oscillator. Ann. Phys. 1996, 246, 133–165. [Google Scholar] [CrossRef]
  72. Weniger, E.J. Construction of the strong coupling expansion for the ground state energy of the quartic, sextic and octic anharmonic oscillator via a renormalized strong coupling expansion. Phys. Rev. Lett. 1996, 77, 2859–2862. [Google Scholar] [CrossRef] [PubMed]
  73. Weniger, E.J.; Čížek, J.; Vinette, F. Very accurate summation for the infinite coupling limit of the perturbation series expansions of anharmonic oscillators. Phys. Lett. A 1991, 156, 169–174. [Google Scholar] [CrossRef]
  74. Weniger, E.J.; Čížek, J.; Vinette, F. The summation of the ordinary and renormalized perturbation series for the ground state energy of the quartic, sextic, and octic anharmonic oscillators using nonlinear sequence transformations. J. Math. Phys. 1993, 34, 571–609. [Google Scholar] [CrossRef]
  75. Graffi, S.; Grecchi, V. Borel summability and indeterminacy of the Stieltjes moment problem: Application to the anharmonic oscillators. J. Math. Phys. 1978, 19, 1002–1006. [Google Scholar] [CrossRef]
  76. Stirling, J. The Differential Method, or, a Treatise Concerning the Summation and Interpolation of Infinite Series; Strahan: London, UK, 1749. [Google Scholar]
  77. Meschkowski, H. Differenzengleichungen; Vandenhoek & Rupprecht: Göttingen, Germany, 1959. [Google Scholar]
  78. Nörlund, N.E. Leçons sur les Séries d’Interpolation; Gautier-Villars: Paris, France, 1926. [Google Scholar]
  79. Nörlund, N.E. Leçons sur les Équations Linéaires aux Différences Finies; Gautier-Villars: Paris, France, 1929. [Google Scholar]
  80. Knopp, K. Theorie und Anwendung der Unendlichen Reihen, 5th ed.; Springer: Berlin, Germany, 1964. [Google Scholar]
  81. Nielsen, N. Lehrbuch der Unendlichen Reihen; Teubner: Leipzig/Berlin, Germany, 1909. [Google Scholar]
  82. Landau, E. Über die Grundlagen der Theorie der Fakultätenreihen. Sitzungsb. Königl. Bay. Akad. Wissensch. München Math.-Phys. Kl. 1906, 36, 151–218. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Borghi, R. Factorial Series Representation of Stieltjes Series Converging Factors. Mathematics 2024, 12, 2330. https://doi.org/10.3390/math12152330

AMA Style

Borghi R. Factorial Series Representation of Stieltjes Series Converging Factors. Mathematics. 2024; 12(15):2330. https://doi.org/10.3390/math12152330

Chicago/Turabian Style

Borghi, Riccardo. 2024. "Factorial Series Representation of Stieltjes Series Converging Factors" Mathematics 12, no. 15: 2330. https://doi.org/10.3390/math12152330

APA Style

Borghi, R. (2024). Factorial Series Representation of Stieltjes Series Converging Factors. Mathematics, 12(15), 2330. https://doi.org/10.3390/math12152330

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop