KEVOS
ArticlesServicesCase studiesAboutContact
ArticlesServicesCase studiesAboutContact
← ArticlesMulti-index Notation for Differential Operators | KEVOS® MathematicsEngineering · Engineering MathematicsLesson 177/883← PrevNext →
ArticlePublished 15 Aug 2026Updated 9 Aug 202618 min readBy KEVOS® Knowledge Library
On this page

Ask about this page

KEVOS AIMulti-index Notation for Differential Operators | KEVOS® Mathematics

KEVOS knowledge first · trusted web sources when needed

Skip to content
KEVOS® Engineering · Mathematics Knowledge Library
  1. Home
  2. Engineering
  3. Mathematics
  4. Multi-index Notation for Differential Operators

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 x13x20x35∂12∂24∂30 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α2⋯xnαn∈K[X];
  • the iterated derivative ∂α=∂1α1∂2α2⋯∂nα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α1⋯xnα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,g∈K[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)α!(x−a)α.
(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+n−1n−1),#{α∈ℕ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α1⋯xnαn, or the operator multiplying by it
∂β
the iterated partial derivative ∂1β1⋯∂nβ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+degD′−2.

Examples and Special Cases

Reading a multi-index

With n=3 and α=(2,0,1): |α|=3, α!=2!⋅0!⋅1!=2, xα=x12x3, and ∂α=∂12∂3. 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 ∂1∂22(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 ∂1∂2x12x2 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 ∂β=∂1∂2. 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 x12x2∂1∂2;
    • γ=(1,0): (α−γ)!=1!⋅1!=1, factor 2, term 2x1x2∂2;
    • γ=(0,1): (α−γ)!=2!⋅0!=2, factor 1, term x12∂1;
    • γ=(1,1): (α−γ)!=1!⋅0!=1, factor 2, term 2x1.
    ∂1∂2x12x2=x12x2∂1∂2+2x1x2∂2+x12∂1+2x1.
  3. Step 3 - verify on test polynomials

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

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

    Apply both sides to x1x2: the left side is ∂1∂2(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

∂1∂2x12x2=x12x2∂1∂2+2x1x2∂2+x12∂1+2x1, an element of degree 3 in a space of dimension 35. Every correction term has degree 3−2|γ|, 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 ℕ2n→K. 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(αi−1)⋯(α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.

Related Engineering Topics

  • The Weyl Algebra Defined as a Ring of Operators

    The algebra whose elements this notation describes.

    Prerequisite
  • Canonical Form of an Element of the Weyl Algebra

    Where the pair of multi-indices becomes an index for a basis.

    Next step
  • The Degree of a Differential Operator

    Degree as the maximum length of an occurring multi-index.

    Next step
  • What Algebraic D-modules Are

    The wider setting in which this bookkeeping pays off.

    Background
  • The Weyl Algebra and the Birth of Quantum Mechanics

    Why operators rather than polynomials are the objects of study.

    Background

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.

On this page

  1. Overview
  2. Definition
  3. Core Concepts
  4. Key Equations
  5. Variable Definitions
  6. Properties and Behaviour
  7. Examples and Special Cases
  8. Worked Example
  9. Applications and Industry Use
  10. Standards and Codes
  11. Computational Notes
  12. Failure Modes and Common Mistakes
  13. Historical Notes
  14. Comparison
  15. Key Takeaways
  16. FAQs
  17. Related Engineering Topics
  18. References
  19. AI Suggested Questions

Part of the KEVOS® Engineering › Mathematics knowledge library, collection Algebraic D-modules.

This page is original explanatory prose. Source results are attributed in the References section; the underlying textbook is in copyright and is cited, not reproduced.

Page ID KVS-ENG-MATH-0332 Taxonomy /engineering/mathematics Category ID ENG / ENG-MATH Level Foundation Reading time 20 min Page version 1.0.0 Content version 2026.08 Reviewed 2026-08-09

Continue learning

A Roadmap Through Algebraic D-module Theory | KEVOS® MathematicsArticle · Engineering MathematicsAutomorphisms of the Weyl Algebra | KEVOS® MathematicsArticle · Engineering MathematicsBernstein's Inequality | KEVOS® MathematicsArticle · Engineering MathematicsCanonical Form of an Element of the Weyl Algebra | KEVOS® MathematicsArticle · Engineering Mathematics
KEVOS · Engineering, manufacturing and project improvement
ArticlesServicesCase studiesAboutContact
© 2026 KEVOS®