Optimized Self-Similar Borel Summation

The method of Fractional Borel Summation is suggested in conjunction with self-similar factor approximants. The method used for extrapolating asymptotic expansions at small variables to large variables, including the variables tending to infinity, is described. The method is based on the combination of optimized perturbation theory, self-similar approximation theory, and Borel-type transformations. General Borel Fractional transformation of the original series is employed. The transformed series is resummed in order to adhere to the asymptotic power laws. The starting point is the formulation of dynamics in the approximations space by employing the notion of self-similarity. The flow in the approximation space is controlled, and ``deep'' control is incorporated into the definitions of the self-similar approximants. The class of self-similar approximations, satisfying, by design, the power law behavior, such as the use of self-similar factor approximants, is chosen for the reasons of transparency, explicitness, and convenience. A detailed comparison of different methods is performed on a rather large set of examples, employing self-similar factor approximants, self-similar iterated root approximants, as well as the approximation technique of self-similarly modified Pade - Borel approximations.


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.

Optimization and Self-Similarity
In this section, we present the main ideas of the approach in order for the reader to grasp the basis of the particular techniques that are used.The two pivotal points are the notions of convergence optimization and self-similarity.

Asymptotic Series
The starting point of the consideration concerns the well-known fact that, in practical applications, we often satisfy the necessity of solving problems by applying some kind of perturbation theory, resulting in expansions to powers of a parameter or a variable.Such expansions usually represent asymptotic series diverging for finite values of the expansion variable.
To be specific, let us be interested in a real function f (x) of a real positive variable x, defined by rather complicated equations that, because of their complexity, can be solved only by means of a kind of perturbation theory.The latter would result in an expansion in powers of asymptotically small x → 0: having the form of a truncated series where f 0 (x) is a known function.Our main concern is the summation of the power expansion, because of which we shall concentrate on the expansion of the form As is evident, the more general form (2) can be easily reduced to (3).
In practically all interesting problems, expression (3) represents an asymptotic series diverging for a finite value of x.At the same time, the quantities of interest usually correspond to finite and sometimes to quite large values of the variable x.Moreover, in many cases, the most interesting region is the region of very large x → ∞, where the sought function behaves according to the power law Thus, the principal problem is how it would be possible to reconstruct the large-variable behavior of the sought function (4), knowing only the small-variable expansion (3)?

Control Functions
It seems clear that, in order to define an effective sum of the divergent series (3), it is necessary to rearrange the divergent series into a convergent series.When one is interested in rather large values of x, then one cannot resort to Padé approximants, since for asymptotically large x values, the limit of a Padé approximant is proportional to x M −N , depending on which of the approximants from the Padé table is accepted.In that sense, the large-variable behavior of Padé approximants is not defined, since x M −N can tend to either +∞, −∞, or zero, depending on whether M > N, M < N, or M = N.This implies that Padé approximants do not converge for large x → ∞.
To force a sequence to converge, it is necessary to introduce some quantities that regulate the convergence.Such quantities, governing the sequence convergence, can be named control functions The following highlights are the same.The idea of introducing control functions so that the reorganized sequence converges was advanced in Refs.[15,16].Control functions are to be defined by optimization conditions regulating the sequence convergence, so that the initial perturbation theory is reorganized to an optimized perturbation theory yielding a converging sequence of optimized approximants [15,16].All mathematical details can be found in recent reviews [17,18].
The introduction of control functions u k transforms the initial sequence {f k (x)} into a sequence {F k (x, u k )}.The latter can be symbolically denoted as ( This introduction can be conducted in several ways: by incorporating control functions into an initial approximation of perturbation theory, through a change in variables, or by a direct sequence transformation [17,18].The derivation of equations defining control functions starts with the Cauchy criterion of convergence.The criterion tells that a sequence converges if and only if, for any ε > 0, there exists a number k ε such that for all k > k ε and p > 0. From the Cauchy criterion, one can derive (see details in [17,18]) the optimization conditions defining the control functions.One possibility is the minimal difference condition whose simplest variant is The other possibility is the minimal derivative condition Several other representations of optimization conditions are admissible [17,18].In some cases, control functions can be defined by the boundary conditions implying asymptotic equivalence at a small variable of the renormalized and initial terms of the approximations.The latter condition is also called the accuracy-through-order procedure.
After the control functions u k (x) are found from the optimization conditions, the optimized approximants read as

Self-Similar Relation
Self-similar relations are known to arise in renormalization group theory, where one considers scaling with respect to spatial or momentum variables, as in the real-space decimation procedure [19,20] or in quantum field theory [21,22].Then, self-similar relations connect characteristic quantities, such as effective Hamiltonians, Lagrangians, or correlation functions, with different spatial or momentum scales.In that sense, renormalization groups in statistical physics or quantum field theory provide relations between the characteristic quantities with scaled variables.The notion of self-similarity in approximation theory, introduced in Refs.[23,24], does not concern a scaling of variables, but rather, the scaling of approximation orders.The number labeling the approximation order plays the role of discrete time.
It seems natural then to construct a dynamical system in the space of approximations.To this end, it is necessary to introduce an endomorphism into the approximation space.Let us start with an initial approximation defining an expansion function x k (f ).Then, it is possible to define the approximation function The approximation space is given by the closed linear envelope over all admissible approximation functions: This is similar to approximation spaces in approximation theory [25] or to phase spaces in physics [26][27][28].Thus, we obtain a dynamical system in discrete time, called the approximation cascade: The points of the cascade form the cascade trajectory.The Cauchy criterion (6) now acquires the form If control functions have been chosen so that the sequence of the optimized approximants is convergent, then, in terms of dynamical theory, this implies the existence of a fixed point, where Conditions ( 16) and ( 17) lead [17,18,23,24] to the self-similar relation Fixed-point solutions y * k (f ) to this equation, using the relation (13), give the self-similar approximants ).The stability of the method can be checked by investigating map multipliers, similarly to the stability analysis conducted in numerical calculations [29,30].
The self-similar relation (18) describes the motion of a dynamical system (cascade).Generally, a dynamical system does not necessarily tend to a fixed point, but it can display chaotic motion [31].This is why the incorporation into the approximation cascade of control functions, governing sequence convergence, is so important.

Self-Similar Root Approximants
Generally, the self-similar relation, depending on additional imposed constraints playing the role of boundary conditions, can lead to different solutions.For instance, when the fixed-point solution is required to satisfy the prescribed asymptotic expansions at small as well as at large variables, we obtain [32] (see also [33,34]) the self-similar root approximants When the large-variable behavior is of the law with the known power β, then we have which results in the approximant All control parameters A j are uniquely defined by the asymptotic boundary condition The root approximant (22) at large x values behaves as which gives the amplitude Below, we refer to the root approximants introduced above as iterated root approximants or just as iterated roots.

Self-Similar Factor Approximants
Accepting the asymptotic expansion (3) as a boundary condition and representing this expansion in the form yields [35,36] the fixed-point solution in the form of self-similar factor approximants where For odd k values, scaling relations [17] allow us to set A 1 = 1.All other control parameters A j and n j are uniquely defined by the boundary condition requiring the asymptotic expansion of f * k (x) at small x values to coincide with f k (x).The factor approximants serve as a generalized representation of the products of functions [37,38].It has been shown that the self-similar factor approximants provide good accuracy for a wide class of problems [35,36,39,40].
At large x values, the approximant (26) gives with the amplitude and the large-variable exponent It is appropriate to notice that Padé approximants are a particular kind of factor approximants, where a portion of the powers equals +1 and the other portion is −1, so that = P M/N (x).

Self-Similar Combined Approximants
It is possible to combine different types of approximants.Consider the situation when one type of approximation better suits the small-variable limit, while the other type of approximation is better in the description of the large-variable limit.For this purpose, one can consider the first q terms of the expansion f k (x) with q < k, When constructing a self-similar approximant f * q (x) on the basis of expansion (31), we define the correcting ratio and expand it in powers of x to obtain On the basis of the latter expansion, we construct a self-similar approximant C * k/q (x).The final combined approximant is As has been mentioned above, Padé approximants are a particular kind of self-similar approximant.So, for the correcting function C * k/q (x), it is possible to take a Padé approximant P M/N (x) with M + N = k − q.In that case, the final approximant is Of course, the necessary boundary conditions have to be preserved when choosing q, M, and N.
3 Self-Similar Borel-Type Transforms The convergence of series is known to be improved by resorting to Borel summation [2,41,42].Borel, or Borel-type summation, can be combined by employing self-similar approximants [43].

Self-Similar Borel Transform
The Borel transform for expansion (3) is with the inverse transformation The series (36) can be summed by employing self-similar approximants, obtaining a selfsimilar Borel transform B * k (x).Then, the inverse transformation (37) gives The most difficult problem is the study of the large-variable limit.Therefore, we pay more attention to this limiting behavior.In the present case, by substituting the large-variable form of the self-similar approximant, into the integral (37), we obtain where the large-variable amplitude is

Self-Similar Borel-Leroy Transform
The Borel-Leroy transform reads as where λ is a control parameter that has to be defined from optimization conditions, providing convergence to the sequence of approximants.The inverse transformation is Constructing a self-similar approximation B * k (x, λ k ) yields the self-similar optimized Borel-Leroy approximant where the control parameter λ = λ k can be defined, e.g., from either the minimal difference optimization condition (8) or the minimal derivative optimization condition ( 9) for the largevariable amplitude.In the large-variable limit, the form results in the approximant (40) with the amplitude

Self-Similar Iterated Borel Transform
The method used for the Borel and Borel-Leroy transforms can be repeated several times, resulting in iterated transforms [9].For example, after accomplishing the first step and obtaining the transform it is straightforward to repeat the procedure to obtain the double Borel-Leroy transform After p iterations, one obtains For simplicity, one assumes that the control parameters at different steps are the same, so that λ (j) = λ, which gives The inverse transformation is By constructing a self-similar approximant B * kp (x, λ kp ) on the basis of the p-iterated transform (50), one obtains [44] the self-similar iterated approximants where the control parameters λ kp are defined from the optimization conditions.For the largevariable limit of the transform the large-variable behavior of the sought function becomes with the amplitude

Self-Similar Mittag-Leffler Transform
A generalization of the Borel summation can be conducted by means of the Mittag-Leffler [45] summation.By introducing the Mittag-Leffler transform the inverse transformation reads as The series (56) can be represented by a self-similar approximant, giving M * k (x, µ k ), where the control parameter µ k is found from the optimization conditions [46].Then, the inverse transformation yields the self-similar approximant For large variables, the transform leads to large-variable behavior of the self-similar approximant with the amplitude For the Mittag-Leffler transform, one can define iterated transforms, similarly to the iterated Borel-Leroy transforms.

Self-Similar Modified Borel-Leroy Transform
The Borel-Leroy transform can be modified [47,48] to the form with the inverse transformation being in which ν is an integer, so that Constructing a self-similar transform B * k (x, λ k , ν), where λ k is a control parameter, gives [49] the inverse transformation With the large-variable behavior of the transform the large-variable limit of the self-similar approximant (65) reads as shown in Equation ( 60) with the amplitude As is evident, setting λ k = 0 leads to the modified Borel transform.

Self-Similar Fractional Iterated Transform
The modified Borel-Leroy transform shown in the previous section can be further generalized in two aspects: First, it is straightforward to iterate it.And, second, it is possible to treat the parameters ν as well as the iteration order p, not as integers, but as fractional quantities [44,49].
Starting with the iterated modified Borel-Leroy transform we have the inverse transformation Then, by constructing a self-similar approximant B * kp (x, λ kp , u kp ), where λ kp and u kp and p = p k are control parameters, we obtain a self-similar fractional iterated transform Strictly speaking, the differential operator in action ( 64) is defined only for the integer ν.However, it is admissible to formally treat ν as an arbitrary real number u when it acts on power law expressions, considering Equation ( 64) as a definition.
Then, for the large-variable behavior of the transform we find the large-variable limit for the sought function in which the amplitude is In this, and in the previous sections, the control parameters are defined by the optimization conditions, as discussed in Section 2.2, as the minimal difference or minimal derivative conditions.In particular, when one is interested in the correct evaluation of the large-variable amplitudes, these conditions should be used with respect to the amplitudes B kp (u) considered as functions of the control parameters u.For brevity, here, we use a single control parameter u, although there can be several of them.It is possible to use the minimal difference condition in the form or in the form Otherwise, one can resort to the minimal derivative condition For both conditions, the minimal difference and minimal derivative are equivalent.Previously, in the paper [44], it was suggested that the number of iterations p should be considered a continuous control parameter.The multidimensional integrals are relatively easy to define for the integer p. Introducing a continuous p implies a smooth interpolation between the values of the integrals for discrete p values.The approach is applicable only in the limiting case of large x values.
While setting u = 0, the parameter p = p k can be found from the minimal derivative condition as the unique solution to the equation Alternatively, we can solve the minimal difference equation with a positive integer k = 1, 2, 3, . ... For the special singular cases of β = −1, −2 . .., one can study the inverse quantity.One of the main advantages of the method involving the combination of the methodology shown in Section 2.6 and the Borel-light summation shown in Section 4.3 is that they do not involve any explicit singular terms.This property allows one to avoid an extra transformation and allows one to work with the original truncations directly.In such cases, one is bound to the optimizations with marginal amplitudes, since the complete expression for the amplitudes involves singularities.

Critical Amplitudes from Fractional Iterations
Consider the case where one has to find explicitly a real, sign-definite, positive-valued function f (x) of a real variable x.In addition, the function possesses the power law asymptotic behavior (4).The critical exponent β > 0 is known.The case of a negative β will be considered by studying the inverse of f (x).Let us find the amplitude B.
We consider the case where it is impossible to find the sought function f (x) 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).
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 λ kp 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.
1.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].
2. 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.
3. The third method, Fractional Borel-light or self-similar combined approximants is explained in great detail in Sections 2.6 and 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 B kp (or with the marginal amplitudes C kp ), 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 Sections 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 β = 0 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.

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 n j 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 A j 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.

Modified Padé Approximations
Alternatively, in place of the factors, one can apply the well-known Padé approximants P n/m [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 k = 1, 3 . .., one should study the following transformed series for the sought function f (x), T (x) = f (x) −1/β .The following modified Padé approximants are defined for even cases, where n = 0, 1 . . ., n max is a non-negative integer.
In the odd case, (2n + 1) max = k.The approximants (82) do comply with the power law (4) at x → ∞.One can relatively easily find the sequence of approximations for the amplitudes.In the case of even k = 2, 4, . .., the following modified-even Padé approximants are defined [58].Here, the corrector was introduced to ensure the correct asymptotic behavior [58].Also, G(x) = f (x) K(x) 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, n = 1, 2 . . ., n max , is a positive integer.In the even case, n max = k/2.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.

Critical Amplitudes from Hybrids of Factors and Padé Approximations
The sought function can be reconstructed from the transform f * i (x) 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 i = 1, 2 . . .k [49], while n = 1, 2 . . .k/2 for even k, and n = 1, 2 . . .(k − 1)/2 for odd k.Formula ( 34) is built only on the diagonal sequences in adherence with [67].
Of course, when only the single approximant f * q (x) of the order q with q < k 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 x → ∞, f * i (x, u) ≃ C i (u)x β , the sought amplitude is approximated by the table of values expressed by means of the hybrid formula The parameter u = u i 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 C i (u i ).

Examples
In the following sections, we are primarily concerned with the comparison of different factorbased approximations with root-based [49] and Padé-based approaches [58,68].

Quantum Quartic Oscillator
For the quantum anharmonic oscillator [69], perturbation theory yields a rather long expansion for the ground-state energy Here, the parameter g ≥ 0 measures a deviation of the anharmonic potential from the quantum harmonic oscillations.The coefficients a n 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 g → ∞, is also known, and the parameters are B = 0.667986 and β = 1/3.Factor approximants, when applied to the truncation E k (g), give the following sequence of approximate amplitudes: 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, (1) = P 13/13 (∞) C 9 (1) = 0.67397.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, B = 0.669356, 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, B = 0.655086.

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 ∆(z) between the two states at small z values can be represented (in low orders), according to the paper [72], as 2∆(z) ≃ 1 + 6z − 26z 2 + 190.6666666667z 3 − 1756.666666667z 4 + 18048.33650794z 5 , (90) with the variable z = (1/ga) 4 .Here, g stands for a coupling parameter and a is the lattice spacing.The coefficients a n are known up to the 13th order and can be found in [72].
In the continuous limit, ∆(z) follows the power law [72] ∆ where B = 1.1284 and β = 1/4.Factor approximants, when applied to the original expansion, give the following sequence: with a reasonable value given by B 8 .The complex values are omitted again here and below.Borel factor approximants, calculated with u = 0 and p = 1, give the following sequence: 05511, with a reasonable result that is very similar to the previous result B 8 .
The most consistent results are found by applying the Borel-light approximations with a fixed marginal amplitude and varying correcting terms, 15488.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 p = 1, 2 amounts to solving the equation It brings a sensible result for the sought amplitude B = 1.21118 with the control parameter u = u 3 = 0.286426.
The Borel-light approximations with the same marginal amplitude and varying correcting terms give even better results: 18951.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: B = B 6,10 = 1.1452.It is rather close to the value B = 1.14(3) obtained with finite-lattice calculations.Bear in mind that different advanced series methods give B = 1.25 (15) [72].The best result is marked in bold in Table 2.

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 B = 0.6418, β = −1/3.Sometimes, an addition of one more trial term with a 3 = 0 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 B ≈ 0.5456.The application of fractional Borel u-optimization with the factor approximants for k = 3, p = 1, 2 amounts to solving the equation It brings, after the inversion, a very good result for the sought amplitude with the uniquely determined control parameter u = u 3 = 0.0961685.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 0.5827.
For p-optimization with u = 0, the parameter p = p 3 is found from the minimal derivative condition as the unique solution to the equation As expected, the results found with optimizations and presented above represent an improvement over the non-optimized Borel factor techniques corresponding to u = 0 and p = 1, which also produce the very reasonable B ≈ 0.6351.
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.

Anomalous Dimension
Consider the cusp anomalous dimension Ω(g) of a light-like Wilson loop in the n = 4 supersymmetric Yang-Mills theory [79,80].It depends only on the variable x = g 2 expressed though the coupling g.The problem can be written down in terms of the function f (x) = Ω(x) x .The latter function has the following weak-coupling expansion, f (x) ≃ 4 − 13.1595x + 95.2444x 2 − 937.431x 3 , x → 0. In the strong-coupling limit, f (x) takes the form shown in (4), Direct application of the factor approximants in the second order gives B ≈ 2.1307, while in the third order, it gives B ≈ 1.8389.The Borel summation with factor approximants in the second order gives B ≈ 1.8798 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, B ≈ 2.0233, is found with p-optimization with u = 0, with the optimum obtained at p = p 2 ≈ 0.43297.
The Borel-light technique, when applied to inverse quantities, gives the inverse amplitude and the total critical amplitude after inversion is B = 1.92018.meaning is retained.The result B ≈ 2.1115 is found for the u-optimization with the optimal value u = u 2 ≈ 0.263034.The fractional Borel technique with iterated roots [49] gives a unique solution B = 1.90291 in the case of the u-optimization problem.
Some results obtained by different methods are shown in Table 4.The best result, B = 2.0118, is obtained from the methodology described above in Section 3.2.
The best result is marked in bold in Table 4.

Two-Dimensional Polymer: Swelling
For a two-dimensional polymer one can theoretically study the so-called swelling factor Υ [81].For Υ, perturbation theory yields expansions in the powers of the dimensionless coupling parameter g.The swelling factor can be represented as the following truncation, As g → ∞, the swelling factor is expressed in the form of (4), i.e., Υ(g) ≃ Bg β .
Here, the critical index β = 1/2 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 u = 0, the parameter p can be found from the minimal derivative condition  The best result is marked in bold in Table 5.

Three-Dimensional Polymer: Swelling
For a three-dimensional polymer, one can find the swelling factor Υ(g) [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 B ≈ 1.531, and β ≈ 0.3544.Factor approximants give the following sequence: Borel factor approximants give the sequence For p-optimization with u = 0, the parameter p = p 3 can be found from the condition ∂B 3 (p) ∂p = 0.The unique solution to the latter equation, p = p 3 = 0.358042, is found and the critical amplitude is B = B 3 = 1.53441.
The Borel-light technique gives with the optimal value u = u 2 = 0.369247.However, the Borel-light technique gives 53475, and the latter estimate for the critical amplitude is rather good.
The results are shown in Table 6.The best result is marked in bold in Table 6.
The fractional p-optimal Borel technique with roots [49] gives a fairly reasonable value, B = 1.53523.Consistent results have also been found with the Padé-Borel techniques shown in the paper [58].

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 f (g), which can be expanded with a dimensionless coupling parameter g.For the fifth order, f 5 (g) = 1 + 5 n=1 a n z n , where all a n 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), f (g) ≃ 3  2 g 2/3 .Factor approximants give the following, apparently convergent, sequence for the amplitude: 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.

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 B = 5 4 , β = 2/5 [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 B ≈ 1.25983.
The use of fractional Borel u-optimization with the factor approximants for k = 3, p = 1, 2 amounts to solving the equation The optimum is found for p = p 2 = 0.705782.The Borel-light technique gives the amplitude and the total critical amplitude, B = 1.28598, 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.

Bose Temperature Shift
The ideal Bose gas is unstable below the condensation temperature T 0 [86].Atomic interactions induce the shift ∆T c ≡ T c − T 0 to the realistic T c of a non-ideal Bose system.The shift is characterized by the ratio ∆T c T 0 ≃ c 1 γ , for the asymptotically small gas parameter γ → 0. Monte Carlo simulations [87][88][89] suggest that In order to calculate c 1 theoretically, it has been suggested that one should first calculate an auxiliary function c 1 (g) [90][91][92].Then, one can find c 1 as follows: The latter limit is found from the expansion over an effective coupling parameter, In the fourth order of factor approximants, we find c 1 ≈ 1.1, which is much smaller than that expected from the simulations.
In what follows, we work with the original c 1 (g).For instance, in p-optimization with u = 0, we can set, by analogy with Equation (78), The latter equation gives the optimal solution p = p 2 = 0.41657.Using the Borel-light technique with the optimal p 2 , we have a good result 28421.Even by just setting p = 1, u = 0, we arrive at 17351 .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 u = u 2 = −0.4351064.
The correction with the Borel-light technique gives the amplitude , where C 2 (u 2 ) ≡ C 2,1 (u 2 ), and the sought amplitude is B = B * 2,2 = 1.27224.Note that there is another solution to the equation (101), u 2 = −0.56557,and it gives B * 2,2 = 1.28553.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, c 1 = 1.339 [44]; with the optimal Mittag-Leffler summation with iterated roots with c 1 = 1.3397 [46]; and with the corrected iterated roots, c 1 = 1.3092 (see [49], and references therein).Kastening [90][91][92], using the optimized variational perturbation theory, estimated c 1 as 1.27 ± 0.11.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.
1.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.
2. 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.
3. 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 Sections 2.5 and 4.1.
4. 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.
5. 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 Sections 2.2 and 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.
6.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 Sections 2.6 and 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 a n are better treated by means of Borel summations with iterated roots.The cases with slowly diminishing by magnitude coefficients a n 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.

Table 1 :
Critical amplitude for the one-dimensional quartic oscillator

Table 3 :
Schwinger model: amplitude optimum located at u = u 3 ≈ 0.443811.The Borel-light summation technique fails to produce a holomorphic diagonal Padé approximation.The results are presented in Table5.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.

Table 7 :
Non-linear quantum model

Table 9 :
Shift of the Bose-Einstein condensation temperature and analogous models

Table 10 :
[49]i gas energy in the unitary limit 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.