Preprint
Article

This version is not peer-reviewed.

The Sequential Radii of Convergence for Power Series

Submitted:

01 June 2026

Posted:

03 June 2026

You are already at the latest version

Abstract
The classical radius of convergence provides a definitive geometric boundary for the validity of a power series, yet it offers little insight into the actual efficiency of the approximation at a specific point. This paper introduces the concept of sequential radii of convergence to bridge the gap between the qualitative "where" of convergence and the quantitative "how fast" of partial-sum approximation. By defining the indexed radius—which measures the number of terms required to meet a prescribed error tolerance—and the geometric radius, we provide a framework for evaluating the computational cost of series representations. We categorize power series into three fundamental tail-decay regimes: sublinear, linear, and superlinear. Within this taxonomy, we derive explicit asymptotic formulas that reveal how the speed of convergence is influenced by both the coefficient growth and the evaluation point’s proximity to the boundary of the convergence disk. Furthermore, we investigate the stability of these sequential radii under standard algebraic operations, including sums, products, and quotients. To demonstrate the practical utility of this refinement, we present a comprehensive catalog of classical series—ranging from the logarithmic and binomial series to more complex forms like the Lambert W0 and Jacobi theta-type series—analyzing their behavior through this new lens. Ultimately, this work offers a more granular perspective on power-series analysis, transforming the classical radius from a static boundary into a dynamic measure of analytic efficiency.
Keywords: 
;  ;  ;  

1. Introduction

1.1. The Classical Radius of Convergence

A power series centered at a R has the form
f ( x ) = n = 0 c n ( x a ) n .
Its classical radius of convergence R [ 0 , ] determines the interval | x a | < R on which the series converges absolutely, while the boundary points require separate analysis. Thus (R) gives the first geometric description of the region where the formal series represents a genuine function. Various mathematicians have contributed to its formal definition and computation as follows:
Brook Taylor (1685–1731) played a central role in the early development of power-series representation through Taylor’s theorem, which expresses a sufficiently regular function by polynomial approximations near a point. Although the modern theory of convergence came later, Taylor’s work established the local expansion viewpoint that remains fundamental in analysis [1,2,3].
Jean le Rond d’Alembert (1717–1783) refined the study of infinite series through the ratio test. In the setting of a power series, this test gives the familiar computation
R = lim n c n c n + 1 ,
whenever the limit exists. This provided a practical way to determine the interval of convergence directly from successive coefficients [4].
Augustin-Louis Cauchy (1789–1857) gave the theory of infinite series a rigorous foundation through systematic use of limits and convergence tests. His root-test viewpoint connects convergence to the exponential size of the coefficients and leads naturally to formulas involving | c n | n . This shifted the theory from formal manipulation of series toward asymptotic control of coefficient growth [5,6].
Niels Henrik Abel (1802–1829) clarified one of the most delicate issues in the theory of power series: the relation between convergence inside the interval and limiting behavior at the boundary. Abel’s theorem shows, in particular, that when a power series converges at an endpoint, its sum is approached by the values of the series from inside the interval [7]. This result gave a precise form to the boundary question that remains outside the direct information provided by the radius R itself.
Jacques Hadamard (1865–1963) completed the modern coefficient-growth formulation through what is now called the Cauchy–Hadamard formula:
R = 1 lim sup n | c n | n ,
with the conventions 1 / 0 = + and 1 / = 0 . This is the form most commonly used today when the ordinary ratio limit does not exist, since the lim sup captures the upper exponential growth rate of the coefficient sequence. In this sense, the classical radius is determined by the asymptotic size of the coefficients, and the Cauchy–Hadamard formula gives the standard bridge between coefficient growth and the geometry of convergence [5,8,9].

1.2. Motivation

The classical radius of convergence is one of the central organizing notions in the theory of power series. It gives a sharp geometric description of the interval, or disk in the complex setting, on which the series converges. In this sense, the classical radius answers the first natural question: where does the power series converge?
However, convergence is not only a question of location. Once a point (x) lies inside the interval of convergence, a second question becomes natural: how many terms are required before the partial sums approximate f ( x ) within a prescribed tolerance ε ? The classical radius R does not answer this question. It determines the domain of convergence, but it does not quantify the computational cost of convergence inside that domain.
This distinction is visible even among standard power series. Two series may have the same classical radius, while their tails decay at very different speeds. A geometric tail, an algebraically corrected geometric tail, and a factorial tail may all arise from familiar analytic functions, but they lead to different approximation costs. Thus, the usual classification by the classical radius alone does not fully describe the practical behavior of partial-sum approximation.
Standard references on analysis and power series usually emphasize convergence tests, endpoint behavior, uniform convergence on compact subintervals, and analytic continuation. These tools are fundamental, but they mainly describe qualitative convergence. They do not systematically attach to a power series a pointwise quantitative radius measuring the entry index into an ε -tube around the limiting value [4,5,8,9].
The purpose of this paper is to fill this gap by introducing sequential radii of convergence for power series. The indexed radius ρ f , a ( x ; ε ) measures the number of terms required to reach accuracy ε at a fixed point x, while the geometric radius ρ f , a * ( x ; ε ) gives a reciprocal efficiency scale. Together, these quantities refine the classical radius by separating where a power series converges from how efficiently its partial sums converge there.

1.3. Organization of the Paper

This paper is organized as follows. Section 2 reviews the mathematical background needed for the sequel, including classical convergence, tail remainders, indexed and geometric radii for sequences, and the associated tail-classification theorem. Section 3 develops the theory of sequential radii of convergence for power series in the space P a , viewed as a natural extension of the corresponding theory for numerical sequences. First, Section 3.1 introduces the sequential radii of convergence for power series, proves the main classification theorem, and defines the three associated tail-class subspaces of P a . Second, Section 3.2 describes the geometry of the sequential radii, emphasizing both the indexed and geometric viewpoints. Third, Section 3.3 studies the behavior of sequential radii under algebraic operations. Fourth, Section 3.4 presents a catalog of examples and explicit computations for two important tail-class subclasses of P a : the linear and superlinear regimes. Finally, Section 4 concludes the paper by summarizing its main contributions, discussing its limitations, and outlining possible directions for future work.

2. Preliminaries

The reader who has studied real analysis is well acquainted with the following background theory of sequences and series [4,5,6,10]:

2.1. Classical Convergence and Tail Remainders

Definition 1 
(Series)). Given a sequence { c n } n 1 of real numbers, its associated series is the symbolic sum n = 1 c n = c 1 + c 2 + c 3 + . Its Nth partial sum is s N = n = 1 N c n and the tail reminder is r N = n = N + 1 + c n . The series converges to s if lim N s N = s , and in this case we write n = 1 c n = s ; otherwise, the series diverges.
Definition 2 
(Power Series). Given a sequence of real coefficients { c n } n 0 and a center a R , a power series around a is a series of the form f ( x ) = n = 0 c n ( x a ) n , x R . Furthermore, its tail is defined by R N ( x ) = n = N + 1 c n ( x a ) n . Its convergence depends on x, and the set of convergence is described by a radius of convergence R [ 0 , ] around the center a.
Definition 3 
(Taylor Series). A Taylor series is a power series whose coefficients are determined by the derivatives of a infinitely differentiable function f at a center x = a given by T f , a ( x ) = n = 0 f ( n ) ( a ) n ! ( x a ) n .
Definition 4 
(Space of Real Power Series at a Center). For a fixed center a R , let P a denote the set of all real power series f ( x ) = n = 0 c n ( x a ) n with classical radius of convergence R f > 0 . For f , g P a with coefficient sequences { c n } n 0 and { d n } n 0 , respectively, and for λ R , the following operations are defined on the common open interval of convergence | x a | < min { R f , R g } :
1.
Scalar multiplication. ( λ f ) ( x ) : = n = 0 ( λ c n ) ( x a ) n .
2.
Sum and difference. ( f ± g ) ( x ) : = n = 0 ( c n ± d n ) ( x a ) n .
3.
Pointwise product. ( f g ) ( x ) : = f ( x ) g ( x ) , defined as the ordinary product of values on the common interval of convergence. As a function, f g coincides with the Cauchy product f * g (cf. item 5 below).
4.
Synchronized product approximants. For each n 0 , set p n ( x ) : = s n f ( x ) s n g ( x ) , and p ( x ) : = f ( x ) g ( x ) . The sequence { p n ( x ) } n 0 is a polynomial approximant of degree 2 n to f g and is, in general, distinct from the Cauchy partial sums of f * g , which have degree n.
5.
Cauchy product. ( f * g ) ( x ) : = n = 0 e n ( x a ) n , where e n : = k = 0 n c k d n k .
6.
Hadamard product. ( f g ) ( x ) : = n = 0 c n d n ( x a ) n .
7.
Quotient. When g ( a ) = d 0 0 , the quotient ( f / g ) ( x ) : = f ( x ) / g ( x ) is well-defined on a neighborhood of a.
8.
Termwise differentiation and integration. f ( x ) : = n = 1 n c n ( x a ) n 1 and ( I f ) ( x ) : = n = 0 c n n + 1 ( x a ) n + 1 .
9.
Argument transformations. For λ R { 0 } and k N , ( σ λ f ) ( x ) : = f a + λ ( x a ) and ( τ k f ) ( x ) : = f a + ( x a ) k .
Remark 1 
(Distinction among the three products). The pointwise product f g , the Cauchy product f * g , and the Hadamard product f g coincide as formal power series only in degenerate cases, and in general they produce distinct coefficient structures and distinct sequential radii of convergence [4,11]. This distinction is essential in Section 3.3 (cf. Remark 6).
Theorem 1 
(Classical Convergence Tests for Series). Let n = 1 c n be a series.
1.
Cauchy’s Root Test. Define
α = lim sup n | c n | n .
If α < 1 , then c n converges absolutely; if α > 1 , then c n diverges; and if α = 1 , the test is inconclusive.
2.
d’Alembert’s Ratio Test. Suppose c n 0 eventually and define
β = lim sup n c n + 1 c n .
If β < 1 , then c n converges absolutely; if β > 1 , then c n diverges; and the borderline case β = 1 is inconclusive.
Theorem 2 
(Classical Radius of Convergence for a Power Series). Let f ( x ) = n = 0 c n ( x a ) n , a R , be a real power series. Define α = lim sup n | c n | n . Then the classical radius of convergence is R = 1 α , with the conventions 1 / 0 = and 1 / = 0 . Hence f ( x ) converges absolutely for | x a | < R and diverges for | x a | > R . The endpoint cases x = a ± R when 0 < R < , must be checked separately. If, moreover, the ratio limit β = lim n c n + 1 c n , exists, then equivalently R = 1 β .
Remark 2 
(Classical radius dependency on spatial position ( x ) ). Let f ( x ) = n = 0 c n ( x a ) n be a real power series, with associated tail R n ( x ) and classical radius of convergence R . Then, inside the convergence interval, the asymptotic cost of making the tail small depends on how fast the remainder R n ( x ) decays. In particular, the position x plays an essential role only in the genuinely geometric case, where the effective ratio is comparable to | x a | / R (See Table 1).
The following asymptotic results are later frequently used in examples in this paper [12,13,14,15]:
Theorem 3 
(The Stirling’s Approximation). For any positive integer n, as n :
n ! 2 π n n e n ,
ln ( n ! ) = n ln n n + O ( ln n ) .
Theorem 4 
(Asymptotic Expansion of Harmonic Numbers). Let H n = k = 1 n 1 k denote the n-th harmonic number, and let γ 0.57721 be the Euler–Mascheroni constant. As n , the partial sums satisfy the logarithmic approximation:
H n = ln ( n ) + γ + O 1 n .
Theorem 5 
(Euler Number Asymptotic Formula). Let E 2 n denote the even-indexed Euler numbers. As n , the growth rate of the normalized coefficients satisfies the asymptotic behavior:
| E 2 n | ( 2 n ) ! 4 n + 1 π 2 n + 1
Theorem 6 
(Bernoulli Number Asymptotic Formula). Let B 2 n denote the even-indexed Bernoulli numbers, and let n . Then, the asymptotic behavior is given by:
B 2 n = ( 1 ) n 1 2 ( 2 n ) ! ( 2 π ) 2 n ζ ( 2 n ) : ζ ( 2 n ) 1 .
Theorem 7 
(Hardy–Ramanujan Asymptotic Formula). Let p ( n ) denote the number of partitions of a positive integer n. As n , the partition function satisfies the asymptotic behavior:
p ( n ) 1 4 3 n exp ( π 2 n 3 ) .

2.2. The Indexed Radius of Convergence and Tail Errors

The following notions are adapted from the indexed, geometric, and Cauchy radii introduced for numerical series in [10,12]. They are recalled here in a self-contained form, since the power-series construction in Section 3 is obtained by applying them pointwise to the associated tail sequence R N ( x ) .
Definition 5 
(Indexed Radius of convergence of a series). Let s = n = 1 c n be a real series, and let s N : = n = 1 N c n denote its N-th partial sum. Suppose that the partial sums have limit S ( s ) R { , + } . For ε > 0 , define the radius of convergence of s at level ε by
ρ s ( ε ) : = inf R ( s , ε ) ,
where
R ( s , ε ) : = N N : n N , s n < ε , if S ( s ) = , N N : n N , | s n S ( s ) | < ε , if S ( s ) R , N N : n N , s n > ε , if S ( s ) = + .
Definition 6 
(Geometric radius of convergence of a series). Let s = n = 1 c n be a real series, and let ρ s ( ε ) denote its index-radius of convergence at level ε > 0 . To emphasize the analogy with geometric notions of radius, one may equivalently work with the rescaled radius
ρ s * ( ε ) : = 1 ρ s ( ε ) 1 , ε > 0 ,
with the convention 1 / 0 : = + .
Definition 7 
(Cauchy radius of a series). Let s = n = 1 c n be a real series, and let s N : = k = 1 N c k denote its N-th partial sum. For ε > 0 , define the Cauchy radius of s at level ε by
ρ s C ( ε ) : = inf C ( s , ε ) ,
where
C ( s , ε ) : = N N : m , n N , | s n s m | < ε .
Theorem 8 
(Tail-Inversion Principle for Series). Let s = n = 1 c n be a convergent real series with finite sum S R . Let s n = k = 1 n c k , R n = S s n = k = n + 1 c k , n 1 , denote the partial sums and the corresponding tails. Suppose that there exists a strictly decreasing, continuous, and invertible function f : ( 0 , ) ( 0 , ) such that | R n | f ( n ) ( n ) . Then,
(i)
ρ s ( ε ) f 1 ( ε ) , ε 0 ,
(ii)
ρ s * ( ε ) 1 f 1 ( ε ) , ε 0 .
Corollary 1 
(Tail-decay classes and induced convergence radii). Let s = n = 1 c n be a convergent real series with finite sum S R . For n 1 , write s n = k = 1 n c k , R n = S s n = k = n + 1 c k . Let ρ s ( ε ) and ρ s * ( ε ) denote, respectively, the index-radius and the geometric radius of convergence of the partial-sum sequence s = ( s n ) . As ε 0 , the asymptotic decay of the tail R n determines the corresponding radius class as summarized in Table 2.

3. Theory of Sequential Radii of Convergence for Power Series

The preceding section recalled the indexed and geometric radii for numerical series. We now transfer this point of view to power series by fixing an evaluation point ( x ) inside the classical interval of convergence and applying the numerical theory to the corresponding tail sequence ( R N ( x ) ) . This pointwise reduction is the basic mechanism behind the sequential radii of convergence for power series.

3.1. The General Classification Theorem

Definition 8 
(Indexed radius of convergence of a power series). Let f ( x ) = n = 0 c n ( x a ) n , a R , be a power series, and for each N 0 define the N-th partial sum and tail by s N ( x ) : = n = 0 N c n ( x a ) n , R N ( x ) : = f ( x ) s N ( x ) = n = N + 1 c n ( x a ) n . For a fixed x at which the power series converges to f ( x ) , and for ε > 0 , define the indexed radius of convergence at x by
ρ f , a ( x ; ε ) : = inf N N : | R n ( x ) | < ε for all n N .
Remark 3 
(Geometric and Cauchy radii for a power series). Let f ( x ) = n = 0 c n ( x a ) n , and s N ( x ) : = n = 0 N c n ( x a ) n . For a fixed x R and ε > 0 , the following two associated radii are useful.
1.
Geometric radius. If ρ f , a ( x ; ε ) denotes the indexed radius of convergence of the power series at x, define
ρ f , a * ( x ; ε ) : = 1 ρ f , a ( x ; ε ) 1 ,
with the convention 1 / 0 : = + .
2.
Cauchy radius. Define the Cauchy radius of the power series at x by
ρ f , a C ( x ; ε ) : = inf C f , a ( x ; ε ) ,
where
C f , a ( x ; ε ) : = N N : m , n N , | s n ( x ) s m ( x ) | < ε .
Definition 9 
(Sequential Radii of Convergence for power series). Let f ( x ) = n = 0 c n ( x a ) n , a R , be a power series, and ε > 0 . We refer to the pair ρ f , a ( x ; ε ) , ρ f , a * ( x ; ε ) , as its sequential radii at given fixed point x a . .
Theorem 9 
(Sequential Radii of Convergence for Power Series). Let f ( x ) = n = 0 c n ( x a ) n , a R , be a power series with classical radius of convergence R ( 0 , ] . Fix x with | x a | < R , and define its associated nth partial sum and tail s n ( x ) : = k = 0 n c k ( x a ) k , R n ( x ) : = f ( x ) s n ( x ) = k = n + 1 c k ( x a ) k , respectively. Assume that, as n , the tail R n ( x ) belongs to one of the following decay regimes. Then, as ε 0 , the corresponding indexed and geometric radii are given by
Tail regime Tail asymptotic Indexed radius ρ f , a ( x ; ε ) Geometric radius ρ f , a * ( x ; ε ) Sublinear | R n ( x ) | C 1 n ( p 1 ) , C 1 > 0 , p > 1 C 1 ε 1 / ( p 1 ) ε C 1 1 / ( p 1 ) Linear | R n ( x ) | C 2 r x n , r x : = | x a | R ( 0 , 1 ) log ( C 2 / ε ) log ( R / | x a | ) log ( R / | x a | ) log ( C 2 / ε ) Superlinear | R n ( x ) | C 3 ( n + 1 ) ! , C 3 > 0 log ( C 3 / ε ) log log ( C 3 / ε ) log log ( C 3 / ε ) log ( C 3 / ε )
In particular, in the linear regime, the spatial dependence of the indexed radius enters through the classical relative position r x = | x a | R . If x = a , then R n ( a ) = 0 for all n 0 , and therefore ρ f , a ( a ; ε ) = 0 , ρ f , a * ( a ; ε ) = + , respectively.
Proof. 
Fix x with | x a | < R , and set b n ( x ) : = c n ( x a ) n . Then n = 0 b n ( x ) = n = 0 c n ( x a ) n = f ( x ) , and its n-th tail is exactly k = n + 1 b k ( x ) = k = n + 1 c k ( x a ) k = R n ( x ) . Hence the indexed radius ρ f , a ( x ; ε ) is precisely the indexed radius of the fixed numerical series n = 0 b n ( x ) . The three asymptotic formulas for ρ f , a ( x ; ε ) therefore follow directly from Corollary 1 applied to the numerical tail R n ( x ) . The formulas for the geometric radius follow from the definition in Remark 3. Finally, in the linear power-series regime, the decay ratio is the relative distance of x from the center of expansion, namely r x = | x a | R , which gives the stated spatial dependence. The case x = a is immediate, since all powers ( x a ) n vanish for n 1 , so the tail is identically zero. □
Theorem 9 identifies three model regimes for the decay of the tail ( R N ( x ) ) . It is therefore natural to separate those power series whose tails fall into one of these regimes at a fixed point ( x ) . The following definition records this separation as a tail-class decomposition inside the space ( P a ) .
Definition 10 
(Tail-class subspaces of P a ). Fix a center a R , and let P a be the space of real power series centered at a in the sense of Definition 4. Thus P a : = f ( x ) = n = 0 c n ( x a ) n : R f > 0 , where R f denotes the classical radius of convergence of f. For f P a , define the N-th partial sum and tail by s N f ( x ) : = n = 0 N c n ( x a ) n , and R N f ( x ) : = f ( x ) s N f ( x ) = n = N + 1 c n ( x a ) n , respectively. For a fixed point x at which the power series converges, define the following tail-class subspaces of P a :
1 . Sublinear class . P a sub ( x ) : = C 1 > 0 p > 1 f P a : | R N f ( x ) | C 1 N ( p 1 ) ( N ) .
2 . Linear class . P a lin ( x ) : = C 2 > 0 0 < r x < 1 f P a : | R N f ( x ) | C 2 r x N ( N ) .
3 . Superlinear class . P a sup ( x ) : = C 3 > 0 f P a : | R N f ( x ) | C 3 ( N + 1 ) ! ( N ) .
The corresponding tail-class decomposition is therefore:
P a seq ( x ) : = P a sub ( x ) P a lin ( x ) P a sup ( x ) .
Here, P a seq ( x ) denotes the subclass of power series in P a whose tail at the fixed evaluation point x belongs to one of the three model regimes appearing in Theorem 9.
Remark 4 
(Structural Consequences). Theorem 9 has several immediate consequences as follows:
(a)
Interpretation of the tail constant C: The constant C is the leading amplitude of the tail asymptotic | R n ( x ) | C ϕ n , where ϕ n = n ( p 1 ) , r x n , or ( n + 1 ) ! 1 in the sublinear, linear, and superlinear regimes, respectively. While C is crucial for calculating the exact error at a specific index n , the Theorem shows that C is irrelevant for the asymptotic growth rate of the sequential radii of convergence. Figure 1(a) visualizes that changing C = C 2 shifts the indexed radius curve but does not change the asymptotic growth rate order.
(b)
Convergence-speed Classification: The Theorem classifies sequential convergence by the decay of the geometric radius ρ f , a * ( x ; ε ) : sublinear tails give a polynomial scale ε 1 / ( p 1 ) , linear tails give the logarithmic scale log ( R / | x a | ) / log ( 1 / ε ) , and superlinear tails give the faster efficiency scale log log ( 1 / ε ) / log ( 1 / ε ) . Thus the table separates polynomial, geometric, and factorial tail behavior into three distinct radius classes, while in all nondegenerate cases ρ f , a * ( x ; ε ) 0 as ε 0 , reflecting the continuous narrowing of the admissible geometric window. Figure 1(b) compares the three radius classes: polynomial/sublinear, logarithmic/linear, and factorial/superlinear.
(c)
Geometry of Sequential Convergence: The spatial geometry of the Theorem appears explicitly in the linear regime through r x = | x a | / R , where ρ f , a ( x ; ε ) log ( 1 / ε ) / log ( R / | x a | ) . Thus convergence is fastest near the center x = a , where the tail is degenerate, and becomes arbitrarily slow as | x a | R , reflecting the increasing cost of approximation near the boundary of the classical disk of convergence. Figure 1(c) shows the explosion of the indexed radius as r = | x a | approaches the boundary R .

3.2. The Geometry of the Sequential Radii

The preceding section summarized the structural consequences of Theorem 9, including the role of the tail constant, the separation of convergence-speed classes, and the spatial dependence of the linear regime. The last of these is intrinsically geometric, since the sequential radii depend on the relative position r x = | x a | / R inside the classical interval of convergence. While Figure 1(c) showed this dependence along a one-dimensional radial slice, the present section lifts the same phenomenon to a two-dimensional surface over x and log ( 1 / ε ) .
Figure 2 further illustrates the two-dimensional geometry induced by the linear regime of Theorem 9. In this regime, with C 2 = 1 , R = 1 , and a = 0 , the indexed radius has the asymptotic form
ρ f , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / | x | ) , 0 < | x | < 1 .
Panel (a) shows the corresponding indexed-radius surface over the spatial variable x and the precision scale log ( 1 / ε ) . The surface has a funnel-like geometry: near the center x = 0 , the approximation cost is small because the tail is degenerate at the center of expansion, while the indexed radius grows increasingly fast as | x | approaches the boundary R = 1 . Thus the same tolerance ε requires substantially more terms near the edge of the classical interval of convergence than near the center. This visualizes the spatial dependence encoded in the denominator log ( 1 / | x | ) , which becomes small as | x | 1 .
Panel (b) displays the associated geometric radius
ρ f , 0 * ( x ; ε ) log ( 1 / | x | ) log ( 1 / ε ) , 0 < | x | < 1 .
This surface gives the reciprocal geometric interpretation of the same phenomenon. While the indexed radius measures the number of terms required for a given accuracy, the geometric radius measures the corresponding normalized efficiency scale. Hence the surface is largest near the center of expansion, where convergence is most efficient, and it collapses toward zero both as the precision demand increases and as x moves toward the boundary of convergence. Together, the two panels show that the linear regime is not only a tail-decay class, but also a spatial geometry: convergence becomes progressively more expensive as the evaluation point approaches the classical boundary | x a | = R .

3.3. Algebraic Operations

The preceding sections studied the sequential radii of a single power series through the decay of its tail sequence. We now record how these radii behave under elementary algebraic operations. The point of view is local in the evaluation point x: after fixing x inside the common interval of convergence, each power series gives a convergent sequence of partial sums, and the algebraic operations reduce to corresponding operations on these partial-sum sequences. Thus the estimates below should be read as stability bounds for the computational cost of approximation under scalar multiplication, addition, multiplication, and division.
Throughout this subsection, fix a center a R and consider two power series f , g P a , with f ( x ) = n = 0 c n ( x a ) n , and g ( x ) = n = 0 d n ( x a ) n , and classical radii of convergence R f and R g , respectively. For each n 0 , denote their n-th partial sums and tails by s n f ( x ) : = k = 0 n c k ( x a ) k , and R n f ( x ) : = f ( x ) s n f ( x ) , and analogously s n g ( x ) , R n g ( x ) , and let ρ f , a ( x ; ε ) , ρ g , a ( x ; ε ) denote the associated indexed radii in the sense of Definition 5. Throughout, x is fixed in the common open interval of convergence | x a | < min { R f , R g } . Under this notation, the following theorem records the stability of the sequential radii under the algebraic operations of Definition 4. We use this notation in the following Theorem:
Theorem 10 
(Algebraic stability of sequential radii). Let f and g be two real power series with common center a, and fix x in the common open interval of convergence. Then the following estimates hold for every ε > 0 .
(i) 
Scalar multiplication. Let λ R and h = λ f . If λ 0 , then ρ h , a ( x ; ε ) = ρ f , a x ; ε | λ | . If λ = 0 , then h 0 and ρ h , a ( x ; ε ) = 0 .
(ii) 
Sum and difference. Let h = f ± g . Then ρ h , a ( x ; ε ) max ρ f , a x ; ε 2 , ρ g , a x ; ε 2 . More generally, for every t ( 0 , 1 ) , we have ρ f ± g , a ( x ; ε ) max ρ f , a ( x ; t ε ) , ρ g , a ( x ; ( 1 t ) ε ) .
(iii) 
Pointwise product of synchronized approximants. Define the synchronized product approximants by p n ( x ) : = s n f ( x ) s n g ( x ) , and p ( x ) : = f ( x ) g ( x ) . Let r × ( ε ; x ) : = min 1 , ε 1 + | f ( x ) | + | g ( x ) | . Then ρ f g , a pt ( x ; ε ) max ρ f , a ( x ; r × ( ε ; x ) ) , ρ g , a ( x ; r × ( ε ; x ) ) , where ρ f g , a pt denotes the indexed radius associated with the approximating sequence p n ( x ) = s n f ( x ) s n g ( x ) .
(iv) 
Quotient. Assume g ( x ) 0 . Define the synchronized quotient approximants by q n ( x ) : = s n f ( x ) s n g ( x ) , and q ( x ) : = f ( x ) g ( x ) . Then s n g ( x ) 0 for all sufficiently large n, and q n ( x ) q ( x ) . For every ε > 0 , set r ÷ ( ε ; x ) : = min { | g ( x ) | 2 , ε | g ( x ) | 2 2 | f ( x ) | + | g ( x ) | } . Then ρ f / g , a pt ( x ; ε ) max ρ f , a ( x ; r ÷ ( ε ; x ) ) , ρ g , a ( x ; r ÷ ( ε ; x ) ) .
Proof. 
The proof is divided according to the listed assertions. Fix x in the common open interval of convergence. Then s n f ( x ) f ( x ) , and s n g ( x ) g ( x ) , and the indexed radii ρ f , a ( x ; ε ) and ρ g , a ( x ; ε ) are precisely the entry indices of these two partial-sum sequences into their corresponding ε -tubes.
(i). For scalar multiplication, if h = λ f , then R n h ( x ) = h ( x ) s n h ( x ) = λ f ( x ) s n f ( x ) = λ R n f ( x ) . Hence | R n h ( x ) | < ε if and only if | R n f ( x ) | < ε | λ | when λ 0 , which gives the stated identity. If λ = 0 , the tail is identically zero, so the radius is 0.
(ii). For sums and differences, write h = f ± g . Then, R n h ( x ) = R n f ( x ) ± R n g ( x ) . Therefore, if n max ρ f , a ( x ; t ε ) , ρ g , a ( x ; ( 1 t ) ε ) , then | R n h ( x ) | | R n f ( x ) | + | R n g ( x ) | < t ε + ( 1 t ) ε = ε . Taking t = 1 / 2 gives the symmetric estimate.
(iii). For the synchronized product, write
s n f ( x ) s n g ( x ) f ( x ) g ( x ) = s n f ( x ) f ( x ) s n g ( x ) g ( x ) + f ( x ) s n g ( x ) g ( x ) + g ( x ) s n f ( x ) f ( x ) .
If n max ρ f , a ( x ; r × ) , ρ g , a ( x ; r × ) , then | s n f ( x ) f ( x ) | < r × , and | s n g ( x ) g ( x ) | < r × . Since r × 1 , we obtain
| s n f ( x ) s n g ( x ) f ( x ) g ( x ) | r × 2 + | f ( x ) | + | g ( x ) | r × 1 + | f ( x ) | + | g ( x ) | r × ε .
This proves the product estimate.
(iv). Finally, assume g ( x ) 0 . We accomplish the proof in few steps as in follows.
Step 1. We consider the synchronized quotient approximants q n ( x ) : = s n f ( x ) s n g ( x ) , and q ( x ) : = f ( x ) g ( x ) . Since s n g ( x ) g ( x ) 0 , the denominators s n g ( x ) are eventually nonzero. More precisely, if | s n g ( x ) g ( x ) | < | g ( x ) | 2 , then, by the reverse triangle inequality, | s n g ( x ) | | g ( x ) | | s n g ( x ) g ( x ) | > | g ( x ) | | g ( x ) | 2 = | g ( x ) | 2 . Thus, whenever the tail error of g is smaller than | g ( x ) | / 2 , the denominator s n g ( x ) is bounded away from zero. Moreover, 1 s n g ( x ) 1 g ( x ) = | s n g ( x ) g ( x ) | | s n g ( x ) | | g ( x ) | . Using | s n g ( x ) | | g ( x ) | / 2 , we obtain 1 | s n g ( x ) | 2 | g ( x ) | , and therefore
1 s n g ( x ) 1 g ( x ) 2 | g ( x ) | 2 | s n g ( x ) g ( x ) | .
Step 2. Now write
s n f ( x ) s n g ( x ) f ( x ) g ( x ) = s n f ( x ) f ( x ) s n g ( x ) + f ( x ) 1 s n g ( x ) 1 g ( x ) .
Taking absolute values gives
s n f ( x ) s n g ( x ) f ( x ) g ( x ) | s n f ( x ) f ( x ) | | s n g ( x ) | + | f ( x ) | 1 s n g ( x ) 1 g ( x ) .
Hence, by inequality (22)
s n f ( x ) s n g ( x ) f ( x ) g ( x ) 2 | g ( x ) | | s n f ( x ) f ( x ) | + 2 | f ( x ) | | g ( x ) | 2 | s n g ( x ) g ( x ) | .
Step 3. We now choose a common error tolerance that guarantees both the denominator condition and the final ε -bound. Define r ÷ ( ε ; x ) : = min | g ( x ) | 2 , ε | g ( x ) | 2 2 ( | f ( x ) | + | g ( x ) | ) . Let N : = max ρ f , a ( x ; r ÷ ( ε ; x ) ) , ρ g , a ( x ; r ÷ ( ε ; x ) ) . Then for every n N , we have | s n f ( x ) f ( x ) | < r ÷ ( ε ; x ) , and | s n g ( x ) g ( x ) | < r ÷ ( ε ; x ) . Accordingly, by inequality (23):
s n f ( x ) s n g ( x ) f ( x ) g ( x ) < 2 | g ( x ) | + 2 | f ( x ) | | g ( x ) | 2 r ÷ ( ε ; x ) = 2 ( | f ( x ) | + | g ( x ) | ) | g ( x ) | 2 r ÷ ( ε ; x ) < ε .
Therefore, N R f / g , a pt ( x ; ε ) , and hence ρ f / g , a pt ( x ; ε ) max ρ f , a ( x ; r ÷ ( ε ; x ) ) , ρ g , a ( x ; r ÷ ( ε ; x ) ) . This proves the quotient estimate. □
Remark 5 
(Interpretation of the algebraic bounds). The estimates in Theorem 10 show that algebraic operations do not destroy sequential convergence inside the common interval of convergence. Scalar multiplication simply rescales the tolerance. Sums and differences require both input tails to be small, so the resulting radius is controlled by the worse of the two entry indices. Products and quotients introduce dependence on the limiting values f ( x ) and g ( x ) , because the error propagation contains mixed terms. In this sense, the sequential radii behave like computational stability profiles for algebraic operations on power-series approximations.
Remark 6 
(Pointwise vs. Cauchy vs. Hadamard multiplication). The following observations distinguish three notions of multiplication on power series.
(a) 
The pointwise product f g and the Cauchy product f * g coincide as functions on the common open interval of convergence, with shared coefficient structure e n = k = 0 n c k d n k . The Hadamard product f g is a genuinely distinct power series, with coefficients c n d n and limit generally different from f ( x ) g ( x ) .
(b) 
Although f g and f * g agree as functions, their natural approximating sequences differ. The synchronized pointwise approximants are p n ( x ) : = s n f ( x ) s n g ( x ) , of degree 2 n , whereas the Cauchy-product partial sums are S n f * g ( x ) : = j = 0 n e j ( x a ) j , of degree n. Both sequences converge to f ( x ) g ( x ) on the common interval of convergence, but enter the ε-tube at generally distinct indices. Hence their indexed and geometric radii are generally distinct.
(c) 
For the Cauchy and Hadamard cases, the tails R n f * g ( x ) and R n f g ( x ) both vanish on the respective intervals of convergence, so ρ f * g , a ( x ; ε ) < and ρ f g , a ( x ; ε ) < for every ε > 0 . In general, however, no universal formula expresses ρ f * g , a ( x ; ε ) or ρ f g , a ( x ; ε ) in terms of ρ f , a ( x ; ε ) and ρ g , a ( x ; ε ) alone; once an asymptotic estimate for the relevant tail is available, the corresponding radius follows from Theorem 9. Table 3 presents the list of counterexamples.
Table 3. Counterexamples supporting Remark 3.8(c): twin pairs ( f 1 , g 1 ) and ( f 2 , g 2 ) exhibiting identical individual indexed radii but distinct product radii. Parameters: a = 0 , x 0 = 1 / 2 ; asymptotics are stated as ε 0 .
Table 3. Counterexamples supporting Remark 3.8(c): twin pairs ( f 1 , g 1 ) and ( f 2 , g 2 ) exhibiting identical individual indexed radii but distinct product radii. Parameters: a = 0 , x 0 = 1 / 2 ; asymptotics are stated as ε 0 .
(a) Cauchy product fg.
# f ( x ) g ( x ) ( f * g ) ( x ) ρ f , 0 ( x 0 ; ε ) ρ g , 0 ( x 0 ; ε ) ρ f * g , 0 ( x 0 ; ε )
1 n = 0 x n n ! n = 0 ( x ) n n ! 1 log ( C f / ε ) log log ( C f / ε ) log ( C g / ε ) log log ( C g / ε ) 0
2 n = 0 x n n ! n = 0 x n n ! n = 0 ( 2 x ) n n ! log ( C f / ε ) log log ( C f / ε ) log ( C g / ε ) log log ( C g / ε ) log ( C f * g / ε ) log log ( C f * g / ε )
(b) Hadamard product fg.
# f ( x ) g ( x ) ( f g ) ( x ) ρ f , 0 ( x 0 ; ε ) ρ g , 0 ( x 0 ; ε ) ρ f g , 0 ( x 0 ; ε )
1 n = 0 x 2 n n = 0 x 2 n + 1 0 log ( C f / ε ) log 2 log ( C g / ε ) log 2 0
2 n = 0 x 2 n n = 0 x 2 n n = 0 x 2 n log ( C f / ε ) log 2 log ( C g / ε ) log 2 log ( C f g / ε ) log 2
Notes. Asymptotic equivalences (~) are taken in the sense of Theorem 9 applied to the appropriate tail; scaled exponential series fw(x) = ewx =   n = 0 ( w x ) n n ! , (w ≠ 0), has superlinear tail regime; lacunary-type geometric series fa,b(x) =  x a 1 - x b  =   n = 0 x a + b n , (a, b > 0), has linear tail regime.

3.4. A Catalog of Classical Power Series: Examples and Explicit Sequential-Radius Computations

This section records a selected catalog of classical power series and computes their sequential radii in the sense of Definition 8 and Remark 3. This list is not exhaustive rather it is chosen to represent often historically important series expansions, different classical radii of convergence, and different coefficient structures. Throughout this subsection, the center is taken a = 0 unless otherwise stated, and x is fixed in the interior of the classical interval of convergence. In the next two subsections, we use the preceding lemmas to identify general patterns in the tail behavior of our examples [6,13,14].

3.4.1. Linear/Geometric Series: Class 2

The series discussed in this section have finite classical radius 0 < R < . we write r x = | x | R ( 0 , 1 ) . The phrase "linear class" means that the tail is geometric up to subexponential or algebraic factors in the index. This assumption slightly broadens the strict hypothesis R n ( x ) C r x n in Theorem 9, but only at lower order: factors such as n α , ( log n ) β , exp ( C n ) , or the sparsity of nonzero coefficients may contribute secondary corrections, while the leading indexed-radius remains as in Theorem 9. We start with the first key pattern lemma [12,16]:
Lemma 1 
(Series Tail Pattern(1)). Let f ( x ) = n = 0 + c n x n be a real power series with finite positive classical convergence radius 0 < R < + . Suppose, its tail satisfy:
| R N ( x ) | C ( x ) N α ( log ( N ) ) β exp ( A N ) r N ( N + )
where C ( x ) > 0 is a function, 0 < r < 1 , and α , β , N 0 . Then:
ρ f , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) : r = | x | R , ( ϵ 0 ) .
Proof. 
Set L : = log ( 1 / ε ) , and γ : = log ( 1 / r ) > 0 . Since r N = e γ N , taking logarithms gives:
log | R N ( x ) | = log C ( x ) + α log N + β log log N + A N γ N + o ( 1 ) .
Moreover, log N = o ( N ) , log log N = o ( N ) , and N = o ( N ) , and therefore
log | R N ( x ) | = γ N + o ( N ) .
We now invert this estimate in both directions. Fix δ ( 0 , 1 ) and define N + ( ε ) : = ( 1 + δ ) L γ . At this scale N + = O ( L ) , hence o ( N + ) = o ( L ) . Consequently, log | R N + ( x ) | = γ N + + o ( N + ) = ( 1 + δ ) L + o ( L ) < L for sufficiently small ε . Thus | R N + ( x ) | < e L = ε . Since the geometric factor dominates the lower-order factors, the tail is eventually decreasing in magnitude, and therefore N + is an admissible entry index. Hence
ρ f , 0 ( x ; ε ) ( 1 + δ ) L γ + O ( 1 ) .
For the reverse estimate, define N ( ε ) : = ( 1 δ ) L γ . Again N = O ( L ) , and hence log | R N ( x ) | = γ N + o ( N ) = ( 1 δ ) L + o ( L ) > L for sufficiently small ε . Thus | R N ( x ) | > ε , so no index N N can be admissible. Therefore
ρ f , 0 ( x ; ε ) ( 1 δ ) L γ + O ( 1 ) .
Combining the upper and lower estimates in (27) and (28) gives
1 δ lim inf ε 0 ρ f , 0 ( x ; ε ) L / γ lim sup ε 0 ρ f , 0 ( x ; ε ) L / γ 1 + δ .
Letting δ 0 in (29) yields
ρ f , 0 ( x ; ε ) L γ = log ( 1 / ε ) log ( 1 / r ) .
Lemma 2 
(Tail pattern for geometric series). Let f ( x ) = n = 0 c n x n be a real power series with finite positive classical radius of convergence 0 < R < . Assume that the coefficients satisfy c n K n α ( log n ) β exp { A n } R n , K 0 , for some parameters α , β , A R . Fix an interior point 0 < | x | < R . Then, as N ,
| R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N : r = | x | R , q = x R , C ( x ) = | K | r | 1 q | .
Proof. 
The proof is divided into several steps.
Step 1. For the fixed point x, write the n-th term of the power series as
u n ( x ) : = c n x n ( K n α ( log n ) β exp { A n } R n ) x n ( K n α ( log n ) β exp { A n } q n = K L ( n ) q n : L ( n ) : = n α ( log n ) β exp { A n } .
Step 2. Since | q | = r < 1 , the geometric factor q n dominates the algebraic, logarithmic, and subexponential factors in L ( n ) . In particular, for every fixed k 0 ,
L ( N + 1 + k ) L ( N + 1 ) = N + 1 + k N + 1 α log ( N + 1 + k ) log ( N + 1 ) β exp A N + 1 + k N + 1 1 × 1 × 1 = 1 , ( N ) .
Step 3. We now use (31) and (32) to express the tail in terms of the first omitted term:
R N ( x ) = n = N + 1 u n ( x ) = u N + 1 ( x ) k = 0 u N + 1 + k ( x ) u N + 1 ( x ) u N + 1 ( x ) k = 0 ( K L ( N + 1 + k ) q N + 1 + k K L ( N + 1 ) q N + 1 ) u N + 1 ( x ) k = 0 ( L ( N + 1 + k ) L ( N + 1 ) ) ( q k ) u N + 1 ( x ) k = 0 ( q k ) u N + 1 ( x ) 1 q : | q | < 1 .
Step 4. It remains to rewrite the first omitted term in the desired N-scale. Since by (31) we have
u N + 1 ( x ) K ( N + 1 ) α ( log ( N + 1 ) ) β exp { A N + 1 } q N + 1 , K ( ( N ) α ) ( ( log ( N ) ) β ) ( exp { A N } ) q N + 1 , = K q N α ( log N ) β exp { A N } q N .
Therefore, by (33) and (34) we have:
| R N ( x ) | K q 1 q N α ( log N ) β exp { A N } | q | N .
This proves the claim. □
Example 1 
(Lacunary geometric series). The lacunary geometric function has the power-series expansion x a 1 x b = n = 0 x b n + a , a 0 , b N . Fix 0 < | x | < 1 , and write y = x b . Then x a 1 x b = x a n = 0 y n . This is the geometric pattern of Lemma 1 in the variable y, with ( K = 1 , α = 0 , β = 0 , A = 0 , R y = 1 ) . Therefore ( q = y = x b , r = | y | = | x | b ) . By Lemma 2, the lacunary-index tail satisfies | R N ( x ) | | x | a | x | b | 1 x b | | x | b N . Equivalently, | R N ( x ) | | x | a + b | 1 x b | | x | b N . Thus, by Lemma 1, ρ f , 0 lac ( x ; ε ) log ( 1 / ε ) log ( 1 / | x | b ) = 1 b log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 2 
(Reciprocal square series). The reciprocal square function has the power-series expansion 1 ( 1 x ) 2 = n = 1 n x n 1 = m = 0 ( m + 1 ) x m , | x | < 1 . Thus the coefficient sequence is c m = m + 1 , and therefore c m m . This is of the form in Lemma 2 with ( K = 1 , α = 1 , β = 0 , A = 0 , R = 1 ) . For fixed 0 < | x | < 1 , we have ( q = x R = x , r = | x | R = | x | ) , and C ( x ) = | K | r | 1 q | = | x | | 1 x | . Hence, by Lemma 2, | R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N . Therefore | R N ( x ) | | x | | 1 x | N | x | N . This is a linear/geometric tail with one algebraic factor. Applying Lemma 1, we obtain ρ f , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 3 
(Logarithm series). The logarithm function has the power-series expansion log ( 1 + x ) = n = 1 ( 1 ) n 1 n x n , | x | < 1 . Equivalently, setting y = x , we may write log ( 1 + x ) = n = 1 y n n . Thus, in the variable y, the coefficient sequence satisfies c n = 1 n n 1 . This is of the form in Lemma 2 with ( K = 1 , α = 1 , β = 0 , A = 0 , R = 1 ) . For fixed 0 < | x | < 1 , we have ( q = y R = x , r = | y | R = | x | ) Therefore C ( x ) = | K | r | 1 q | = | x | | 1 + x | . Hence, by Lemma 2, | R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N . Therefore | R N ( x ) | | x | | 1 + x | N 1 | x | N . This is a linear/geometric tail with one algebraic factor N 1 . Applying Lemma 1, we obtain ρ log ( 1 + x ) , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 4 
(Inverse hyperbolic tangent series). The inverse hyperbolic tangent function has the power-series expansion arctanh ( x ) = n = 0 x 2 n + 1 2 n + 1 , | x | < 1 . Fix 0 < | x | < 1 , and write y = x 2 . Then arctanh ( x ) = x n = 0 y n 2 n + 1 . The coefficient sequence in the variable y satisfies c n = 1 2 n + 1 1 2 n 1 . Therefore this is of the form in Lemma 2 with ( K = 1 2 , α = 1 , β = 0 , A = 0 , R y = 1 ) . For fixed 0 < | x | < 1 , we have ( q = y R y = x 2 , r = | y | R y = | x | 2 ) . Thus C y ( x ) = | K | r | 1 q | = | x | 2 2 | 1 x 2 | . After multiplying by the outer factor x, Lemma 2 gives | R N ( x ) | | x | C y ( x ) N 1 r N . Hence | R N ( x ) | | x | 3 2 | 1 x 2 | N 1 | x | 2 N . This is a lacunary linear/geometric tail with one algebraic factor N 1 . Applying Lemma 1, we obtain ρ arctanh , 0 lac ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | 2 ) = 1 2 log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 5 
(Harmonic generating function series). The harmonic generating function has the power-series expansion log ( 1 x ) 1 x = n = 1 H n x n , | x | < 1 , where H n = k = 1 n 1 k . By Theorem 4, H n = log n + γ + O 1 n , and therefore H n log n . Thus the coefficient sequence satisfies c n = H n log n . This is of the form in Lemma 2 with ( K = 1 , α = 0 , β = 1 , A = 0 , R = 1 ) . For fixed 0 < | x | < 1 , we have ( q = x R = x , r = | x | R = | x | ) . Therefore C ( x ) = | K | r | 1 q | = | x | | 1 x | . Hence, by Lemma 2, | R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N . Consequently, | R N ( x ) | | x | | 1 x | ( log N ) | x | N . This is a linear/geometric tail with one logarithmic correction. Applying Lemma 1, we obtain ρ f , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 6 
(Legendre chi series). The Legendre chi function has the power-series expansion χ s ( x ) = n = 0 x 2 n + 1 ( 2 n + 1 ) s , s = 2 , 3 , | x | < 1 . Fix 0 < | x | < 1 , and write y = x 2 . Then χ s ( x ) = x n = 0 y n ( 2 n + 1 ) s . The coefficient sequence in the variable y satisfies c n = 1 ( 2 n + 1 ) s 2 s n s . Therefore this is of the form in Lemma 2 with ( K = 2 s , α = s , β = 0 , A = 0 , R y = 1 ) . For fixed 0 < | x | < 1 , we have ( q = y R y = x 2 , r = | y | R y = | x | 2 ) . Thus, C y ( x ) = | K | r | 1 q | = | x | 2 2 s | 1 x 2 | . After multiplying by the outer factor x, Lemma 2 gives | R N ( x ) | | x | C y ( x ) N s r N . Hence | R N ( x ) | | x | 3 2 s | 1 x 2 | N s | x | 2 N . This is a lacunary linear/geometric tail with algebraic factor N s . Applying Lemma 1, we obtain ρ χ s , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | 2 ) = 1 2 log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 7 
(Arcsine series). The arcsine function has the power-series expansion arcsin ( x ) = n = 0   2 n n 4 n ( 2 n + 1 ) x 2 n + 1 , | x | < 1 . Fix 0 < | x | < 1 , and write y = x 2 . Then arcsin ( x ) = x n = 0 2 n n 4 n ( 2 n + 1 ) y n . Using the central-binomial asymptotic 2 n n 4 n π n , we obtain 2 n n 4 n ( 2 n + 1 ) 1 π n ( 2 n + 1 ) 1 2 π n 3 / 2 . Thus the coefficient sequence in the variable y satisfies c n = 2 n n 4 n ( 2 n + 1 ) 1 2 π n 3 / 2 . Therefore this is of the form in Lemma 2 with ( K = 1 2 π , α = 3 2 , β = 0 , A = 0 , R y = 1 ) . For fixed 0 < | x | < 1 , we have ( q = y R y = x 2 , r = | y | R y = | x | 2 ) . Thus C y ( x ) = | K | r | 1 q | = | x | 2 2 π | 1 x 2 | . After multiplying by the outer factor x, Lemma 2 gives | R N ( x ) | | x | C y ( x ) N 3 / 2 r N . Hence | R N ( x ) | | x | 3 2 π | 1 x 2 | N 3 / 2 | x | 2 N . This is a lacunary linear/geometric tail with algebraic factor N 3 / 2 . Applying Lemma 1, we obtain ρ arcsin , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | 2 ) = 1 2 log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 8 
(Binomial series). The generalized binomial function has the power-series expansion f λ ( x ) = ( 1 + x ) λ = n = 0 λ n x n , | x | < 1 , where λ Z . Using the standard coefficient asymptotic λ n ( 1 ) n Γ ( λ ) n λ 1 , n , we write y = x . Then the coefficient sequence in the variable y satisfies ( 1 ) n λ n 1 Γ ( λ ) n λ 1 . Therefore this is of the form in Lemma 2 with ( K = 1 Γ ( λ ) , α = λ 1 , β = 0 , A = 0 , R y = 1 ) . For fixed 0 < | x | < 1 , we have ( q = y R y = x , r = | y | R y = | x | ) . Thus C ( x ) = | K | r | 1 q | = | x | | Γ ( λ ) | | 1 + x | . Hence, by Lemma 2, | R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N . Consequently, | R N ( x ) | | x | | Γ ( λ ) | | 1 + x | N λ 1 | x | N . This is a linear/geometric tail with algebraic factor N λ 1 . Applying Lemma 1, we obtain ρ f λ , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 9 
(Polylogarithm series). The polylogarithm function has the power-series expansion Li s ( x ) = n = 1 x n n s , s > 0 , | x | < 1 . Here the summation starts from n = 1 , since n s is not defined at n = 0 . The coefficient sequence satisfies c n = 1 n s = n s . Therefore this is of the form in Lemma 2 with ( K = 1 , α = s , β = 0 , A = 0 , R = 1 ) . For fixed 0 < | x | < 1 , we have ( q = x R = x , r = | x | R = | x | ) . Thus C ( x ) = | K | r | 1 q | = | x | | 1 x | . Hence, by Lemma 2, | R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N . Consequently, | R N ( x ) | | x | | 1 x | N s | x | N . This is a linear/geometric tail with algebraic factor N s . Applying Lemma 1, we obtain ρ Li s , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 10 
(Catalan generating function series). The Catalan generating function has the power-series expansion f ( x ) = 1 1 4 x 2 x = n = 0 2 n n n + 1 x n , | x | < 1 4 . The coefficient sequence is c n = 2 n n n + 1 [17]. Using the central-binomial asymptotic 2 n n 4 n π n , we obtain c n = 2 n n n + 1 4 n π n ( n + 1 ) 1 π n 3 / 2 4 n . Therefore this is of the form in Lemma 2 with ( K = 1 π , α = 3 2 , β = 0 , A = 0 , R = 1 4 ) . For fixed 0 < | x | < 1 / 4 , we have ( q = x R = 4 x , r = | x | R = 4 | x | ) . Thus C ( x ) = | K | r | 1 q | = 4 | x | π | 1 4 x | . Hence, by Lemma 2, | R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N . Consequently, | R N ( x ) | 4 | x | π | 1 4 x | N 3 / 2 ( 4 | x | ) N . This is a linear/geometric tail with algebraic factor N 3 / 2 . Applying Lemma 1, we obtain ρ f , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / ( 4 | x | ) ) ( ε 0 ) .
Example 11 
(Lambert W 0 series). The principal Lambert W-function has the power-series expansion W 0 ( x ) = n = 1 ( n ) n 1 n ! x n , | x | < 1 e [18]. Since ( n ) n 1 n ! = ( 1 ) n 1 n n 1 n ! , we write y = x . Then W 0 ( x ) = n = 1 n n 1 n ! y n . Using Stirling’s approximation we obtain n n 1 n ! n n 1 2 π n ( n / e ) n = e n 2 π n 3 / 2 . Thus the coefficient sequence in the variable y satisfies c n = n n 1 n ! 1 2 π n 3 / 2 e n . Therefore this is of the form in Lemma 2 with ( K = 1 2 π , α = 3 2 , β = 0 , A = 0 , R y = 1 e ) . For fixed 0 < | x | < 1 / e , we have ( q = y R y = e x , r = | y | R y = e | x | ) . Thus C ( x ) = | K | r | 1 q | = e | x | 2 π | 1 + e x | . Hence, by Lemma 2, | R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N . Consequently, | R N ( x ) | e | x | 2 π | 1 + e x | N 3 / 2 ( e | x | ) N . This is a linear/geometric tail with algebraic factor N 3 / 2 . Applying Lemma 1, we obtain ρ W 0 , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / ( e | x | ) ) ( ε 0 ) .
Example 12 
(Complete elliptic integral series). The complete elliptic integral of the first kind has the power-series expansion K ( x ) = π 2 n = 0 ( 2 n n 4 n ) 2 x 2 n , | x | < 1 . Fix 0 < | x | < 1 , and write y = x 2 . Then K ( x ) = π 2 n = 0 ( 2 n n 4 n ) 2 y n . Using the central-binomial asymptotic 2 n n 4 n π n , we obtain ( 2 n n 4 n ) 2 1 π n . Thus, before multiplying by the outer factor π / 2 , the coefficient sequence in the variable y satisfies c n = ( 2 n n 4 n ) 2 1 π n 1 . Therefore this is of the form in Lemma 2 with ( K = 1 π , α = 1 , β = 0 , A = 0 , R y = 1 . ) For fixed 0 < | x | < 1 , we have ( q = y R y = x 2 , r = | y | R y = | x | 2 ) . Thus C y ( x ) = | K | r | 1 q | = | x | 2 π | 1 x 2 | . After multiplying by the outer factor π / 2 , Lemma 2 gives | R N ( x ) | π 2 C y ( x ) N 1 r N . Hence | R N ( x ) | | x | 2 2 | 1 x 2 | N 1 | x | 2 N . This is a lacunary linear/geometric tail with algebraic factor N 1 . Applying Lemma 1, we obtain ρ K , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | 2 ) = 1 2 log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 13 
(Secant function series). The secant function has the power-series expansion sec x = n = 0 | E 2 n | ( 2 n ) ! x 2 n , | x | < π 2 , where E 2 n denotes the Euler numbers. Using the Euler-number asymptotic Theorem 5 we obtain | E 2 n | ( 2 n ) ! 4 π 4 π 2 n . Fix 0 < | x | < π / 2 , and write y = x 2 . Then sec x = n = 0 | E 2 n | ( 2 n ) ! y n . Thus the coefficient sequence in the variable y satisfies c n = | E 2 n | ( 2 n ) ! 4 π 4 π 2 n . Therefore this is of the form in Lemma 2 with ( K = 4 π , α = 0 , β = 0 , A = 0 , R y = π 2 4 ) . For fixed 0 < | x | < π / 2 , we have ( q = y R y = 4 x 2 π 2 , r = | y | R y = 4 | x | 2 π 2 ) . Thus C y ( x ) = | K | r | 1 q | = 16 | x | 2 π 3 1 4 x 2 π 2 . Hence, by Lemma 2, | R N ( x ) | C y ( x ) N α ( log N ) β exp { A N } r N . Consequently, | R N ( x ) | 16 | x | 2 π 3 1 4 x 2 π 2 4 | x | 2 π 2 N . This is a lacunary linear/geometric tail. Applying Lemma 1, we obtain ρ sec , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log π 2 4 | x | 2 ( ε 0 ) .
Example 14 
(Tangent function series). The tangent function has the power-series expansion tan x = n = 1 B 2 n ( 4 ) n ( 1 4 n ) ( 2 n ) ! x 2 n 1 , | x | < π 2 . The n = 0 term is zero, so the effective summation begins at n = 1 . By Theorem 6, B 2 n = ( 1 ) n 1 2 ( 2 n ) ! ( 2 π ) 2 n ζ ( 2 n ) , ζ ( 2 n ) 1 . Therefore B 2 n ( 4 ) n ( 1 4 n ) ( 2 n ) ! 2 · 4 2 n ( 2 π ) 2 n = 2 4 π 2 n . Fix 0 < | x | < π / 2 , and write y = x 2 . Reindexing by m = n 1 , we have tan x = x m = 0 B 2 m + 2 ( 4 ) m + 1 ( 1 4 m + 1 ) ( 2 m + 2 ) ! y m . The coefficient sequence in the variable y satisfies c m = B 2 m + 2 ( 4 ) m + 1 ( 1 4 m + 1 ) ( 2 m + 2 ) ! 2 4 π 2 m + 1 = 8 π 2 4 π 2 m . Thus this is of the form in Lemma 2 with ( K = 8 π 2 , α = 0 , β = 0 , A = 0 , R y = π 2 4 ) . For fixed 0 < | x | < π / 2 , we have ( q = y R y = 4 x 2 π 2 , r = | y | R y = 4 | x | 2 π 2 ) . Thus C y ( x ) = | K | r | 1 q | = 32 | x | 2 π 4 1 4 x 2 π 2 . After multiplying by the outer factor x, Lemma 2 gives | R N ( x ) | | x | C y ( x ) r N . Consequently, | R N ( x ) | 32 | x | 3 π 4 1 4 x 2 π 2 4 | x | 2 π 2 N . This is a lacunary linear/geometric tail. Applying Lemma 1, we obtain ρ tan , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log π 2 4 | x | 2 ( ε 0 ) .
Example 15 
(Partition generating function series). The partition generating function has the power-series expansion f ( x ) = n = 1 p ( n ) x n , | x | < 1 , where p ( n ) denotes the number of partitions of n [15,19]. By Theorem 7, p ( n ) 1 4 3 n exp π 2 n 3 . Therefore the coefficient sequence satisfies c n = p ( n ) 1 4 3 n 1 exp π 2 3 n . This is of the form in Lemma 2 with ( K = 1 4 3 , α = 1 , β = 0 , A = π 2 3 , R = 1 ) . For fixed 0 < | x | < 1 , we have ( q = x R = x , r = | x | R = | x | ) . Thus C ( x ) = | K | r | 1 q | = | x | 4 3 | 1 x | . Hence, by Lemma 2, | R N ( x ) | C ( x ) N α ( log N ) β exp { A N } r N . Consequently, | R N ( x ) | | x | 4 3 | 1 x | N 1 exp π 2 N 3 | x | N . This is a linear/geometric tail with a subexponential correction. Applying Lemma 1, we obtain ρ f , 0 ( x ; ε ) log ( 1 / ε ) log ( 1 / r ) = log ( 1 / ε ) log ( 1 / | x | ) ( ε 0 ) .
Example 16 
(Jacobi theta-type lacunary series). The Jacobi theta-type function has the lacunary power-series expansion θ s ( x ) = 1 + 2 n = 1 s n x n 2 , | s | = 1 , | x | < 1 [6,14]. Here the nonzero powers occur only at the square indices n 2 . Therefore the series is not directly of the coefficient-asymptotic form in Lemma 2, since the ordinary coefficient sequence has many zeros. Instead, the tail is computed from the square-lacunary structure. For fixed 0 < | x | < 1 , define the lacunary-index tail by R N ( x ) = 2 n = N + 1 s n x n 2 . Since | s | = 1 , the absolute size of the terms is governed by | s n x n 2 | = | x | n 2 . Moreover, | x | ( n + 1 ) 2 | x | n 2 = | x | ( n + 1 ) 2 n 2 = | x | 2 n + 1 0 . Thus the terms decrease faster than geometrically in the lacunary index. Consequently, the tail is asymptotic to its first omitted term: | R N ( x ) | 2 | x | ( N + 1 ) 2 . Equivalently, | R N ( x ) | 2 | x | | x | N 2 + 2 N . The lacunary-index radius is therefore determined by 2 | x | ( N + 1 ) 2 < ε . Taking logarithms gives log 2 ( N + 1 ) 2 log ( 1 / | x | ) < log ( 1 / ε ) . Hence ( N + 1 ) 2 log ( 2 / ε ) log ( 1 / | x | ) . Therefore ρ θ s , 0 ( x ; ε ) log ( 2 / ε ) log ( 1 / | x | ) 1 / 2 ( ε 0 ) .

3.4.2. Superlinear/Entire Series: Class 3

The series discussed in this section have infinite classical radius R = . These series are entire and their tail is superlinear, usually factorial or stronger [11,14]. We start with the second key pattern lemma:
Lemma 3 
(Series Tail Pattern(2)). Let f ( x ) = n = 0 + c n x n be a real power series with infinite classical convergence radius R = + . Suppose, its tail satisfy:
| R N ( x ) | C ( x ) N α Λ ( x ) N N μ N ( N ) ,
where C ( x ) > 0 , Λ ( x ) > 0 , α R , and μ > 0 . Then:
ρ f , 0 ( x ; ε ) 1 μ log ( 1 / ε ) log log ( 1 / ε ) , ( ϵ 0 ) .
Proof. 
Set L : = log ( 1 / ε ) . Taking logarithms in the assumed tail asymptotic gives
log | R N ( x ) | = μ N log N + N log Λ ( x ) + α log N + log C ( x ) + o ( 1 ) .
Since N = o ( N log N ) , and log N = o ( N log N ) , the dominant term is μ N log N . Thus the inequality | R N ( x ) | < ε is governed asymptotically by
μ N log N L .
We now invert this relation. Fix δ ( 0 , 1 ) , and define N + ( ε ) : = ( 1 + δ ) L μ log L . At this scale, μ N + log N + = ( 1 + δ ) L ( 1 + o ( 1 ) ) . Moreover, N + log Λ ( x ) = o ( L ) , α log N + = o ( L ) , and log C ( x ) = o ( L ) . Therefore, log | R N + ( x ) | = ( 1 + δ ) L + o ( L ) < L for sufficiently small ε . Hence | R N + ( x ) | < e L = ε . Since the model tail C ( x ) N α Λ ( x ) N N μ N is eventually decreasing, N + ( ε ) is an admissible entry index. Consequently,
ρ f , 0 ( x ; ε ) ( 1 + δ ) L μ log L + O ( 1 ) .
For the reverse estimate, define N ( ε ) : = ( 1 δ ) L μ log L . Then μ N log N = ( 1 δ ) L ( 1 + o ( 1 ) ) , and again the remaining terms are o ( L ) . Hence log | R N ( x ) | = ( 1 δ ) L + o ( L ) > L for sufficiently small ε . Therefore | R N ( x ) | > ε , so no index below this scale can be admissible. Thus
ρ f , 0 ( x ; ε ) ( 1 δ ) L μ log L + O ( 1 ) .
Combining the two estimates in (37) and (38) gives
1 δ lim inf ε 0 ρ f , 0 ( x ; ε ) L / ( μ log L ) lim sup ε 0 ρ f , 0 ( x ; ε ) L / ( μ log L ) 1 + δ .
Letting δ 0 in (39), we obtain
ρ f , 0 ( x ; ε ) L μ log L = 1 μ log ( 1 / ε ) log log ( 1 / ε ) .
This proves the result. □
Lemma 4 
(Tail pattern for factorial-power entire series). Let f ( x ) = n = 0 u 1 ( n ) u 2 ( n ) x a n + b ( c n + d ) ! e , x R , where a , c , e N , and b , d N 0 . Assume that u 1 ( n ) u 2 ( n ) K θ n n q , where K > 0 , θ > 0 , q R . Then, for fixed x 0 , as N , we have
| R N ( x ) | C ( x ) N α Λ ( x ) N N μ N : μ = e c , Λ ( x ) = e e c θ | x | a c e c , α = q e c + d + 1 2 , C ( x ) = K θ | x | a + b ( 2 π ) e / 2 c e ( c + d + 1 / 2 ) .
Proof. 
The proof is constituted of several steps as follows:
Step 1. Write
R N ( x ) = n = N + 1 v n ( x ) : v n ( x ) : = u 1 ( n ) u 2 ( n ) x a n + b ( c n + d ) ! e .
By assumption, u 1 ( n ) u 2 ( n ) K θ n n q . Hence,
| v n ( x ) | K θ n n q | x | a n + b ( c n + d ) ! e .
Step 2. We first show that the tail R N ( x ) is asymptotic to its first omitted term. For x 0 , by (42) we have:
| v n + 1 ( x ) | | v n ( x ) | θ | x | a n + 1 n q ( c n + d ) ! ( c ( n + 1 ) + d ) ! e = θ | x | a n + 1 n q ( c n + d ) ! ( c n + d + c ) ! e = θ | x | a n + 1 n q 1 k = 1 c ( c n + d + k ) e θ | x | a n + 1 n q 1 c c n c e = θ | x | a n + 1 n q 1 c c e n c e θ | x | a c c e n c e 0 .
Consequently, by (43) the terms decrease faster than geometrically, and therefore by (41):
| R N ( x ) | | v N + 1 ( x ) | .
Step 3. It remains to compute the asymptotic form of v n ( x ) . By approximation (42) and Stirling’s formula,
| v n ( x ) | K θ n n q | x | a n + b ( c n + d ) ! e = K θ n n q | x | a n + b 2 π c c n + d + 1 / 2 n c n + d + 1 / 2 e c n e K θ n n q | x | a n + b e e c n ( 2 π ) e / 2 c e ( c n + d + 1 / 2 ) n e ( c n + d + 1 / 2 ) = K | x | b ( 2 π ) e / 2 c e ( d + 1 / 2 ) n q e ( d + 1 / 2 ) e e c θ | x | a c e c n n e c n .
Step 4. By (44) and (45) we have:
| R N ( x ) | K | x | b ( 2 π ) e / 2 c e ( d + 1 / 2 ) ( N + 1 ) q e ( d + 1 / 2 ) Λ ( x ) N + 1 ( N + 1 ) e c ( N + 1 ) : Λ ( x ) : = e e c θ | x | a c e c K | x | b ( 2 π ) e / 2 c e ( d + 1 / 2 ) ( N q e ( d + 1 / 2 ) ) Λ ( x ) N + 1 ( e e c N e c N e c N ) K | x | b Λ ( x ) e e c ( 2 π ) e / 2 c e ( d + 1 / 2 ) N q e ( d + 1 / 2 ) e c Λ ( x ) N N e c N : Λ ( x ) e e c = θ | x | a c e c , = K θ | x | a + b ( 2 π ) e / 2 c e ( c + d + 1 / 2 ) N q e ( c + d + 1 / 2 ) Λ ( x ) N N e c N .
This proves the stated tail-format estimate. □
Example 17 
( e w x ). The scaled exponential function has the power-series expansion e w x = n = 0 w n x n n ! . This is of the form in Lemma 4 with ( a = 1 , b = 0 , c = 1 , d = 0 , e = 1 , K = 1 , θ = w , q = 0 ) . Therefore ( μ = e c = 1 , Λ ( x ) = e | x | | w | , α = 3 2 , C ( x ) = | x | 2 π | w | ) . Hence, for fixed x 0 , we have | R N ( x ) | | x w | 2 π N 3 / 2 ( e | x w | ) N N N . By Lemma 3, ρ e x , 0 ( x ; ε ) log ( 1 / ε ) log log ( 1 / ε ) ( ε 0 ) .
Example 18 
( sin x ). The sine function has the power-series expansion sin x = n = 0 ( 1 ) n x 2 n + 1 ( 2 n + 1 ) ! . This corresponds to ( a = 2 , b = 1 , c = 2 , d = 1 , e = 1 , K = 1 , θ = 1 , q = 0 ) . Therefore ( μ = e c = 2 , Λ ( x ) = e 2 | x | 2 4 , α = 7 2 , C ( x ) = | x | 3 2 π 2 7 / 2 . ) Hence | R N ( x ) | | x | 3 2 π 2 7 / 2 N 7 / 2 e 2 | x | 2 4 N N 2 N . By Lemma 3, ρ sin , 0 ( x ; ε ) 1 2 log ( 1 / ε ) log log ( 1 / ε ) ( ε 0 ) .
Example 19 
( cos x ). The cosine function has the power-series expansion cos x = n = 0 ( 1 ) n x 2 n ( 2 n ) ! . This corresponds to ( a = 2 , b = 0 , c = 2 , d = 0 , e = 1 , K = 1 , θ = 1 , q = 0 ) . Therefore ( μ = e c = 2 , Λ ( x ) = e 2 | x | 2 4 , α = 5 2 , C ( x ) = | x | 2 2 π 2 5 / 2 ) . Hence | R N ( x ) | | x | 2 2 π 2 5 / 2 N 5 / 2 e 2 | x | 2 4 N N 2 N . By Lemma 3, ρ cos , 0 ( x ; ε ) 1 2 log ( 1 / ε ) log log ( 1 / ε ) ( ε 0 ) .
Example 20 
( sinh x ). The hyperbolic sine function has the power-series expansion sinh x = n = 0 x 2 n + 1 ( 2 n + 1 ) ! . This corresponds to ( a = 2 , b = 1 , c = 2 , d = 1 , e = 1 , K = 1 , θ = 1 , q = 0 ) . Therefore ( μ = e c = 2 , Λ ( x ) = e 2 | x | 2 4 , α = 7 2 , C ( x ) = | x | 3 2 π 2 7 / 2 ) . Hence | R N ( x ) | | x | 3 2 π 2 7 / 2 N 7 / 2 e 2 | x | 2 4 N N 2 N . By Lemma 3, ρ sinh , 0 ( x ; ε ) 1 2 log ( 1 / ε ) log log ( 1 / ε ) ( ε 0 ) .
Example 21 
( cosh x ). The hyperbolic cosine function has the power-series expansion cosh x = n = 0 x 2 n ( 2 n ) ! . This corresponds to ( a = 2 , b = 0 , c = 2 , d = 0 , e = 1 , K = 1 , θ = 1 , q = 0 ) . Therefore ( μ = e c = 2 , Λ ( x ) = e 2 | x | 2 4 , α = 5 2 , C ( x ) = | x | 2 2 π 2 5 / 2 ) . Hence | R N ( x ) | | x | 2 2 π 2 5 / 2 N 5 / 2 e 2 | x | 2 4 N N 2 N . By Lemma 3, ρ cosh , 0 ( x ; ε ) 1 2 log ( 1 / ε ) log log ( 1 / ε ) ( ε 0 ) .
Example 22 
( erf ( x ) ). The error function has the power-series expansion erf ( x ) = 2 π n = 0 ( 1 ) n x 2 n + 1 n ! ( 2 n + 1 ) . This corresponds to ( a = 2 , b = 1 , c = 1 , d = 0 , e = 1 ) . Moreover, u 1 ( n ) u 2 ( n ) = 2 π 1 2 n + 1 1 π n 1 . Thus ( K = 1 π , θ = 1 , q = 1 ) . Therefore ( μ = e c = 1 , Λ ( x ) = e | x | 2 , α = 5 2 , C ( x ) = | x | 3 π 2 ) . Hence | R N ( x ) | | x | 3 π 2 N 5 / 2 ( e | x | 2 ) N N N . By Lemma 3, ρ erf , 0 ( x ; ε ) log ( 1 / ε ) log log ( 1 / ε ) ( ε 0 ) .
Example 23 
( J 0 ( x ) ). The Bessel function J 0 has the power-series expansion J 0 ( x ) = n = 0 ( 1 ) n x 2 n 4 n ( n ! ) 2 . This corresponds to ( a = 2 , b = 0 , c = 1 , d = 0 , e = 2 ) . Moreover, u 1 ( n ) u 2 ( n ) = 4 n . Thus ( K = 1 , θ = 1 4 , q = 0 . ) Therefore ( μ = e c = 2 , Λ ( x ) = e 2 | x | 2 4 , α = 3 , C ( x ) = | x | 2 8 π ) . Hence | R N ( x ) | | x | 2 8 π N 3 e 2 | x | 2 4 N N 2 N . By Lemma 3, ρ J 0 , 0 ( x ; ε ) 1 2 log ( 1 / ε ) log log ( 1 / ε ) ( ε 0 ) .

3.4.3. Catalog Summary Table

Table 4 presents the summary of examples presented in previous sections. We identify the patterns as follows:
3.4.3.1. Classical-radius ranges.
The catalog separates the examples into finite-radius and infinite-radius families. Within the linear/geometric class, the classical radius may be smaller than one, equal to one, or larger than one: the Catalan and Lambert W 0 series have R < 1 , most logarithmic, polylogarithmic, binomial, lacunary, elliptic, partition, and theta-type examples have R = 1 , while the trigonometric tangent and secant examples have R = π / 2 > 1 . In contrast, the factorial-type examples in the superlinear class all have R = + . Thus the classical radius records the domain of convergence, but it does not by itself determine the sequential cost; the same finite radius range may contain several different coefficient structures and lacunarity patterns.
3.4.3.2. Speed of the indexed radius.
For the linear/geometric examples, the indexed radius grows on the logarithmic scale ρ ( x ; ε ) log ( 1 / ε ) log ( R / | x | ) , up to constant factors caused by lacunarity, such as 1 / b or 1 / 2 . Algebraic, logarithmic, and subexponential coefficient corrections, such as N α , log N , or exp ( A N ) , affect the tail amplitude and lower-order behavior but do not change the leading logarithmic order. By contrast, the superlinear entire examples have factorial-type tails and therefore require only ρ ( x ; ε ) log ( 1 / ε ) log log ( 1 / ε ) , with constants such as 1 / μ distinguishing full factorial, even/odd factorial subsequences, and squared-factorial structures.
3.4.3.3. Overall findings.
The main lesson of the catalog is that the sequential radius refines the classical radius by measuring the approximation cost inside the convergence interval. Finite-radius series in the linear class can have very different coefficient origins—algebraic singularities, central-binomial coefficients, Bernoulli or Euler numbers, partition asymptotics, or lacunary supports—yet their leading indexed-radius behavior is governed by the same geometric mechanism. Entire factorial-type series form a separate, faster class: their infinite classical radius is accompanied by substantially smaller indexed radii. Hence the table shows that R tells where the power series converges, while ρ describes how many terms are needed to converge at a prescribed accuracy.
0pt 0pt [longtable]width=
The preceding catalog shows that many classical examples fall into two tractable model families: geometric tails with lower-order corrections and factorial-type entire tails. For later tail-transfer computations, it is useful to isolate these two families as explicit subclasses of the broader linear and superlinear classes. The following definition records the subclasses on which the subsequent transformation analysis is intended to operate.
Definition 11 
(Model subclasses for tail-transfer computations). Let P 0 be the space of real power series centered at 0, as in Definition 4. For f P 0 , write f ( x ) = n = 0 c n x n , s N f ( x ) : = n = 0 N c n x n , and R N f ( x ) : = f ( x ) s N f ( x ) . The tail-transfer computations below are carried out on the following two model subclasses of the tail classes introduced in Definition 10:
1 . Model Linear class . P 0 , mod lin ( x ) : = α , β , A R R > 0 f P 0 : | R N f ( x ) | C f ( x ) N α ( log N ) β exp ( A N ) r f , x N . P 0 , coef lin ( x ) : = K 0 , α , β , A R R > 0 f P 0 : f ( x ) = n = 0 c n x n , c n K n α ( log n ) β exp ( A n ) R n . 2 . Model Superlinear class . P 0 , mod sup ( x ) : = μ > 0 α R f P 0 : | R N f ( x ) | C f ( x ) N α Λ f ( x ) N N μ N . P 0 , fac sup ( x ) : = a , c , e N , b , d N 0 K > 0 , θ > 0 , q R f P 0 : f ( x ) = n = 0 u 1 ( n ) u 2 ( n ) x a n + b ( ( c n + d ) ! ) e , u 1 ( n ) u 2 ( n ) K θ n n q .
Remark 7 (Class Subsets).
The following chain of subsets hold using Lemma 1, Lemma 2, Lemma 3, Lemma 4 and Definition 10:
P 0 , coef lin ( x ) P 0 , mod lin ( x ) P 0 lin ( x ) P 0 seq ( x ) P 0 .
P 0 , fac sup ( x ) P 0 , mod sup ( x ) P 0 sup ( x ) P 0 seq ( x ) P 0 .
Figure 3 presents the Venn diagram of the discussed sets of series in inequalities (47) and () in this section:

4. Discussion

The preceding sections introduced sequential radii for power series, established their basic classification, and computed them for a range of standard examples. We now summarize the main contributions of the paper, indicate the limitations of the present framework, and describe several directions in which the theory may be extended.

4.1. Summary and Contributions

This paper contributes to the theory of powers series in three dimensions as follows:
The first and the most prominent contribution of this paper is to complement the classical radius with the sequential radii ρ f , a ( x ; ε ) and ρ f , a * ( x ; ε ) , which measure how many terms are needed to approximate the value of the series within a prescribed tolerance. In this sense, the sequential radius records a quantitative convergence profile rather than only a qualitative convergence domain. The classification theorem shows that the asymptotic behavior of the tail remainder R n ( x ) induces three basic convergence-cost regimes. Sublinear tails lead to polynomial indexed radii, linear or geometric tails lead to logarithmic indexed radii, and superlinear or factorial tails lead to still faster logarithmic-over-logarithmic growth. Thus, the paper separates power series not only according to where they converge, but also according to how efficiently their partial sums approach the limiting value.
The second contribution is the geometric interpretation of the linear regime. In this case, the effective ratio r x = | x a | / R shows explicitly how the evaluation point x affects the approximation cost. The indexed radius is small (i.e., low cost computationally) near the center of expansion and increases (i.e., higher cost computationally) as x approaches the boundary of the classical interval of convergence. This gives a geometric refinement of the classical radius: points inside the same convergence interval may have substantially different computational costs.
Finally, the paper records stability properties of the sequential radii under elementary algebraic operations and provides a catalog of classical examples. These examples show that many familiar power series fall naturally into the linear or superlinear tail classes. The resulting computations illustrate that the proposed radii are not merely formal definitions, but can be explicitly evaluated for standard analytic functions.
The results above show that sequential radii provide a useful quantitative refinement of the classical radius. At the same time, the present treatment is intentionally model-based, and its scope is limited by the tail regimes and pointwise setting considered here.

4.2. Limitations

The limitations of this work are clear. First, the present work focuses on model tail regimes. The three classes considered here capture several important patterns, but they do not exhaust all possible tail behaviors of convergent power series. More irregular tails, oscillatory remainders, mixed asymptotic scales, or slowly varying corrections may require refinements beyond the three basic regimes studied in this paper.
The second limitation is that most explicit computations are carried out under asymptotic assumptions strong enough to permit direct tail inversion. In particular, the catalog examples rely on tails that can be reduced to geometric, algebraic-geometric, or factorial-type forms. For more general coefficient structures, the existence of the classical radius does not automatically provide enough information to determine the sequential radius.
The third limitation is that the algebraic stability results should also be understood as bounds rather than exact formulas. For sums, products, and quotients, the sequential radius of the resulting expression may depend on cancellations, coefficient interactions, or the particular approximating sequence used. This is especially visible in the distinction among pointwise, Cauchy, and Hadamard products, where identical individual radii need not determine the product radius.
Finally, the paper is restricted to real power series at a fixed center and to pointwise evaluation inside the interval of convergence. there is no guarantee that the discussed properties of the sequential radii of series remain invariant outside the scope of the present treatment.
The limitations above should be read as boundaries of the present formulation rather than obstructions to the approach. They identify the points at which the current pointwise and model-tail theory would need to be refined in order to support broader versions of the sequential-radius framework.

4.3. Future Work

This work has several natural continuations. First, extensions to complex power series, and multivariable power series is plausible. Second, extensions to uniform radii over compact subsets, and endpoint behavior is also plausible. Finally, a natural continuation of this work is to develop a systematic theory of tail-transfer principles for transformed series [9,11,20,21]. Here, given a power series A ( x ) = n = 0 a n x n and a transformation h = ( ( f n ) n = 0 , g ) , one may define a transformed series T h ( A ( x ) ) = B ( x ) : = n = 0 f n ( a n ) g ( x n ) . The central question is not only how T h changes the classical radius of convergence, but how it transfers the original series into a new convergence-cost class. Future work may therefore classify transformations according to their effect on both the classical radius R and the sequential radii ρ and ρ * . Such a theory would connect coefficient transformations, argument transformations, and tail asymptotics in a unified framework. It would also help explain how operations on a power series modify not only its domain of convergence, but also the computational cost of approximating its values.

Conclusions

The sequential-radius viewpoint developed in this paper gives a quantitative language for comparing power-series approximations after convergence has already been established. Its value lies in separating the existence of a limiting value from the cost of reaching it by partial sums. In this form, the classical radius remains the geometric boundary of convergence, while the sequential radii describe the internal rate structure of the convergent region. This distinction provides a basis for further study of power series through their tails, transformations, and approximation efficiency.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Hairer, E.; Wanner, G. Analysis by its history; Springer, 2008. [Google Scholar]
  2. Grabiner, J. V. The origins of Cauchy’s rigorous calculus; MIT Press, 1981. [Google Scholar]
  3. Stillwell, J. Mathematics and its history, 3rd ed.; Springer, 2010. [Google Scholar]
  4. Knopp, K. Theory and application of infinite series; Dover Publications, 1990. [Google Scholar]
  5. Rudin, W. Principles of mathematical analysis, 3rd ed.; McGraw–Hill, 1976. [Google Scholar]
  6. Whittaker, E. T.; Watson, G. N. A course of modern analysis; Cambridge University Press, 1996. [Google Scholar]
  7. Abel, N. H. Untersuchungen über die Reihe 1+m1x+m(m-1)1·2x2+m(m-1)(m-2)1·2·3x3+…. J. Für Die Reine Und Angew. Math. 1826, 1, 311–339. [Google Scholar]
  8. Ahlfors, L. V. Complex analysis, 3rd ed.; McGraw–Hill, 1979. [Google Scholar]
  9. Conway, J. B. Functions of one complex variable I, 2nd ed.; Springer, 1978. [Google Scholar]
  10. Soltanifar, M. What is the radius of convergence in the sequence space Seq(R)? Preprints 2025. [Google Scholar] [CrossRef]
  11. Boas, R. P. Entire functions; Academic Press, 1954. [Google Scholar]
  12. Olver, F. W. J. Asymptotics and special functions; A K Peters, 1997. [Google Scholar]
  13. Olver, F. W. J.; Lozier, D. W.; Boisvert, R. F.; Clark, C. W. (Eds.) NIST handbook of mathematical functions; Cambridge University Press, 2010. [Google Scholar]
  14. NIST Digital Library of Mathematical Functions. NIST digital library of mathematical functions; Olver, F. W. J., Olde Daalhuis, A. B., 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.; Release 1.2.1 of 2024-06-15, 2024; Available online: https://dlmf.nist.gov/.
  15. Andrews, G. E. The theory of partitions; Cambridge University Press, 1984. [Google Scholar]
  16. Bingham, N. H.; Goldie, C. M.; Teugels, J. L. Regular variation; Cambridge University Press, 1987. [Google Scholar]
  17. Flajolet, P.; Sedgewick, R. Analytic combinatorics; Cambridge University Press, 2009. [Google Scholar]
  18. Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; Jeffrey, D. J.; Knuth, D. E. On the Lambert (W) function. Adv. Comput. Math. 5 1996, 329–359. [Google Scholar] [CrossRef]
  19. Hardy, G. H.; Ramanujan, S. Asymptotic formulae in combinatory analysis. Proc. Lond. Math. Soc. 1918, 17(1), 75–115. [Google Scholar] [CrossRef]
  20. Silverman, R. A. Calculus with Analytic Geometry; Prentice-Hall, Inc.: New Jersey, USA, 1985; ISBN 0-13-111634-7. [Google Scholar]
  21. Henrici, P. Applied and computational complex analysis; Power series, integration, conformal mapping, location of zeros; Wiley, 1974; Volume 1. [Google Scholar]
Figure 1. Structural consequences of Theorem 9. Panel (a) illustrates the role of the tail constant C: changing C = C 2 for the linear tail regime (with fixed parameters ( R = 1 , | x a | = 0.5 ) ) vertically shifts the indexed-radius curve but does not change its asymptotic growth rate. Panel (b) compares the three sequential convergence classes induced by sublinear (with parameters: ( p = 2 , C 1 = 1 )), linear (with parameters ( R / | x a | = 1 , C 2 = 1 )), and superlinear (with parameters ( C 3 = 1 )) tail decay, showing the separation between polynomial, logarithmic, and log / log log -type indexed-radius growth. Panel (c) illustrates the spatial geometry of the linear regime (with parameters ( C 2 = 1 , R = 1 )), where the indexed radius grows rapidly as r = | x a | / R approaches 1. The informal shaded “Near Boundary Zone” denotes points close to the boundary of the classical interval of convergence, where 0.9 < r < 1 , equivalently | x a | R , and hence log ( R / | x a | ) becomes small, causing the number of required terms to become very large.
Figure 1. Structural consequences of Theorem 9. Panel (a) illustrates the role of the tail constant C: changing C = C 2 for the linear tail regime (with fixed parameters ( R = 1 , | x a | = 0.5 ) ) vertically shifts the indexed-radius curve but does not change its asymptotic growth rate. Panel (b) compares the three sequential convergence classes induced by sublinear (with parameters: ( p = 2 , C 1 = 1 )), linear (with parameters ( R / | x a | = 1 , C 2 = 1 )), and superlinear (with parameters ( C 3 = 1 )) tail decay, showing the separation between polynomial, logarithmic, and log / log log -type indexed-radius growth. Panel (c) illustrates the spatial geometry of the linear regime (with parameters ( C 2 = 1 , R = 1 )), where the indexed radius grows rapidly as r = | x a | / R approaches 1. The informal shaded “Near Boundary Zone” denotes points close to the boundary of the classical interval of convergence, where 0.9 < r < 1 , equivalently | x a | R , and hence log ( R / | x a | ) becomes small, causing the number of required terms to become very large.
Preprints 216520 g001
Figure 2. Two-dimensional geometry of the sequential radii in the linear regime. The parameters are C 2 = 1 , R = 1 , and a = 0 . Panel (a) shows the indexed radius ρ f , 0 ( x ; ε ) log ( 1 / ε ) / log ( 1 / | x | ) , illustrating the growth in the required number of terms as the point x approaches the boundary | x | = 1 . Panel (b) shows the associated geometric radius ρ f , 0 * ( x ; ε ) log ( 1 / | x | ) / log ( 1 / ε ) , which records the reciprocal efficiency scale and decreases toward zero near the boundary and under increasingly strict precision demands.
Figure 2. Two-dimensional geometry of the sequential radii in the linear regime. The parameters are C 2 = 1 , R = 1 , and a = 0 . Panel (a) shows the indexed radius ρ f , 0 ( x ; ε ) log ( 1 / ε ) / log ( 1 / | x | ) , illustrating the growth in the required number of terms as the point x approaches the boundary | x | = 1 . Panel (b) shows the associated geometric radius ρ f , 0 * ( x ; ε ) log ( 1 / | x | ) / log ( 1 / ε ) , which records the reciprocal efficiency scale and decreases toward zero near the boundary and under increasingly strict precision demands.
Preprints 216520 g002
Figure 3. The Venn diagram for the nested tail-regime subclasses inside the sequential-radius power-series classes for the space of Power series P 0 .
Figure 3. The Venn diagram for the nested tail-regime subclasses inside the sequential-radius power-series classes for the space of Power series P 0 .
Preprints 216520 g003
Table 1. Classical radius ranges associated with the three tail-decay regimes for power series.
Table 1. Classical radius ranges associated with the three tail-decay regimes for power series.
Class Tail regime Typical tail Classical Radius R Position x
Class 1 Sublinear | R n ( x ) | C 1 n ( p 1 ) , C 1 > 0 , p > 1 0 < R 1 < Usually | x a | = R 1
Class 2 Linear | R n ( x ) | C 2 r n , 0 < r < 1 0 < R 2 < Usually | x a | < R 2
Class 3 Superlinear | R n ( x ) | C 3 ( n + 1 ) ! , C 3 > 0 R 3 = Any fixed x R
Table 2. Tail-decay regimes and the induced radii ρ s ( ε ) and ρ s * ( ε ) .
Table 2. Tail-decay regimes and the induced radii ρ s ( ε ) and ρ s * ( ε ) .
# Tail regime Tail asymptotic Indexed-radius ρ s ( ε ) Geometric radius ρ s * ( ε )
1 Sublinear | R n | C 1 n ( p 1 ) , C 1 > 0 , p > 1 C 1 ε 1 / ( p 1 ) ε C 1 1 / ( p 1 )
2 Linear | R n | C 2 r n , C 2 > 0 , 0 < r < 1 log ( C 2 / ε ) log ( 1 / r ) log ( 1 / r ) log ( C 2 / ε )
3 Superlinear | R n | C 3 ( n + 1 ) ! , C 3 > 0 log ( C 3 / ε ) log log ( C 3 / ε ) log log ( C 3 / ε ) log ( C 3 / ε )
Table 4. Summary of catalog classical power series with linear/geometric/superlinear-type tails.
Table 4. Summary of catalog classical power series with linear/geometric/superlinear-type tails.
Function(formula) Series Classic R Indexed radius ρ Features
Function(formula) Series Classic R Indexed radius ρ Features
Continued on next page
Linear: Class 2
Lacunary geometric: x a 1 x b n = 0 x b n + a 1 1 b log ( 1 / ε ) log ( 1 / | x | ) Lacunary geometric; nonzero exponents b n + a .
Reciprocal square: 1 ( 1 x ) 2 n = 1 n x n 1 1 log ( 1 / ε ) log ( 1 / | x | ) Linear tail with algebraic factor N.
Logarithm: log ( 1 + x ) n = 1 ( 1 ) n 1 n x n 1 log ( 1 / ε ) log ( 1 / | x | ) Alternating geometric tail with algebraic factor N 1 .
Inverse hyperbolic tangent: arctanh ( x ) n = 0 x 2 n + 1 2 n + 1 1 1 2 log ( 1 / ε ) log ( 1 / | x | ) Odd lacunary logarithmic tail; coefficient factor N 1 .
Harmonic generating function: log ( 1 x ) 1 x n = 1 H n x n 1 log ( 1 / ε ) log ( 1 / | x | ) Linear tail with logarithmic correction log N .
Legendre chi: χ s ( x ) , s = 2 , 3 n = 0 x 2 n + 1 ( 2 n + 1 ) s 1 1 2 log ( 1 / ε ) log ( 1 / | x | ) Odd lacunary polylog-type tail; algebraic factor N s .
Arcsine: arcsin ( x ) n = 0 2 n n 4 n ( 2 n + 1 ) x 2 n + 1 1 1 2 log ( 1 / ε ) log ( 1 / | x | ) Odd lacunary; central-binomial coefficient gives N 3 / 2 .
Binomial: ( 1 + x ) λ , λ Z n = 0 λ n x n 1 log ( 1 / ε ) log ( 1 / | x | ) Algebraic singularity at x = 1 ; coefficient factor N λ 1 .
Polylogarithm: Li s ( x ) , s > 0 n = 1 x n n s 1 log ( 1 / ε ) log ( 1 / | x | ) Linear tail with algebraic factor N s .
Catalan generating function: 1 1 4 x 2 x n = 0 2 n n n + 1 x n 1 4 log ( 1 / ε ) log ( 1 / ( 4 | x | ) ) Finite radius 1 / 4 ; coefficient factor N 3 / 2 4 N .
Lambert W 0 : W 0 ( x ) n = 1 ( n ) n 1 n ! x n 1 e log ( 1 / ε ) log ( 1 / ( e | x | ) ) Finite radius 1 / e ; Stirling gives N 3 / 2 e N .
Complete elliptic integral K ( x ) π 2 n = 0 2 n n 4 n 2 x 2 n 1 1 2 log ( 1 / ε ) log ( 1 / | x | ) Even lacunary; coefficient factor N 1 .
Secant: sec x n = 0 | E 2 n | ( 2 n ) ! x 2 n π 2 log ( 1 / ε ) log ( π 2 / ( 4 | x | 2 ) ) Even lacunary; Euler-number asymptotic gives geometric factor ( 4 / π 2 ) N .
Tangent: tan x n = 1 B 2 n ( 4 ) n ( 1 4 n ) ( 2 n ) ! x 2 n 1 π 2 log ( 1 / ε ) log ( π 2 / ( 4 | x | 2 ) ) Odd lacunary; Bernoulli-number asymptotic gives geometric factor ( 4 / π 2 ) N .
Partition generating function: n = 1 1 1 x n n = 1 p ( n ) x n 1 log ( 1 / ε ) log ( 1 / | x | ) Linear tail with subexponential correction exp { π 2 N / 3 } .
Jacobi theta-type: θ s ( x ) 1 + 2 n = 1 s n x n 2 , | s | = 1 1 log ( 2 / ε ) log ( 1 / | x | ) 1 / 2 Square-lacunary; not directly Lemma 2; tail 2 | x | ( N + 1 ) 2 .
Superlinear: Class 3
Scaled exponential: e w x n = 0 w n x n n ! + log ( 1 / ε ) log log ( 1 / ε ) Full factorial; μ = 1 .
Sine: sin x n = 0 ( 1 ) n x 2 n + 1 ( 2 n + 1 ) ! + 1 2 log ( 1 / ε ) log log ( 1 / ε ) Odd factorial subsequence; alternating; μ = 2 .
Cosine: cos x n = 0 ( 1 ) n x 2 n ( 2 n ) ! + 1 2 log ( 1 / ε ) log log ( 1 / ε ) Even factorial subsequence; alternating; μ = 2 .
Sinh: sinh x n = 0 x 2 n + 1 ( 2 n + 1 ) ! + 1 2 log ( 1 / ε ) log log ( 1 / ε ) Odd factorial subsequence; positive signs; μ = 2 .
Cosh: cosh x n = 0 x 2 n ( 2 n ) ! + 1 2 log ( 1 / ε ) log log ( 1 / ε ) Even factorial subsequence; positive signs; μ = 2 .
Error: erf ( x ) 2 π n = 0 ( 1 ) n x 2 n + 1 n ! ( 2 n + 1 ) + log ( 1 / ε ) log log ( 1 / ε ) Factorial with polynomial denominator; μ = 1 .
Bessel: J 0 ( x ) n = 0 ( 1 ) n x 2 n 4 n ( n ! ) 2 + 1 2 log ( 1 / ε ) log log ( 1 / ε ) Squared factorial denominator; μ = 2 .
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.