Math 2.0

Numerical analysis Theorem

Dahlquist barriers

Builds on Linear multistep methodBackward differentiation formula

Statement

Consider a linear kk-step method ∑i=0kαiyj+i=h∑i=0kβifj+i\sum_{i=0}^k \alpha_i y_{j+i} = h \sum_{i=0}^k \beta_i f_{j+i} with αk=1\alpha_k = 1 and characteristic polynomials ρ(ζ)=∑iαiζi\rho(\zeta) = \sum_i \alpha_i \zeta^i, σ(ζ)=∑iβiζi\sigma(\zeta) = \sum_i \beta_i \zeta^i, where ρ\rho and σ\sigma have no common root. The method has order pp if ρ(ez)−zσ(ez)=O(zp+1)\rho(e^z) - z \sigma(e^z) = O(z^{p+1}) as z→0z \to 0. It is zero-stable if every root of ρ\rho has modulus at most 11 and the roots of modulus 11 are simple. It is A-stable if, for every complex λ\lambda with Re⁡λ<0\operatorname{Re} \lambda < 0 and every h>0h > 0, every sequence produced by the method for y′=λyy' = \lambda y tends to 00.

  1. First Dahlquist barrier: a zero-stable linear kk-step method has order p≤k+2p \le k + 2 if kk is even, and p≤k+1p \le k + 1 if kk is odd. If it is moreover explicit, then p≤kp \le k.
  2. Second Dahlquist barrier: an A-stable linear multistep method has order p≤2p \le 2, and no explicit linear multistep method is A-stable.

Remarks

For the first barrier, the substitution ζ=1+w1−w\zeta = \frac{1 + w}{1 - w} maps the left half plane Re⁡w<0\operatorname{Re} w < 0 onto the open unit disk. It turns the order conditions into conditions on the power series of a rational function of ww. Zero-stability forces sign conditions on the coefficients that are incompatible with a higher order. For the second barrier, A-stability means that ρ(ζ)−zσ(ζ)\rho(\zeta) - z \sigma(\zeta) has all its roots in the open unit disk whenever Re⁡z<0\operatorname{Re} z < 0. Hence Re⁡ρ(ζ)σ(ζ)≥0\operatorname{Re} \frac{\rho(\zeta)}{\sigma(\zeta)} \ge 0 for ∣ζ∣>1\abs{\zeta} > 1, and this positivity is incompatible with order greater than 22. Example: Simpson's two-step method yj+2−yj=h3(fj+2+4fj+1+fj)y_{j+2} - y_j = \frac{h}{3} (f_{j+2} + 4 f_{j+1} + f_j) has order 4=k+24 = k + 2 and is zero-stable, so the first barrier is sharp for k=2k = 2. The trapezoidal rule yj+1−yj=h2(fj+1+fj)y_{j+1} - y_j = \frac{h}{2} (f_{j+1} + f_j) has order 22 and is A-stable, so the second barrier is sharp; among A-stable methods of order 22 it has the smallest error constant. The backward differentiation formulas of order 33 to 66 are zero-stable but not A-stable, as the second barrier requires.

Going deeper

For a student

A kk-step method can reach order 2k2k, but such methods are not zero-stable for k≥3k \ge 3. The first barrier explains why practical codes use Adams methods instead, of order kk or k+1k + 1 for kk steps. The second barrier explains why stiff solvers use BDF methods. Above order 22 these methods give up full A-stability, but they keep stability in a large sector of the left half plane. It also motivates implicit Runge–Kutta methods, which escape the barrier because they are not linear multistep methods: the Gauss–Legendre methods are A-stable of every even order. Both barriers were proved by Germund Dahlquist.

Reading path

To understand this page from the ground up, read these 3 pages in order. The foundations come first.

Show the path as a graph
Each page uses the pages above it. The colour of the bar on a box shows how strong the evidence is, and the box also gives it as text. The list has the same pages.
  • no evidence level
  1. Backward differentiation formula concept
  2. Linear multistep method concept
  3. Dahlquist barriers statement

Further reading

  1. Web https://en.wikipedia.org/wiki/Linear_multistep_method
About this page
TagBQBZ7QMJE2
URIm2://BQBZ7QMJE2
Kindstatement, theorem
WritingW0 stub
MSC65L06, 65L20
StatusThis copy of the site has no build data. Run m2 status BQBZ7QMJE2.
AboutLinear multistep method
UsesLinear multistep method, Backward differentiation formula
Changed2 November 2026

Markdown · JSON · Edit on GitHub

What this page needs for W1

  • an inventory row: the row of the inventory that the node accounts for (0 of 1)
  • a source with a locator (0 of 1)

The writing levels say how complete the text of a page is. They do not measure its truth.

What links here

No page links to this page yet.