← LibraryMulti-index Notation for Differential Operators | KEVOS® MathematicsProject Delivery · Project ManagementLesson 44/72← PrevNext →
ArticlePublished 9 Aug 202618 min readBy Kevin Jogin
Skip to content
KEVOS® Engineering · Mathematics Knowledge Library

EngineeringMathematicsFoundation

Multi-index Notation for Differential Operators

A multi-index is a vector αn of exponents. Writing xα and β turns unwieldy products into single symbols and makes the canonical basis of An indexable by pairs (α,β)2n.

Collection Algebraic D-modulesTopic stream weyl-algebraSource Ch. 1 §2Reading time 20 minPage ID KVS-ENG-MATH-0332

Overview

Working in the Weyl algebra An means writing products such as x13x20x35122430 over and over. Written out, these expressions are unreadable and the patterns in them are invisible. Multi-index notation replaces each block of exponents by a single vector: α=(3,0,5), β=(2,4,0), and the operator becomes xαβ.

This is more than an abbreviation. Once exponents are vectors, statements about differential operators become statements about n, and the combinatorics can be done in the exponent lattice. The canonical basis of An is indexed by n×n=2n; the degree of an operator is a maximum of lengths of multi-indices; counting basis elements becomes counting lattice points in a simplex; and the general Leibniz rule becomes a sum over the componentwise interval below a multi-index.

Two facts about the notation deserve emphasis at the outset, because both are common sources of confusion. First, multi-indices carry two different orderings - the partial order by components and the total preorder by length - and almost every statement in the subject depends on which is meant. Second, xαβ denotes a composite operator, not a function: it means multiply by xα after differentiating, whereas β(xα) means differentiate the polynomial xα. Those two expressions are different objects and are almost never equal.

This page collects the definitions, the identities that get used constantly, and the counting formulas that later dimension arguments rely on.

Definition

Multi-index and its invariants

A multi-index in n variables is an element α=(α1,,αn) of n, where ={0,1,2,}. Attached to it are

  • its length |α|=α1+α2++αn;
  • its factorial α!=α1!α2!αn!;
  • the monomial xα=x1α1x2α2xnαnK[X];
  • the iterated derivative α=1α12α2nαn.

Because the xi commute with each other and the i commute with each other, neither xα nor α depends on the order in which the factors are written.

Order, sum and binomial coefficient

For multi-indices α,γn write γα when γiαi for every i. Sums and differences are componentwise, with αγ defined only when γα. The multi-index binomial coefficient is

(αγ)=i=1n(αiγi)=α!γ!(αγ)!(γα).

Finally, ei denotes the multi-index with 1 in position i and 0 elsewhere, so that xei=xi and ei=i.

Note

A pair (α,β) of multi-indices in n is itself a multi-index in 2n, of length |α|+|β|. That observation is what lets the degree of an element of An be defined by a single length, and it is why the counting formulas below are stated in 2n variables when they are applied to An.

Core Concepts

The exponent lattice

Think of a monomial not as a product but as the point of n recording its exponents. Multiplying monomials adds points; the total degree is the 1 length; divisibility of one monomial by another is the componentwise order. Under this dictionary, questions about polynomials become questions about lattice points, and a surprising number of them become elementary counting problems.

Two orders, doing different jobs

Ordering by length is total but coarse: (2,0) and (0,2) and (1,1) all have length 2. Ordering componentwise is fine but partial: (2,0) and (0,2) are incomparable. Each order governs different statements. Derivatives are controlled by the componentwise order - β(xα) is non-zero exactly when βα - while filtrations and degrees are controlled by length. Confusing them produces false statements that look plausible.

Operators versus values

In An the symbol xα means the operator 'multiply by xα', so xαβ is a composition of operators. Applying that composition to a polynomial f gives xαβ(f). The parenthesis is doing real work: βxα without parentheses is the operator obtained by composing in the other order, and by the commutation relations it equals xαβ plus strictly lower-degree corrections.

Why this notation and not tensor indices

Multi-index notation is adapted to commuting symbols. It compresses x1α1xnαn into one symbol precisely because reordering the factors changes nothing. That is why it is used for the variables and for the derivatives separately, but never across the two blocks: x's and 's do not commute, so no single multi-index in 2n can stand for an arbitrary word in the generators. The pair (α,β) stands for the canonical arrangement only.

Key Equations

The identity used most often is the value of an iterated derivative on a monomial. For α,βn,

β(xα)={α!(αβ)!xαβifβα,0otherwise.
(2.1)

The special case β=α gives α(xα)=α!, the fact that drives the proof that the canonical basis is a basis. The general Leibniz rule reads

α(fg)=γα(αγ)γ(f)αγ(g),f,gK[X].
(2.2)

Its operator form, obtained by taking f=xα and reading both sides as elements of An, is the rule for straightening a product:

βxα=γmin(α,β)(βγ)α!(αγ)!xαγβγ.
(2.3)

The multinomial theorem and the Taylor expansion of a polynomial both become one-line statements:

(x1++xn)m=|α|=mm!α!xα,f(x)=ααf(a)α!(xa)α.
(2.4)

Finally the two counting formulas. The number of multi-indices of exact length m and of length at most m in n variables are

#{αn:|α|=m}=(m+n1n1),#{αn:|α|m}=(m+nn).
(2.5)

Applied with 2n in place of n, the second formula counts the canonical basis elements of An of degree at most m, giving (m+2n2n) - the starting point of every dimension computation later in the theory.

Variable Definitions

n
the number of variables
α,β,γ,σ
multi-indices, elements of n
αi
the i-th entry of the multi-index α, a non-negative integer
|α|
the length α1++αn
α!
the factorial α1!αn!
γα
the componentwise order: γiαi for every i
(αγ)
the product of the entrywise binomial coefficients
ei
the multi-index with 1 in position i, zeros elsewhere
xα
the monomial x1α1xnαn, or the operator multiplying by it
β
the iterated partial derivative 1β1nβn

Properties and Behaviour

Basic identities

  • xαxβ=xα+β and αβ=α+β; the exponents add.
  • |α+β|=|α|+|β|, so length is additive.
  • degx(xα)=|α|, and the monomials {xα:|α|m} are a basis of the polynomials of degree at most m.
  • (α+β)!α!β! in general; the correct relation is (α+βα)=(α+β)!/(α!β!), an integer.
  • β(xα)=0 whenever |β|>|α|, but the converse fails: (1,0)(x25)=0 with |β|=1<5.

Lengths bound components

If γα then |γ||α|, but not conversely. Consequently every statement proved using the componentwise order implies the corresponding statement about lengths, and no statement about lengths implies anything componentwise. The one useful partial converse: if γα and |γ|=|α| then γ=α. This is the pigeonhole step at the heart of the canonical basis proof.

Degree in the Weyl algebra

For D=cαβxαβ0 in canonical form, degD=max{|α|+|β|:cαβ0}, with deg0=. Straightening formula (2.3) shows the correction terms all have |αγ|+|βγ|=|α|+|β|2|γ|, strictly smaller for γ0: reordering an operator lowers degree by an even amount, never raises it. This is the multi-index reason behind the commutator bound deg[D,D]degD+degD2.

Examples and Special Cases

Reading a multi-index

With n=3 and α=(2,0,1): |α|=3, α!=2!0!1!=2, xα=x12x3, and α=123. Note x2 does not appear, because its exponent is zero; a multi-index always has n entries even when most are zero.

Formula (2.1) in action

Take n=2, α=(3,2), β=(1,2). Since βα, formula (2.1) gives α!/(αβ)!=(3!2!)/(2!0!)=12/2=6 and 122(x13x22)=6x12. Checking directly: 22(x13x22)=2x13, then 1(2x13)=6x12.

Incomparable indices

α=(2,0) and β=(0,2) have the same length but neither is componentwise below the other. Accordingly β(xα)=22(x12)=0 and α(xβ)=0 as well, while α(xα)=2.

Counting

In two variables there are (2+11)=3 multi-indices of length exactly 2, namely (2,0),(1,1),(0,2), and (2+22)=6 of length at most 2, adding (0,0),(1,0),(0,1). Both match formula (2.5).

Worked Example

Straightening 12x12x2 in A2, and counting a filtration piece

  1. Step 1 - name the indices

    Take n=2, α=(2,1) and β=(1,1), so that xα=x12x2 and β=12. We want the canonical form of the operator βxα. Here min(α,β)=(1,1), so γ in (2.3) ranges over (0,0),(1,0),(0,1),(1,1).

  2. Step 2 - evaluate the four terms

    Every binomial coefficient (βγ) equals 1 here, since each entry of β is 1. With α!=2!1!=2, the factor α!/(αγ)! is:

    • γ=(0,0): (αγ)!=2, factor 1, term x12x212;
    • γ=(1,0): (αγ)!=1!1!=1, factor 2, term 2x1x22;
    • γ=(0,1): (αγ)!=2!0!=2, factor 1, term x121;
    • γ=(1,1): (αγ)!=1!0!=1, factor 2, term 2x1.
    12x12x2=x12x212+2x1x22+x121+2x1.
  3. Step 3 - verify on test polynomials

    Apply both sides to 1: the left side is 12(x12x2)=1(x12)=2x1; on the right only the last term survives, giving 2x1.

    Apply both sides to x2: the left side is 12(x12x22)=1(2x12x2)=4x1x2; the right side gives 0+2x1x2+0+2x1x2=4x1x2.

    Apply both sides to x1x2: the left side is 12(x13x22)=1(2x13x2)=6x12x2; the right side gives x12x2+2x12x2+x12x2+2x12x2=6x12x2. Three independent checks agree.

  4. Step 4 - count the operators of degree at most 3

    The operator just computed has degree |α|+|β|=3. By (2.5) with 2n=4, the space of elements of A2 of degree at most 3 has dimension (3+44)=(74)=35.

    Cross-check by summing over exact degrees, using the first formula in (2.5) with n replaced by 4: degrees 0,1,2,3 contribute (33)=1, (43)=4, (53)=10 and (63)=20. The total is 1+4+10+20=35.

Result

12x12x2=x12x212+2x1x22+x121+2x1, an element of degree 3 in a space of dimension 35. Every correction term has degree 32|γ|, dropping by two at each step, exactly as the degree bound on commutators predicts.

Applications and Industry Use

In a mathematics topic, this section covers downstream use inside mathematics, computing and engineering rather than a manufactured product.

  • Indexing the canonical basis. The basis of An is {xαβ} indexed by 2n, so an operator is a finitely supported function 2nK. Every algorithm in the subject manipulates that function.
  • Filtrations. The Bernstein filtration collects the xαβ with |α|+|β|m; the order filtration collects those with |β|m. Both are conditions on multi-index lengths, which is why their Hilbert functions are lattice point counts.
  • Partial differential equations. A linear PDE with polynomial coefficients is written |β|kaβ(x)βu=0; the notation is standard in analysis for the same reasons it is standard here.
  • Symbols and characteristic varieties. Replacing β by ξβ in the top-degree part turns an operator into a polynomial on the cotangent space; the substitution is literally a relabelling of multi-indices.
  • Computer algebra. Sparse polynomial and operator arithmetic is implemented as arithmetic on exponent vectors with hashed lookup, and term orders are orders on 2n.

Standards and Codes

For mathematics, the relevant standards are notation, numeric and markup standards together with reference implementations.

  • ISO 80000-2 fixes for the non-negative integers - note that the standard includes 0, which is what makes n the right index set here - and for partial differentiation.
  • ISO/IEC 40314 (MathML 3.0) is the encoding used for the expressions on this page, so that a superscripted multi-index remains machine-readable rather than becoming an image.
  • There is no standard for multi-index notation itself; the conventions used here (length written |α|, factorial written α!, componentwise order written ) are the ones common to the analysis and D-module literature and are the ones Coutinho uses.
  • One genuine notational divergence to watch: some texts write Dα or xα for α, and some analysts insert a factor of i or 2πi into Dα to make Fourier transforms cleaner. Those conventions are not used in this collection.

Computational Notes

Read this as the manufacturing section of the template: how the object is actually built by machine, at what cost, and where the computation stops being decidable.

An operator in An is stored as a dictionary keyed by the pair (α,β)2n, with K-coefficients as values. Addition merges dictionaries; multiplication applies (2.3) term by term. Three practical points follow from the notation itself.

  1. Multiplying xαβ by xαβ produces one term for each γmin(α,β), so the number of output terms is i(min(αi,βi)+1). Sparse operators of low degree multiply cheaply; dense high-degree ones do not.
  2. The leading exponent pair of a product is the sum of the leading exponent pairs, because all correction terms have strictly smaller length. This is what makes term orders on 2n compatible with multiplication and non-commutative Groebner bases possible.
  3. Factorial ratios α!/(αγ)! should be computed as falling factorials iαi(αi1)(αiγi+1) rather than as quotients of factorials, to avoid large intermediate integers.

Macaulay2's Dmodules package, Singular's Plural/dmod.lib and SageMath's ore_algebra all use exponent-vector representations of this kind, and their documentation states the same convention: variables left, derivatives right.

Failure Modes and Common Mistakes

Confusing α! with |α|!

For α=(2,1), α!=2 but |α|!=3!=6. Both appear in the subject - α! in (2.1) and (2.2), |α|! in the multinomial coefficient m!/α! of (2.4) - and substituting one for the other silently corrupts every subsequent constant.

Assuming |β||α| implies βα

It does not. With β=(0,1) and α=(2,0) we have |β|=1<2=|α|, yet β(xα)=2(x12)=0. Formula (2.1) is governed by the componentwise order; only the vanishing statement |β|>|α|β(xα)=0 is governed by length.

Reading xαβ as a function value

xαβ is an operator; β(xα) is a polynomial. In A1 with α=β=1, the operator x sends f to xf, whereas (x)=1. Writing the two the same way, or cancelling the xα against the β, is the single most common slip when first working in An.

Dropping zero entries

A multi-index in n variables always has n entries. Writing α=(2,1) when n=3 is undefined; the intended index is (2,1,0). This matters when multi-indices are compared or added across expressions in different numbers of variables - for instance when Am is viewed inside An - where the shorter index must be padded with zeros.

Subtracting outside the allowed range

αγ is a multi-index only when γα; otherwise it has a negative entry and xαγ is not a polynomial. Sums such as (2.2) and (2.3) restrict the summation index for exactly this reason, and an implementation that omits the range check will silently produce Laurent monomials.

Historical Notes

Compressed exponent notation appears wherever many variables are differentiated many times, and it became standard through partial differential equations rather than algebra: Schwartz's theory of distributions in the late 1940s and Hormander's treatments of linear PDE fixed the now-universal conventions xα, α, |α| and α!. Algebraists inherited the notation together with the subject matter.

For the Weyl algebra the notation does something extra. Because the canonical basis is indexed by 2n, questions about An can be transported to the commutative polynomial ring in 2n variables, whose combinatorics is classical. That transport - degrees to lengths, filtration pieces to simplices of lattice points, operators to their symbols - is the technical device behind the dimension theory of Chapters 7 to 9 of Coutinho's book, and it is available only because the notation makes the index set visible.

Comparison

The two orderings of multi-indices and what each controls.
Componentwise order γαOrder by length |γ||α|
Type of orderPartialTotal on lengths, but many indices share a length
ControlsWhen γ(xα)0; divisibility of monomialsDegree, filtrations, Hilbert functions
Typical statementγ(xα)0γαdeg(DD)=degD+degD
Incomparable pairs(2,0) and (0,2)None; lengths are integers
Used inLeibniz rule, straightening formula (2.3)Bernstein filtration, degree, dimension theory
Relationγα implies |γ||α|The converse is false

Key Takeaways

Key points

  • A multi-index is a vector αn; it carries a length |α|, a factorial α!, a monomial xα and a derivative α.
  • The componentwise order decides when β(xα) is non-zero; the length decides degrees and filtrations. They are different orders and are not interchangeable.
  • β(xα)=α!/(αβ)!xαβ when βα, and 0 otherwise; in particular α(xα)=α!.
  • Straightening in An is the single formula (2.3), whose correction terms drop degree by 2|γ|.
  • A pair (α,β) is a multi-index in 2n, so counting elements of An of degree at most m gives (m+2n2n).
  • xαβ is an operator and β(xα) is a polynomial; keeping them apart avoids most beginners' errors.

FAQs

Does include 0 here?

Yes, and it must. The multi-index (0,,0) corresponds to the monomial 1 and to the identity operator, both of which are needed. This also matches ISO 80000-2, which defines to contain 0.

Why is there no single multi-index for a general element of An?

Because the generators do not all commute. A multi-index compresses a product of commuting factors; the x's commute among themselves and the 's do too, so each block gets its own index, but a word mixing them is not determined by its exponent counts. The pair (α,β) specifies a canonically ordered product only.

Is β(xα)=0 the same as |β|>|α|?

No. |β|>|α| is sufficient but not necessary. The exact criterion is that β fails to be componentwise below α; for instance 2(x15)=0 although |β|=1 and |α|=5.

How do I remember formula (2.3)?

Read it as a Leibniz rule. Moving β past xα, each derivative either passes through untouched or lands on the monomial. The index γ records which derivatives land: (βγ) counts the ways of choosing them, and α!/(αγ)! is what they produce.

Why do the counting formulas use 2n?

Because a canonical basis element of An is indexed by a pair of multi-indices in n variables, that is by one multi-index in 2n variables. The degree is its length, so counting operators of degree at most m is counting lattice points of length at most m in 2n, giving (m+2n2n).

Does any of this depend on the characteristic of K?

The notation does not, and formulas (2.1) to (2.3) hold over any commutative ring, with the integers α! and (αγ) read as multiples of 1. What changes in characteristic p is that these integers can become zero in K - for example (p)(x(p))=p!=0 - which is why the canonical basis result needs characteristic zero even though the notation does not.

What is the relationship to the notation used in PDE textbooks?

It is the same notation. A linear operator with polynomial coefficients written |β|kaβ(x)β in an analysis text is exactly an element of An in canonical form, grouped by β. The difference is that analysis usually allows smooth or holomorphic coefficients, whereas An insists on polynomial ones.

References

  1. S. C. Coutinho, A Primer of Algebraic D-modules, London Mathematical Society Student Texts 33, Cambridge University Press, 1995 - Ch. 1 §2, where multi-index notation is introduced and used for the canonical basis (1.2.1) and Lemma (1.2.2).
  2. L. Hormander, The Analysis of Linear Partial Differential Operators I, 2nd edition, Springer, 1990 - §1.1, for the standard analytic conventions.
  3. L. Schwartz, Theorie des distributions, Hermann, Paris, 1950-1951 - the source of much of the modern notation for iterated partial derivatives.
  4. J. C. McConnell and J. C. Robson, Noncommutative Noetherian Rings, Graduate Studies in Mathematics 30, revised edition, American Mathematical Society, 2001 - Ch. 1, for the Weyl algebra in this notation.
  5. M. Saito, B. Sturmfels and N. Takayama, Groebner Deformations of Hypergeometric Differential Equations, Algorithms and Computation in Mathematics 6, Springer, 2000 - Ch. 1, for exponent-vector representations and term orders on 2n.
  6. R. P. Stanley, Enumerative Combinatorics, Volume 1, 2nd edition, Cambridge University Press, 2011 - Ch. 1, for the lattice point counts (2.5).
  7. ISO 80000-2:2019, Quantities and units - Part 2: Mathematics, International Organization for Standardization.
  8. ISO/IEC 40314:2016, Information technology - Mathematical Markup Language (MathML) Version 3.0, International Organization for Standardization.

AI Suggested Questions

  • Verify formula (2.3) for n=1, α=3, β=2, and check the result by applying both sides to 1, x and x2.
  • Count the elements of the canonical basis of A3 of degree exactly 4, and check the answer against the total for degree at most 4.
  • Prove the general Leibniz rule (2.2) by induction on |α|.
  • Show that (α+βα) is an integer for all multi-indices, and interpret it as a count.
  • Explain why the number of terms produced by one application of (2.3) is i(min(αi,βi)+1).
  • Give an example of multi-indices with |γ||α|, γα, and γ(xα)0, and one with γ(xα)=0.
  • Describe a term order on 2n that refines total degree, and check that it is compatible with multiplication in An.

Continue learning

A Roadmap Through Algebraic D-module Theory | KEVOS® MathematicsArticle · Project ManagementAutomorphisms of the Weyl Algebra | KEVOS® MathematicsArticle · Project ManagementBernstein's Inequality | KEVOS® MathematicsArticle · Project ManagementCanonical Form of an Element of the Weyl Algebra | KEVOS® MathematicsArticle · Project Management