1. Introduction
The problem of extrapolating asymptotic series derived for small variables to finite and even very large values of variables is constantly met in various branches of science, such as physics, chemistry, economics, applied mathematics, etc. Different approaches have been suggested to cure this problem, e.g., Padé summation [
1], Borel summation [
2], Shanks transformation [
3], hypergeometric Meijer approximants [
4,
5], the renormalization group [
6,
7], and others [
8,
9,
10,
11,
12,
13,
14]. These techniques can provide good approximations for finite variables, although they are not applicable for the case of variables tending to infinity.
In the present paper, we describe an original approach based on the transparent physical notions of optimization and self-similarity combined with Borel-type summation. The layout of the paper is as follows. In
Section 2, we explain the main ideas: (i) how to make a divergent asymptotic series convergent and how the convergence is optimized by introducing control functions; (ii) how to transform a sequence of approximants into a dynamical system leading to the property of self-similarity; (iii) that the solutions of the self-similar equation of motion extrapolate the given asymptotic series serving as a boundary condition into self-similar approximants that are valid for arbitrary values of variables. In
Section 3, we combine the self-similar approximants with Borel-type summation and demonstrate the ways of introducing control functions or control parameters, defined by the minimal difference and minimal derivative optimization conditions. Several admissible Borel-type transforms are discussed. In
Section 4, we discuss the specifics of the critical amplitude calculations. In
Section 5, we illustrate the discussed methods by concrete examples.
Section 6 concludes the paper.
In our approach to the resummation problem, we are guided by the requirement for the independence of observables from the approximations, transformations, parameters, etc., introduced in the course of the analysis. Nature should not be aware of our difficulties of understanding it. In turn, we should be respectful and introduce minimal necessary assumptions about it. Since the description of the sought function is available to us only by means of truncated series, we need to compensate for the lack of knowledge of the true coefficients by adding some natural assumptions.
First, the initial truncated series is to be transformed so that, instead of the available divergent truncated series, a supposedly better-behaving truncated series is considered. Second, to the transformed series, we apply the so-called approximations with power law behavior at infinity. At this step, we also reconstruct the behavior of the series. Third, the guiding principle of self-similarity leads to self-similar roots and self-similar factor approximants. Fourth, from a technical viewpoint, it is enormously helpful for the expressions for the critical amplitudes to become explicit so that they can be factored into the parts arising from the approximants and from the inverse transformation. Fifth, after an optimization with the minimal difference and (or) minimal derivative conditions, we study the numerical convergence of the sequences of approximations for the sought amplitudes. Besides the numerical convergence, we are also guided, when appropriate, by the rigorous results for the convergence of the diagonal Padé approximants obtained by Gonchar.
4. Critical Amplitudes from Fractional Iterations
Consider the case where one has to find explicitly a real, sign-definite, positive-valued function
of a real variable
x. In addition, the function possesses the power law asymptotic behavior (
4). The critical exponent
is known. The case of a negative
will be considered by studying the inverse of
. Let us find the amplitude
B.
We consider the case where it is impossible to find the sought function
explicitly and exactly from some given equations. We show that it is possible to reconstruct the large-variable amplitudes from the given truncated asymptotic expansions given in the form (
3).
Several methods of finding effective sums of the truncated series (
1)–(
3) exist, based on the ideas of Borel, Mittag–Leffler, and Hardy [
9,
44,
46,
50,
51,
52,
53,
54,
55,
56]. The well-known method is that of Padé [
57,
58,
59]. The hypergeometric approximants [
60,
61,
62,
63] can be employed instead of the Padé approximants. More complex are the hypergeometric Meijer approximants [
62,
63].
Bear in mind that we are interested in accurate analytical calculation of the amplitude
B in the expression (
4) [
64,
65,
66]. To this end, we need, of course, to find effective sums for (
1)–(
3), at least in a general form. The specifics of our techniques and the general knowledge of the large-
x asymptotic behavior (
4) are sufficient in order to find an analytic expression for the amplitudes and to proceed by applying the minimal difference and minimal derivative conditions given by Equations (
74)–(
78).
Let us briefly recapitulate the main tenets of
Section 3.6 needed for the following applications to concrete problems. First, we set all
values to zero. This simplification brings us back to a manageable number of control parameters. Below, we consider
u, which is the order of the operator in (
69), as a continuous control parameter. It is possible to define the multidimensional integrals required for calculating the critical amplitudes for integers
p. Introducing continuous
p values means a smooth interpolation between the values of the integral for discrete
p values [
44,
49]. The number of iterations
p can also be regarded as a continuous control parameter [
44]. The parameters
and
p ought to be found from the optimization conditions. The general-type optimization conditions are the minimal difference and minimal derivative conditions (
74)–(
76), expressed as the conditions for the critical amplitudes. The conditions are equivalent and both are applied below.
We expect that the solution to the optimization problem exists and is unique. There are three main realizations of the Fractional Borel Summation with factor approximants described in
Section 2.5. Often, we simply use the term “Borel summation with factor approximants”, but specify the optimization techniques applied in each particular case. Below we use the following three methods of finding the critical amplitudes.
The first method is based on introducing the fractional order
u of the differentiation operator. The parameter has to be found from the optimization procedure, while
p is fixed. The optimization procedure, where only the parameter
u is considered and
p is fixed, is called Fractional Borel Summation with
u-optimization [
49].
The second method is based on the optimization of the other parameter
p, where the number of iterations extends to arbitrary real numbers from the original integers (see
Section 3.3 and the paper [
44]). This has to be found from
p-optimization, while
u is fixed. The optimization procedure where only the parameter
p is considered and
u is fixed is called Fractional Borel Summation with
p-optimization [
49]. The optimization is performed according to formulas analogous to (
74)–(
76) with straightforward replacement of the parameter
u by the parameter
p.
The third method, Fractional Borel-light or self-similar combined approximants is explained in great detail in
Section 2.6 and
Section 4.3. The method was suggested in the paper [
49], following the main ideas expressed earlier in [
43,
46]. It is based on optimization with minimal derivative or minimal difference conditions of the amplitude
(or with the marginal amplitudes
), either with respect to the fractional
u or fractional
p, with subsequent correction of the marginal amplitudes with the diagonal Padé approximants [
1,
43,
46,
49].
We employ some useful formulas for the factor approximants given above in the general form in
Section 3.6. These are adjusted to the calculations of the amplitudes in
Section 4.1. We also employ the iterated roots explained in
Section 2.4, and diagonal Padé approximants for odd and even number of the terms
k in the truncations [
58].
Factor approximants have the advantage of generality, since the case of is included into the consideration automatically, unlike the case of iterated roots, where such a case should be treated individually. On the other hand, iterated roots are very user-friendly and can be treated analytically in rather high orders. The problem of optimization can be formulated with rather high orders of the perturbation theory, while for the factors, only low-order optimizations can be treated analytically.
The diagonal Padé approximants are routinely extendable to very high orders. Fractional Borel Summation with iterated roots was previously considered in detail [
49], while the diagonal, odd, and even Padé approximants are discussed in [
58].
For completeness and convenience, below, we give some formulas that are required for actual computations of the critical amplitudes when the index
is known. Two types of approximants are discussed in this context, since the iterated roots, presented in
Section 2.4, are well adapted to the calculation of the critical amplitudes and do not require any specific adjustments.
4.1. Critical Amplitudes from Factors and Padé Approximants
It is both convenient and natural to extrapolate the asymptotic series with power law behavior (
4) by means of the self-similar factor approximants [
35,
36]. Let us fix the inner sum of the parameters
in the Formula (
26) to the known index
, so that
in all orders. Then, the critical amplitudes can be found by extrapolating the series (
1)–(
3) to the form (
26). Here
At large
x values, the factor approximants behave as
and the amplitude is given by the Formula (
29). For the even orders
k, one of the
values can be set to one. Such a restriction does not change the critical index and only influences the value of the amplitude. From a technical standpoint, it is rather difficult to optimize factor-based methods in high orders, although lower orders can easily be optimized.
4.2. Modified Padé Approximations
Alternatively, in place of the factors, one can apply the well-known Padé approximants
[
9]. The Padé approximants can be adapted to calculate the amplitudes at infinity in the expression (
4) by means of some transformations. In the case of odd
, one should study the following transformed series for the sought function
,
The following modified Padé approximants
are defined for even cases, where
is a non-negative integer.
In the odd case,
. The approximants (
82) do comply with the power law (
4) at
. One can relatively easily find the sequence of approximations
for the amplitudes. In the case of even
, the following modified-even Padé approximants
are defined [
58]. Here, the corrector
was introduced to ensure the correct asymptotic behavior [
58]. Also,
represents the transformed truncated series, which can be approached again with the diagonal Padé approximants. More details on the application of the modified Padé approximants for the Borel summation can be found in [
58].
The amplitudes sought at infinity can be found as follows:
with
Here,
, is a positive integer. In the even case,
. It is not impossible but is very difficult from a technical standpoint to optimize Padé-based methods, especially at very high orders.
The methods of the factor, root, and Padé approximants can be applied individually or together in some hybrid form to calculate the marginal amplitude
C appearing in the course of the Borel transformation. All of the mentioned approximations can be applied directly to the series (
1)–(
3) under the asymptotic condition (
1) and produce the estimates for the sought amplitude
B. Only the iterated roots shown in
Section 3.3 are seamlessly defined for all
k values. The other two approximations use two different definitions for the odd and even cases.
4.3. Critical Amplitudes from Hybrids of Factors and Padé Approximations
The sought function can be reconstructed from the transform
directly with the help of corrected Padé approximants, following the general idea of the papers [
43,
46]. Such a hybrid approach, coined before Borel-light [
49], is particularly useful when the integral transformation cannot be applied because of singularities in the transformational
-functions or when an explicit formula and not only numerical values are required. In fact, we are dealing with the whole table of approximate values
with
[
49], while
for even
k, and
for odd
k. Formula (
34) is built only on the diagonal sequences in adherence with [
67].
Of course, when only the single approximant
of the order
q with
is employed, and only the highest possible order diagonal Padé approximant is considered, we return to the particular case of the scheme outlined in
Section 2.6. On the other hand, when
n is fixed and
i is varied in the Formula (
87), we are dealing with a different sequence of approximations from that described in
Section 2.6. The whole (almost) table (
87) was employed in the paper [
49].
Since, at
,
the sought amplitude is approximated by the table of values expressed by means of the hybrid formula
The parameter
is found from the optimization conditions. Such conditions are analogous to the Equations (
74)–(
76) [
49] but have to be applied properly to the marginal amplitudes
.
5. Examples
In the following sections, we are primarily concerned with the comparison of different factor-based approximations with root-based [
49] and Padé-based approaches [
58,
68].
5.1. Quantum Quartic Oscillator
For the quantum anharmonic oscillator [
69], perturbation theory yields a rather long expansion for the ground-state energy
Here, the parameter
measures a deviation of the anharmonic potential from the quantum harmonic oscillations. The coefficients
are rapidly growing in magnitude. Their concrete values are known for very high orders and can be retrieved from [
69]. The strong-coupling limit for
,
is also known, and the parameters are
and
.
Factor approximants, when applied to the truncation
, give the following sequence of approximate amplitudes:
Borel factor approximants with
and
give
Complex values are omitted here and in what follows.
The results produced by the following Borel-light approximations with a fixed marginal amplitude and varying correcting terms,
are slightly worse.
The best results are found for rather high orders with the following Borel-light approximations with varying marginal amplitudes and correcting terms of the same order,
Bear in mind that, for the factor approximants, the optimization problems can be solved analytically only for low orders. Only the minimal derivative problem possesses a unique solution. The results of optimization, such as the optimal control parameters, can be used as inputs to construct the sequence of Borel-light approximations. But, even in these cases, the results appear not to be better than those given by the factor approximants by themselves. However, such a strategy can be useful for some other problems and is exploited below.
Various other methods are considered as well. They give even more reasonable results, as shown in
Table 1.
The best result is marked in bold in
Table 1.
It is achieved in the ninth order of perturbation theory with the Fractional Borel Summation with iterated roots [
49]. An even better result,
, is found with the same method in the 10th order [
49].
The method of Corrected Padé approximants needs more terms to obtain the same accuracy [
70] but is able to easily scan very high orders [
70]. But, already in the 10th order, it gives a reasonable estimate,
.
5.2. Schwinger Model: Energy Gap
Let us consider the energy gap between the lowest and second excited states of the scalar boson for the massive Schwinger model in Hamiltonian lattice theory [
71,
72].
The energy gap
between the two states at small
z values can be represented (in low orders), according to the paper [
72], as
with the variable
. Here,
g stands for a coupling parameter and
a is the lattice spacing. The coefficients
are known up to the 13th order and can be found in [
72].
In the continuous limit,
follows the power law [
72]
where
and
.
Factor approximants, when applied to the original expansion, give the following sequence:
with a reasonable value given by
. The complex values are omitted again here and below.
Borel factor approximants, calculated with
and
, give the following sequence:
with a reasonable result that is very similar to the previous result
.
The most consistent results are found by applying the Borel-light approximations with a fixed marginal amplitude and varying correcting terms,
The last two terms can be considered lower and upper bounds, respectively.
As mentioned previously, the optimization problem with factor approximants can be solved analytically for low orders only. The optimization results, such as those for the optimal control parameters, can be used to construct the sequence of Borel-light approximations. Such a strategy appears to be useful, as shown below. Indeed, the application of fractional Borel
u-optimization with factor approximants for
amounts to solving the equation
It brings a sensible result for the sought amplitude
with the control parameter
.
The Borel-light approximations with the same marginal amplitude and varying correcting terms give even better results:
Different resummation methods give the results presented in
Table 2. The method of Fractional Borel-light summation with iterated roots [
49] also gives a rather good result:
. It is rather close to the value
obtained with finite-lattice calculations. Bear in mind that different advanced series methods give
[
72].
The best result is marked in bold in
Table 2.
5.3. Schwinger Model: Critical Amplitude
The ground-state energy
E of the Schwinger model is given by the very short truncated series [
72,
73,
74,
75,
76,
77,
78],
The large-
x behavior is also known. It is expressed in the form shown in (
4),
with
,
. Sometimes, an addition of one more trial term with
may help to improve the results.
Because of the negative
, we work with the inverse truncations when using
u-optimization. The problem appears to be quite difficult for the factor approximations, i.e., the factor approximant of the second order gives only
. The application of fractional Borel
u-optimization with the factor approximants for
,
amounts to solving the equation
It brings, after the inversion, a very good result for the sought amplitude
with the uniquely determined control parameter
. The result is the best among all of the results represented in
Table 3. It is even better than the result found from the self-similar root approximants, which explicitly employs the known subcritical index [
79].
A
u-optimal Borel-light technique with factors applied to the inverse quantities gives the inverse amplitude
and the total critical amplitude
B after inversion equals
.
For
p-optimization with
, the parameter
is found from the minimal derivative condition as the unique solution to the equation
with
and
.
As expected, the results found with optimizations and presented above represent an improvement over the non-optimized Borel factor techniques corresponding to and , which also produce the very reasonable .
The results of calculations by various methods are shown in
Table 3.
The best result is marked in bold in
Table 3.
It is worth stressing that it is always useful to attack the problem using several methods.
5.4. Anomalous Dimension
Consider the cusp anomalous dimension
of a light-like Wilson loop in the
supersymmetric Yang–Mills theory [
79,
80]. It depends only on the variable
expressed though the coupling
g. The problem can be written down in terms of the function
The latter function has the following weak-coupling expansion,
In the strong-coupling limit,
takes the form shown in (
4),
with
and
.
Direct application of the factor approximants in the second order gives , while in the third order, it gives . The Borel summation with factor approximants in the second order gives and is better than the direct factor approximation. In the third order, the Borel summation with factors gives complex results.
Neither of the optimization types brings a unique solution. The best result,
, is found with
p-optimization with
,
with the optimum obtained at
.
The Borel-light technique, when applied to inverse quantities, gives the inverse amplitude
and the total critical amplitude after inversion is
.
The result
is found for the
u-optimization
with the optimal value
.
The fractional Borel technique with iterated roots [
49] gives a unique solution
in the case of the
u-optimization problem.
Some results obtained by different methods are shown in
Table 4. The best result,
, is obtained from the methodology described above in
Section 3.2.
The best result is marked in bold in
Table 4.
5.5. Two-Dimensional Polymer: Swelling
For a two-dimensional polymer one can theoretically study the so-called swelling factor Y [
81]. For Y, perturbation theory yields expansions in the powers of the dimensionless coupling parameter
g. The swelling factor can be represented as the following truncation,
As
, the swelling factor is expressed in the form of (
4), i.e.,
Here, the critical index
is exact [
82,
83]. As for the amplitude
B, one can only say that it is of the order of unity.
Indeed, factor approximants and Borel factor techniques all give results close to unity. For
p-optimization with
, the parameter
p can be found from the minimal derivative condition
with
and
. A much lower result,
, is found with
u-optimization,
with the unique optimum located at
. The Borel-light summation technique fails to produce a holomorphic diagonal Padé approximation.
The results are presented in
Table 5. Most of the methods give results that are close to the conjectured value of unity. Even Padé approximants, factor approximants, and Borel summation with factors all give values close to one.
The best result is marked in bold in
Table 5.
5.6. Three-Dimensional Polymer: Swelling
For a three-dimensional polymer, one can find the swelling factor
[
81,
84], in the form of a truncated series of the type (
3),
The expansion (
94) can be extended to the sixth order [
81,
84]. The strong-coupling behavior of the expansion factor is expressed in the form of (
4),
with
, and
.
Factor approximants give the following sequence:
Borel factor approximants give the sequence
For p-optimization with , the parameter can be found from the condition The unique solution to the latter equation, , is found and the critical amplitude is .
The Borel-light technique gives
A much lower result,
, is found with
u-optimization,
with the optimal value
. However, the Borel-light technique gives
and the latter estimate for the critical amplitude is rather good.
The best result is marked in bold in
Table 6.
The fractional
p-optimal Borel technique with roots [
49] gives a fairly reasonable value,
. Consistent results have also been found with the Padé–Borel techniques shown in the paper [
58].
5.7. One-Dimensional Quantum Nonlinear Model
The ground-state energy of the Bose-condensed atoms in a harmonic trap can be expressed in terms of the function
, which can be expanded with a dimensionless coupling parameter
g. For the fifth order,
where all
values can be found in the paper [
85]. The coefficients rapidly decay by the absolute value. The strong-coupling limit is given in the form (
4),
Factor approximants give the following, apparently convergent, sequence for the amplitude:
Borel factor approximants give the sequence
which defines the upper and lower bounds.
For
p-optimization with
, the parameter
can be found as the unique solution to the equation
with
and
. The Borel-light technique gives
leading to fairly reasonable numbers,
, for the the critical amplitude.
A lower value,
, is found with
u-optimization from the equation
with the optimal value
. The Borel-light technique gives
and the critical amplitude,
is very good.
The results are presented in
Table 7. Almost all of them have a good level of accuracy.
The best result is marked in bold in
Table 7.
Particularly good results are obtained with the Fractional Borel-light optimal techniques. The Padé techniques described in the paper [
58] also give good results, improving the techniques of iterated roots.
5.8. Three-Dimensional Harmonic Trap
The ground-state energy
E of the trapped Bose condensate in the three-dimensional case was investigated in the paper [
86]. The energy can be approximated by the following truncation
with the “trapping” parameter
c.
For a very strong parameter
c, the energy behaves as the power law
And, the critical parameters
,
[
86] are known.
Factor approximation in the second order fails, while in the higher-order factor, the approximants give reasonably good results. However, the Borel summation with the factor approximants in the second order works very well, with .
The use of fractional Borel
u-optimization with the factor approximants for
,
amounts to solving the equation
It gives an amplitude of
with the control parameter
. The corrected Borel-light technique gives the amplitude
and the total critical amplitude
is the best among all estimates.
The result,
, is found with
p-optimization with
, which amounts to solving the equation
The optimum is found for
. The Borel-light technique gives the amplitude
and the total critical amplitude,
, is quite reasonable. The results from the calculations using different methods are presented in
Table 8. All of them give close and rather accurate values.
The best result is marked in bold in
Table 8.
In both cases of trapped Bose condensates discussed above, the methods based on the idea of corrected Borel-factor-light approximants work well. They give numbers that are better than those of the Padé techniques or those of methodologies based on iterated roots.
5.9. Bose Temperature Shift
The ideal Bose gas is unstable below the condensation temperature
[
86]. Atomic interactions induce the shift
to the realistic
of a non-ideal Bose system. The shift is characterized by the ratio
for the asymptotically small gas parameter
.
Monte Carlo simulations [
87,
88,
89] suggest that
In order to calculate
theoretically, it has been suggested that one should first calculate an auxiliary function
[
90,
91,
92]. Then, one can find
as follows:
The latter limit is found from the expansion over an effective coupling parameter,
In the fourth order of factor approximants, we find , which is much smaller than that expected from the simulations.
In what follows, we work with the original
. For instance, in
p-optimization with
, we can set, by analogy with Equation (
78),
The latter equation gives the optimal solution
. Using the Borel-light technique with the optimal
, we have a good result
Even by just setting
,
, we arrive at
The application of fractional Borel
u-optimization with factor approximants amounts to solving the equation
which is written in analogy with (
75). It gives the control parameter
. The correction with the Borel-light technique gives the amplitude
where
, and the sought amplitude is
.
Note that there is another solution to the equation (
101),
, and it gives
. The latter result appears to be identical to the results of Modified Even Padé summation [
58]. The two solutions are very similar to each other, and the non-uniqueness in such a case does not pose a serious problem.
Various results are shown in
Table 9.
The best result is marked in bold in
Table 9.
In the case of the Bose temperature shift, the Borel-light method, based on the idea of corrected Borel-type approximants, optimization, and correction with Padé approximations works well. The Padé method modified for an even number of terms in the expansion gives a result that is close to the latter, even without optimization or Borel transformation. But, it should be noted here that the even approximation is close in spirit to the general idea of the corrected approximants, as described in the
Section 2.6 and in the paper [
70].
Accurate results for the sought parameter are also found with optimal generalized Borel summation with iterated roots,
[
44]; with the optimal Mittag–Leffler summation with iterated roots with
[
46]; and with the corrected iterated roots,
(see [
49], and references therein). Kastening [
90,
91,
92], using the optimized variational perturbation theory, estimated
as
.
5.10. Fermi Gas: Unitary Limit
The ground-state energy
E of a dilute Fermi gas can be obtained from perturbation theory [
93], so that
with the coefficients
The effective coupling parameter
is simply related to the Fermi wave number
and the atomic scattering length
. The limit of very large effective coupling
g is called the unitary Fermi gas [
94]. Monte Carlo simulations for the case of
[
95,
96] yield
The experimental value [
95,
97]
is rather close to the Monte Carlo results.
Factor approximants in the third and fourth orders give rather high estimates for
, as shown in
Table 10, while Borel summation with factor approximants gives rather low estimates.
The best result is marked in bold in
Table 10.
None of the methods, including those using optimizations, bring accurate solutions.
u-Optimization gives
with or without correction terms. The result
is found with
p-optimization by means of the equation
with the optimum at
. However, the Borel-light techniques give a holomorphic correcting approximant with the value of the amplitude in the sixth order being
The corresponding approximant
that leads to a very accurate estimate is shown below:
The lower-order approximants, appear to be non-holomorphic.
From
p-optimization with the minimal derivative condition
we find
. The Borel-light technique then gives a reasonable estimate for the amplitude
We conclude that the Borel-light methods based on Borel summation with optimization and subsequent correction with Padé approximations work well for the problem of a unitary limit, approaching the quality of Monte Carlo and experimental results.
6. In Lieu of Conclusions
For the first time, Fractional Borel Summation was applied in conjunction with self-similar factor approximants. This is the main distinction from the method developed in [
49] in conjunction with the so-called self-similar iterated root approximants. It was found that the technique of Fractional Borel Summation can be most conveniently applied in hybrid form. Such hybrid approximants emerge when the Borel-transformed factor approximations are complemented multiplicatively with Padé approximants to satisfy the original expansions asymptotically. A detailed comparison of different methods is performed on a large set of examples, including the approximation techniques of self-similarly modified Padé-Borel approximations. Such a comparison clearly emphasizes the strong points of each of the techniques.
We confirm that the quality of the analytical reconstructions can match the quality of heavy numerical work. Analytical results are often very similar to the exact numerical data. They follow after a sequence of a few analytical steps, and the most difficult part of solving transcendental equations is performed numerically with any desired accuracy.
The discussed approach to the resummation of the asymptotic series is multi-leveled. General Borel Fractional transformation of the original series was introduced. The found transformed series should be resummed in order to adhere to the asymptotic power laws. One starts with the formulation of dynamics in the approximations space. To this end, self-similarity is employed. The flow in the approximation space should be controlled, and “deep” control is incorporated into the definitions of the self-similar approximants. Certain classes of self-similar approximants satisfying the power law behavior by design, such as root and factor self-similar approximants, are chosen for reasons of transparency, explicitness, and convenience. We also employ properly modified (taking account of power law behavior) Padé approximants by noting that they can be viewed as a particular case of factor approximants. The asymptotic power law properties of the approximations follow directly from the ways in which controls are introduced into the approximation dynamics. The second level of controls concerns the dynamics of Borel-type approximations that emerge after an inverse transformation is accomplished back to the original space of approximations. Self-similarity of the approximants allows us to find the dependencies of the sought critical properties explicitly. The second level of control is then performed by applying the minimal difference and minimal derivative conditions with respect to the parameters that are explicit in the original Borel Fractional transformation.
The main methods compared in the paper each have their own merits. Different methods, based on Padé approximants [
58], can be useful as benchmarks for the evaluation of the results. The standard scheme of odd Padé approximants is not competitive with respect to the other considered methods. Of course, this statement concerns only the physical problems considered above. However, the methods based on Padé approximants are indispensable for computations with very long expansions. For shorter expansions, the other methods discussed in the paper should be used. But, in many cases of such expansions, the diagonal Padé approximants reappear in combination with factors and roots. Such hybrid forms of corrected approximants, when the approximations of different types are applied consequentially, are very useful for the practical purpose of accurate summation.
In conclusion, we recapitulate the main steps of the multi-level methodology developed in the present paper.
The initial truncated series (
2) is transformed into the form of a new, supposedly better behaving series. The chosen transformation is the Fractional Iterated Transform (see
Section 3.6). At this stage, the control parameters are introduced. They have to be found at the final optimization step.
To arrive to the correct asymptotic behavior at infinity, the transformed series is approximated by the class of approximants with power law behavior (
4) at infinity. At this step, we reconstruct the coefficients for an arbitrary
n.
When the guiding principle of self-similarity, described in Section refself, is applied together with the so-called algebraic transformation of the original series [
32], we arrive at two types of approximants with the desired property at infinity. These are the self-similar roots discussed in
Section 2.4 and the self-similar factor approximants described in the
Section 2.5 and
Section 4.1.
In the limit of a large x, the expressions for the critical amplitudes become explicit and factorize into the parts emerging from the self-similar approximants applied to the transformed series and from the -functional terms emerging from the inverse transformation, being dependent on the type of transformation.
The application of self-similar iterated roots was considered previously [
49]. Now, after the optimization with the minimal difference and (or) minimal derivative conditions described in
Section 2.2 and
Section 3.6, we resort to the best known (to us) guiding principle for
numerical convergence of the sequences of the approximations for the sought quantities. The uniqueness of the solution for the given class of approximants is achieved due to careful selection of the transformation to the original truncated series.
In the current paper, we employ the self-similar factor approximants together with the most convenient technically and theoretically sound method of so-called hybrid techniques (see
Section 2.6 and
Section 4.3). Typically, in low orders of the perturbation theory, an optimized factor approximant has to be found. To return to the original series, we restore the factors in the form of the diagonal Padé approximants. Besides numerical convergence, we are also guided by the Gonchar results on the convergence of the diagonal Padé approximants [
67]. Such a selection leads, by design, to a unique limit that can be found numerically as the approximation for the sought quantity.
The Fractional Borel Summation with iterated roots can be successfully applied for various problems [
49], although it is not always the best method. Iterated roots can be exceptionally helpful for problems which appear to be indeterminate and those that are poorly treatable from the standpoint of the Padé approximations [
49] and the factor approximations. The cases with fast-growing and rapidly diminishing by magnitude coefficients
are better treated by means of Borel summations with iterated roots. The cases with slowly diminishing by magnitude coefficients
are best approximated by Borel summations with factors. Fractional Borel Summation with factor approximants in its different realizations can be successfully applied to various problems of the types discussed above. In addition, it is shown to be the best method, or very close to it, for about half of the problems considered in the paper, much more often than the other methods discussed above. Thus, the techniques based on factors are most useful in the context of finding the best method for the particular problems. One can say that the factor approximants are a sharper tool than the iterated roots, but the iterated roots are useful for a wider range of applications.