Les séries divergentes sont en général quelque chose de bien fatal et c’est une honte qu’on ose y fonder aucune démonstration. ("Divergent series are in general something fatal, and it is a disgrace to base any proof on them." Often translated as "Divergent series are an invention of the devil …")

N. H. Abel, letter to Holmboe, January 1826, reprinted in volume 2 of his collected papers.

In mathematics, a **divergent series** is an infinite series that is not convergent, meaning that the infinite sequence of the partial sums of the series does not have a finite limit.

If a series converges, the individual terms of the series must approach zero. Thus any series in which the individual terms do not approach zero diverges. However, convergence is a stronger condition: not all series whose terms approach zero converge. A counterexample is the harmonic series

The divergence of the harmonic series was proven by the medieval mathematician Nicole Oresme.

In specialized mathematical contexts, values can be objectively assigned to certain series whose sequences of partial sums diverge, in order to make meaning of the divergence of the series. A *summability method* or *summation method* is a partial function from the set of series to values. For example, Cesàro summation assigns Grandi's divergent series

the value 1/2. Cesàro summation is an *averaging* method, in that it relies on the arithmetic mean of the sequence of partial sums. Other methods involve analytic continuations of related series. In physics, there are a wide variety of summability methods; these are discussed in greater detail in the article on regularization.

... but it is broadly true to say that mathematicians before Cauchy asked not 'How shall we

define1 − 1 + 1...?' but 'Whatis1 − 1 + 1...?', and that this habit of mind led them into unnecessary perplexities and controversies which were often really verbal.

G. H. Hardy, Divergent series, page 6

Before the 19th century, divergent series were widely used by Leonhard Euler and others, but often led to confusing and contradictory results. A major problem was Euler's idea that any divergent series should have a natural sum, without first defining what is meant by the sum of a divergent series. Augustin-Louis Cauchy eventually gave a rigorous definition of the sum of a (convergent) series, and for some time after this, divergent series were mostly excluded from mathematics. They reappeared in 1886 with Henri Poincaré's work on asymptotic series. In 1890, Ernesto Cesàro realized that one could give a rigorous definition of the sum of some divergent series, and defined Cesàro summation. (This was not the first use of Cesàro summation, which was used implicitly by Ferdinand Georg Frobenius in 1880; Cesàro's key contribution was not the discovery of this method, but his idea that one should give an explicit definition of the sum of a divergent series.) In the years after Cesàro's paper, several other mathematicians gave other definitions of the sum of a divergent series, although these are not always compatible: different definitions can give different answers for the sum of the same divergent series; so, when talking about the sum of a divergent series, it is necessary to specify which summation method one is using.

A summability method *M* is *regular* if it agrees with the actual limit on all convergent series. Such a result is called an *Abelian theorem* for *M*, from the prototypical Abel's theorem. More subtle, are partial converse results, called *Tauberian theorems*, from a prototype proved by Alfred Tauber. Here *partial converse* means that if *M* sums the series *Σ*, and some side-condition holds, then *Σ* was convergent in the first place; without any side-condition such a result would say that *M* only summed convergent series (making it useless as a summation method for divergent series).

The function giving the sum of a convergent series is *linear*, and it follows from the Hahn–Banach theorem that it may be extended to a summation method summing any series with bounded partial sums. This is called the *Banach limit*. This fact is not very useful in practice, since there are many such extensions, inconsistent with each other, and also since proving such operators exist requires invoking the axiom of choice or its equivalents, such as Zorn's lemma. They are therefore nonconstructive.

The subject of divergent series, as a domain of mathematical analysis, is primarily concerned with explicit and natural techniques such as Abel summation, Cesàro summation and Borel summation, and their relationships. The advent of Wiener's tauberian theorem marked an epoch in the subject, introducing unexpected connections to Banach algebra methods in Fourier analysis.

Summation of divergent series is also related to extrapolation methods and sequence transformations as numerical techniques. Examples of such techniques are Padé approximants, Levin-type sequence transformations, and order-dependent mappings related to renormalization techniques for large-order perturbation theory in quantum mechanics.

Summation methods usually concentrate on the sequence of partial sums of the series. While this sequence does not converge, we may often find that when we take an average of larger and larger numbers of initial terms of the sequence, the average converges, and we can use this average instead of a limit to evaluate the sum of the series. A *summation method* can be seen as a function from a set of sequences of partial sums to values. If **A** is any summation method assigning values to a set of sequences, we may mechanically translate this to a *series-summation method* **A**^{Σ} that assigns the same values to the corresponding series. There are certain properties it is desirable for these methods to possess if they are to arrive at values corresponding to limits and sums, respectively.

**Regularity**. A summation method is*regular*if, whenever the sequence*s*converges to*x*,**A**(*s*) =*x*. Equivalently, the corresponding series-summation method evaluates**A**^{Σ}(*a*) =*x*.**Linearity**.**A**is*linear*if it is a linear functional on the sequences where it is defined, so that**A**(*k**r*+*s*) =*k***A**(*r*) +**A**(*s*) for sequences*r*,*s*and a real or complex scalar*k*. Since the terms*a*_{n+1}=*s*_{n+1}−*s*_{n}of the series*a*are linear functionals on the sequence*s*and vice versa, this is equivalent to**A**^{Σ}being a linear functional on the terms of the series.**Stability**(also called*translativity*). If*s*is a sequence starting from*s*_{0}and*s*′ is the sequence obtained by omitting the first value and subtracting it from the rest, so that*s*′_{n}=*s*_{n+1}−*s*_{0}, then**A**(*s*) is defined if and only if**A**(*s*′) is defined, and**A**(*s*) =*s*_{0}+**A**(*s*′). Equivalently, whenever*a*′_{n}=*a*_{n+1}for all*n*, then**A**^{Σ}(*a*) =*a*_{0}+**A**^{Σ}(*a*′).^{[1]}^{[2]}Another way of stating this is that the shift rule must be valid for the series that are summable by this method.

The third condition is less important, and some significant methods, such as Borel summation, do not possess it.^{[3]}

One can also give a weaker alternative to the last condition.

**Finite re-indexability**. If*a*and*a*′ are two series such that there exists a bijection such that*a*_{i}=*a*′_{f(i)}for all*i*, and if there exists some such that*a*_{i}=*a*′_{i}for all*i*>*N*, then**A**^{Σ}(*a*) =**A**^{Σ}(*a*′). (In other words,*a*′ is the same series as*a*, with only finitely many terms re-indexed.) This is a weaker condition than*stability*, because any summation method that exhibits*stability*also exhibits*finite re-indexability*, but the converse is not true.)

A desirable property for two distinct summation methods **A** and **B** to share is *consistency*: **A** and **B** are consistent if for every sequence *s* to which both assign a value, **A**(*s*) = **B**(*s*). (Using this language, a summation method **A** is regular iff it is consistent with the standard sum **Σ**.) If two methods are consistent, and one sums more series than the other, the one summing more series is *stronger*.

There are powerful numerical summation methods that are neither regular nor linear, for instance nonlinear sequence transformations like Levin-type sequence transformations and Padé approximants, as well as the order-dependent mappings of perturbative series based on renormalization techniques.

Taking regularity, linearity and stability as axioms, it is possible to sum many divergent series by elementary algebraic manipulations. This partly explains why many different summation methods give the same answer for certain series.

For instance, whenever *r* ≠ 1, the geometric series

can be evaluated regardless of convergence. More rigorously, any summation method that possesses these properties and which assigns a finite value to the geometric series must assign this value. However, when *r* is a real number larger than 1, the partial sums increase without bound, and averaging methods assign a limit of infinity.

The two classical summation methods for series, ordinary convergence and absolute convergence, define the sum as a limit of certain partial sums. These are included only for completeness; strictly speaking they are not true summation methods for divergent series since, by definition, a series is divergent only if these methods do not work. Most but not all summation methods for divergent series extend these methods to a larger class of sequences.

Absolute convergence defines the sum of a sequence (or set) of numbers to be the limit of the net of all partial sums *a*_{k1} + ... + *a*_{kn}, if it exists. It does not depend on the order of the elements of the sequence, and a classical theorem says that a sequence is absolutely convergent if and only if the sequence of absolute values is convergent in the standard sense.

Cauchy's classical definition of the sum of a series *a*_{0} + *a*_{1} + ... defines the sum to be the limit of the sequence of partial sums *a*_{0} + ... + *a*_{n}. This is the default definition of convergence of a sequence.

Suppose *p _{n}* is a sequence of positive terms, starting from

If now we transform a sequence s by using *p* to give weighted means, setting

then the limit of *t _{n}* as

The Nørlund mean is regular, linear, and stable. Moreover, any two Nørlund means are consistent.

The most significant of the Nørlund means are the Cesàro sums. Here, if we define the sequence *p ^{k}* by

then the Cesàro sum *C*_{k} is defined by *C*_{k}(*s*) = **N**_{(pk)}(*s*). Cesàro sums are Nørlund means if *k* ≥ 0, and hence are regular, linear, stable, and consistent. *C*_{0} is ordinary summation, and *C*_{1} is ordinary Cesàro summation. Cesàro sums have the property that if *h* > *k*, then *C*_{h} is stronger than *C*_{k}.

Suppose *λ* = {*λ*_{0}, *λ*_{1}, *λ*_{2},...} is a strictly increasing sequence tending towards infinity, and that *λ*_{0} ≥ 0. Suppose

converges for all real numbers *x* > 0. Then the *Abelian mean* *A*_{λ} is defined as

More generally, if the series for *f* only converges for large *x* but can be analytically continued to all positive real *x*, then one can still define the sum of the divergent series by the limit above.

A series of this type is known as a generalized Dirichlet series; in applications to physics, this is known as the method of *heat-kernel regularization*.

Abelian means are regular and linear, but not stable and not always consistent between different choices of *λ*. However, some special cases are very important summation methods.

See also: Abel's theorem |

If *λ*_{n} = *n*, then we obtain the method of *Abel summation*. Here

where *z* = exp(−*x*). Then the limit of *f*(*x*) as *x* approaches 0 through positive reals is the limit of the power series for *f*(*z*) as *z* approaches 1 from below through positive reals, and the Abel sum *A*(*s*) is defined as

Abel summation is interesting in part because it is consistent with but more powerful than Cesàro summation: *A*(*s*) = *C*_{k}(*s*) whenever the latter is defined. The Abel sum is therefore regular, linear, stable, and consistent with Cesàro summation.

If *λ*_{n} = *n* log(*n*), then (indexing from one) we have

Then *L*(*s*), the *Lindelöf sum* (Volkov 2001), is the limit of *f*(*x*) as *x* goes to positive zero. The Lindelöf sum is a powerful method when applied to power series among other applications, summing power series in the Mittag-Leffler star.

If *g*(*z*) is analytic in a disk around zero, and hence has a Maclaurin series *G*(*z*) with a positive radius of convergence, then *L*(*G*(*z*)) = *g*(*z*) in the Mittag-Leffler star. Moreover, convergence to *g*(*z*) is uniform on compact subsets of the star.

Several summation methods involve taking the value of an analytic continuation of a function.

If Σ*a*_{n}*x*^{n} converges for small complex *x* and can be analytically continued along some path from *x* = 0 to the point *x* = 1, then the sum of the series can be defined to be the value at *x* = 1. This value may depend on the choice of path. One of the first examples of potentially different sums for a divergent series, using analytic continuation, was given by Callet,^{[4]} who observed that if then

Evaluating at , one gets

However, the gaps in the series are key. For for example, we actually would get

, so different sums correspond to different placements of the 's.

Another example of analytic continuation is the divergent alternating series which is a sum over products of -functions and Pochhammer's symbols. Using the duplication formula of the -function, it reduces to a generalized hypergeometric series

Main article: Euler summation |

Euler summation is essentially an explicit form of analytic continuation. If a power series converges for small complex *z* and can be analytically continued to the open disk with diameter from −1/*q* + 1 to 1 and is continuous at 1, then its value at *q* is called the Euler or (E,*q*) sum of the series Σ*a*_{n}. Euler used it before analytic continuation was defined in general, and gave explicit formulas for the power series of the analytic continuation.

The operation of Euler summation can be repeated several times, and this is essentially equivalent to taking an analytic continuation of a power series to the point *z* = 1.

This method defines the sum of a series to be the value of the analytic continuation of the Dirichlet series

at *s* = 0, if this exists and is unique. This method is sometimes confused with zeta function regularization.

If *s* = 0 is an isolated singularity, the sum is defined by the constant term of the Laurent series expansion.

Main article: Zeta function regularization |

If the series

(for positive values of the *a*_{n}) converges for large real *s* and can be analytically continued along the real line to *s* = −1, then its value at *s* = −1 is called the zeta regularized sum of the series *a*_{1} + *a*_{2} + ... Zeta function regularization is nonlinear. In applications, the numbers *a*_{i} are sometimes the eigenvalues of a self-adjoint operator *A* with compact resolvent, and *f*(*s*) is then the trace of *A*^{−s}. For example, if *A* has eigenvalues 1, 2, 3, ... then *f*(*s*) is the Riemann zeta function, *ζ*(*s*), whose value at *s* = −1 is −1/12, assigning a value to the divergent series 1 + 2 + 3 + 4 + .... Other values of *s* can also be used to assign values for the divergent sums *ζ*(0) = 1 + 1 + 1 + ... = −1/2, *ζ*(−2) = 1 + 4 + 9 + ... = 0 and in general

where *B _{k}* is a Bernoulli number.

If *J*(*x*) = Σ*p*_{n}*x*^{n} is an integral function, then the *J* sum of the series *a*_{0} + ... is defined to be

if this limit exists.

There is a variation of this method where the series for *J* has a finite radius of convergence *r* and diverges at *x* = *r*. In this case one defines the sum as above, except taking the limit as *x* tends to *r* rather than infinity.

In the special case when *J*(*x*) = *e*^{x} this gives one (weak) form of Borel summation.

Valiron's method is a generalization of Borel summation to certain more general integral functions *J*. Valiron showed that under certain conditions it is equivalent to defining the sum of a series as

where *H* is the second derivative of *G* and *c*(*n*) = *e*^{−G(n)}, and *a*_{0} + ... + *a*_{h} is to be interpreted as 0 when *h* < 0.

Suppose that *dμ* is a measure on the real line such that all the moments

are finite. If *a*_{0} + *a*_{1} + ... is a series such that

converges for all *x* in the support of *μ*, then the (*dμ*) sum of the series is defined to be the value of the integral

if it is defined. (If the numbers *μ*_{n} increase too rapidly then they do not uniquely determine the measure *μ*.)

For example, if *dμ* = *e*^{−x} *dx* for positive *x* and 0 for negative *x* then *μ*_{n} = *n*!, and this gives one version of Borel summation, where the value of a sum is given by

There is a generalization of this depending on a variable *α*, called the (B′,*α*) sum, where the sum of a series *a*_{0} + ... is defined to be

if this integral exists. A further generalization is to replace the sum under the integral by its analytic continuation from small *t*.

This summation method works by using an extension to the real numbers known as the hyperreal numbers. Since the hyperreal numbers include distinct infinite values, these numbers can be used to represent the values of divergent series. The key method is to designate a particular infinite value that is being summed, usually , which is used as a unit of infinity. Instead of summing to an arbitrary infinity (as is typically done with ), the BGN method sums to the specific hyperreal infinite value labeled . Therefore, the summations are of the form

This allows the usage of standard formulas for finite series such as arithmetic progressions in an infinite context. For instance, using this method, the sum of the progression is , or, using just the most significant infinite hyperreal part, .^{[6]}

Hardy (1949, chapter 11).

Main article: Hölder summation |

In 1812 Hutton introduced a method of summing divergent series by starting with the sequence of partial sums, and repeatedly applying the operation of replacing a sequence *s*_{0}, *s*_{1}, ... by the sequence of averages *s*_{0} + *s*_{1}/2, *s*_{1} + *s*_{2}/2, ..., and then taking the limit (Hardy 1949, p. 21).

The series *a*_{1} + ... is called Ingham summable to *s* if

Albert Ingham showed that if *δ* is any positive number then (C,−*δ*) (Cesàro) summability implies Ingham summability, and Ingham summability implies (C,*δ*) summability Hardy (1949, Appendix II).

The series *a*_{1} + ... is called Lambert summable to *s* if

If a series is (C,*k*) (Cesàro) summable for any *k* then it is Lambert summable to the same value, and if a series is Lambert summable then it is Abel summable to the same value Hardy (1949, Appendix II).

The series *a*_{0} + ... is called Le Roy^{[citation needed]} summable to *s* if

Hardy (1949, 4.11)

The series *a*_{0} + ... is called Mittag-Leffler (M) summable to *s* if

Hardy (1949, 4.11)

Main article: Ramanujan summation |

Ramanujan summation is a method of assigning a value to divergent series used by Ramanujan and based on the Euler–Maclaurin summation formula. The Ramanujan sum of a series *f*(0) + *f*(1) + ... depends not only on the values of *f* at integers, but also on values of the function *f* at non-integral points, so it is not really a summation method in the sense of this article.

The series *a*_{1} + ... is called (R,*k*) (or Riemann) summable to *s* if

Hardy (1949, 4.17)
The series *a*_{1} + ... is called R_{2} summable to *s* if

Main article: Riesz mean |

If *λ*_{n} form an increasing sequence of real numbers and

then the Riesz (R,*λ*,*κ*) sum of the series *a*_{0} + ... is defined to be

The series *a*_{1} + ... is called VP (or Vallée-Poussin) summable to *s* if

where is the gamma function. Hardy (1949, 4.17).