05.09.10 · symplectic / kam

KP hierarchy, Sato Grassmannian, and tau-functions

shipped3 tiersLean: none

Anchor (Master): Sato 1981 RIMS Kokyuroku 439 (foundational); Date-Jimbo-Kashiwara-Miwa 1981-1983 RIMS Kokyuroku series (transformation groups, vertex operators); Segal-Wilson 1985 Publ. IHES 61 (loop-group reformulation); Mulase 1984 Adv. Math. 54 (algebro-geometric tau-functions); Shiota 1986 Invent. Math. 83 (Schottky-problem characterisation); Kontsevich 1992 Comm. Math. Phys. 147 (intersection theory on moduli spaces of curves)

Intuition Beginner

A soliton keeps its shape because dispersion and nonlinearity exactly cancel. The Korteweg-de Vries equation is the simplest such model in one space dimension, and its periodic siblings (the finite-gap potentials) are completely described by algebraic curves and their Jacobians. The story does not end there: the KdV equation is the first member of an infinite tower of commuting equations, the KdV hierarchy, and the entire tower extends naturally to a two-dimensional analogue called the Kadomtsev-Petviashvili (KP) equation .

The KP hierarchy is built around a single universal object discovered by Mikio Sato in 1981: the Sato Grassmannian, an infinite-dimensional space whose points are subspaces of formal Laurent series. Each point of this Grassmannian determines a solution of the KP hierarchy, and the time evolution of KP is nothing other than the natural group action of a one-parameter family of exponentials on the Grassmannian.

To each point of the Sato Grassmannian Sato attached a single function, the tau-function — the Plücker coordinate of the moving point measured against a fixed reference subspace. The KP solution recovers from by taking two logarithmic derivatives in the space variable. Every soliton solution, every finite-gap solution, the Witten-Kontsevich generating function of intersection numbers on the moduli space of curves — all are tau-functions of one Sato Grassmannian point or another. The KP picture is the universal language in which all of soliton theory is written.

Visual Beginner

Picture three layered objects. At the bottom, a 2-dimensional travelling wave — the KP profile, looking like a KdV soliton extended uniformly in a second direction with a slow modulation. In the middle, the Sato Grassmannian drawn as a vast "sea" of subspaces of formal Laurent series; a single point inside this sea evolves under a one-parameter family of exponential maps, drawing out a smooth trajectory. At the top, the tau-function — a single -valued function of infinitely many variables, the Plücker coordinate of the moving point measured against a fixed reference subspace.

The picture captures Sato's insight: the infinite tower of nonlinear KP equations is, in the right coordinates, nothing but the linear translation of one point in an infinite-dimensional Grassmannian.

Worked example Beginner

The two-soliton solution of KP is the easiest nondegenerate tau-function to write down.

Step 1. Fix two complex numbers , and two phases , . Define , so and .

Step 2. Form the tau-function . This is the simplest rational-exponential tau-function: three terms summing to a positive function.

Step 3. Compute the KP solution (two derivatives in ). One differentiation gives ; differentiating again yields a tanh-squared-shaped 2-soliton, with one wave of crest speed in the diagonal direction and a second wave of crest speed in a different diagonal direction.

Step 4. At , the snapshot is a function of alone with two peaks near and small wings elsewhere. As increases, the two soliton crests separate at their different speeds; for they tilt and exhibit the resonant Y-shaped interaction that is the visual signature of KP.

What this tells us: any finite sum of pure exponentials gives a multi-soliton tau-function. Geometrically, this corresponds to a finite-rank subspace of the Sato Grassmannian — a point spanned by finitely many exponential vectors — passing through a generic stratum.

Check your understanding Beginner

Formal definition Intermediate+

Let be the field of formal Laurent series in , with the standard subspace decomposition where and . Let denote the projection along .

Definition (Sato Grassmannian). The Sato Grassmannian is the set of subspaces such that the restriction is a Fredholm operator of index zero. Equivalently, is commensurable with : the subspaces and are such that has finite codimension in both and , with the same finite codimension on each side.

A point admits an admissible basis with for all but finitely many .

Definition (KP pseudo-differential operator). Let be the algebra of formal pseudo-differential operators in , i.e. expressions with finitely many positive powers and arbitrarily many negative powers, composed using the Leibniz rule . For any write with the differential part and the Volterra part.

Definition (KP hierarchy in Lax-Sato form). The KP hierarchy is the system of evolution equations on the pseudo-differential operator $$ L = \partial_x + u_2 \partial_x^{-1} + u_3 \partial_x^{-2} + u_4 \partial_x^{-3} + \cdots, $$ given by the infinite tower of flows $$ \frac{\partial L}{\partial t_n} = [(L^n)+, L], \qquad n \geq 1, $$ with . The flows commute pairwise: $[\partial{t_n}, \partial_{t_m}] L = 0$.

The KP equation. The first informative flows are

  • : , which is identically zero up to relabelling and gives .
  • : (heat-type), with .
  • : on its restriction to , giving KdV.

Combining the flows with produces the KP equation $$ \frac{3}{4} u_{yy} = \left(u_t - \frac{1}{4}(6 u u_x + u_{xxx})\right)_x, $$ a 2+1-dimensional integrable PDE.

Definition (Tau-function). Fix and consider the family of subspaces , , acting by multiplication on Laurent series. The tau-function of is the Plücker coordinate $$ \tau_W(t_1, t_2, \ldots) = \det(\pi_+|{W(t)} : W(t) \to H+), $$ defined up to overall scalar via the determinantal line bundle on . The function is a solution of the KP hierarchy, and every KP solution decaying at arises this way.

Hirota bilinear form. The KP equation, restated in terms of , is the bilinear KP equation $$ (D_1^4 + 3 D_2^2 - 4 D_1 D_3), \tau \cdot \tau = 0, $$ where the Hirota derivative acts on a pair of functions by $$ D_j^k f \cdot g = \left.\frac{\partial^k}{\partial s^k}\right|_{s=0} f(t_j + s) g(t_j - s). $$ Equivalently, the entire KP hierarchy is the system of bilinear identities (the Hirota equations) obtained from the residue identity for all , where and .

Counterexamples to common slips

  • Not every formal subspace lies in . The commensurability / Fredholm condition is genuine; a generic subspace of has either positive or negative virtual dimension relative to and lies in a different connected component of the universal Grassmannian, indexed by an integer (the Fredholm index).
  • The KP hierarchy is not a single equation. It is an infinite tower of compatible PDEs in infinitely many time variables; the KP equation proper is the compatibility of two specific flows ( and ).
  • Tau-functions are defined up to scalar. Multiplying by a non-vanishing function with — i.e. — does not change . The determinantal line bundle is the geometric reason for the ambiguity.
  • Hirota bilinear equations are not single-function ODEs. Each Hirota equation is a quadratic identity in , never a linear one; the trick is that the bilinear form makes the equation polynomial in the exponential building blocks , which is why generates exact -soliton solutions.

Key theorem with proof Intermediate+

Theorem (Sato 1981 — KP flows linearise on the Sato Grassmannian). Let be a point of the Sato Grassmannian and let be its tau-function. Then satisfies the KP hierarchy. Conversely, every formal-power-series solution of the KP hierarchy decaying at infinity in arises this way from a unique (up to scalar). The KP flows act on as linear translations .

Proof (Sato 1981). The proof has three structural steps.

Step 1. From the Grassmannian to the wave function. Fix with admissible basis . The wave function $$ \psi_W(t, z) = \frac{1}{\tau_W(t)} \exp\Big(\sum_{n \geq 1} t_n z^n\Big) \tau_W\big(t - [z^{-1}]\big) $$ admits the asymptotic expansion at and is, by Sato's construction, the unique generator of normalised so that .

Step 2. From the wave function to the Lax operator. Define the pseudo-differential operator where is the dressing operator, the unique formal with when acts on as . Then is a first-order pseudo-differential operator and . The Lax-Sato flows $$ \frac{\partial L}{\partial t_n} = [(L^n)+, L] $$ are equivalent to the family of evolution equations $\partial{t_n} \psi_W = (L^n)_+ \psi_W$ on the wave function.

Step 3. Plücker identity. The Plücker relations for the determinantal line bundle on the infinite Grassmannian are exactly the Hirota bilinear identities: $$ \mathrm{Res}_{z = \infty} \big(\psi_W(t, z) \cdot \psi^_W(t', z)\big) = 0 \quad \text{for all } t, t', $$ where $\psi^\psi^*W = \tau_W^{-1} e^{-\sum t_n z^n} \tau_W(t + [z^{-1}])\tau$ is $$ (D_1^4 + 3 D_2^2 - 4 D_1 D_3) \tau_W \cdot \tau_W = 0, $$ which by the change of variables $u = 2 \partial{t_1}^2 \log \tau_W$ is the KP equation.

Conversely, given a KP solution , one recovers the unique tau-function from together with the higher-flow constraints, and the point as the column span of the matrix of formal asymptotic coefficients of the wave function. The KP flows are translations on the Grassmannian by Step 2 and the linearity of the action .

Bridge. The Sato theorem identifies an infinite-dimensional nonlinear PDE system (the KP hierarchy in countably many time variables) with a linear translation on a single infinite-dimensional Grassmannian. The dictionary linearises everything: -soliton solutions are rank- subspaces; finite-gap solutions of KP are the subspaces stabilised by a finitely-generated commutative subalgebra of multiplication operators on (Krichever-Mulase 1985), i.e. the affine coordinate rings of compact algebraic curves 05.09.09; the boson-fermion correspondence (Date-Jimbo-Kashiwara-Miwa 1981-1983) reinterprets as the projectivisation of a semi-infinite wedge space, on which acts. The picture extends Liouville-Arnold integrability 05.02.04 to an infinite-dimensional analogue: is the universal moduli space of integrable systems, KP flows are the universal linear evolutions, and tau-functions are the universal sections of the determinantal line bundle.

Exercises Intermediate+

Lean formalization Intermediate+

Mathlib does not yet contain the algebra of formal pseudo-differential operators, the Sato Grassmannian, or the Hirota bilinear formalism. The unit is formalisation-free at the symbolic level; a meaningful Lean statement would require the entire upstream chain documented in Mathlib gap analysis. The theorem statement that would be the target, once the infrastructure exists, has the schematic form:

-- Aspirational, not currently realisable in Mathlib.
theorem sato_kp_correspondence
    (W : SatoGrassmannian) :
    let τ := W.tauFunction
    let u (t : ℕ → ℝ) := 2 * (∂² (Real.log ∘ τ) / ∂t 1 ^ 2)
    KPHierarchy.isFormalSolution u ∧
    (∀ n ≥ 1, ∂ (W.translate (n := n)) / ∂t n = LinearTranslation z^n W) :=
sorry

The statement requires SatoGrassmannian (as an infinite-dimensional ind-scheme with determinantal line bundle), tauFunction (as the Plücker section), KPHierarchy.isFormalSolution (Hirota bilinear equations in infinitely many variables), and pseudo-differential operator algebra PseudoDiffOp with positive-part splitting. None of this exists in current Mathlib. The closest existing infrastructure is the FormalMultilinearSeries machinery, which handles formal power series but not the Laurent / pseudo-differential extension. Tracked as a long-horizon contribution roadmap.

Advanced results Master

The Sato programme (1981). Sato (1981, RIMS Kokyuroku 439 [Sato 1981]) introduced the Sato Grassmannian and the tau-function as the central organising objects of soliton theory. The KP hierarchy is the orbit of a single vector under the affine action with , and the entire infinite tower of nonlinear KP equations reduces to the Plücker identities of the orbit. The originator paper is in Japanese; the standard English exposition is Sato-Sato 1983, with the canonical book-form treatment in Miwa-Jimbo-Date 2000 Solitons: Differential Equations, Symmetries and Infinite-Dimensional Algebras.

The DJKM transformation-group programme. Date, Jimbo, Kashiwara, and Miwa (1981-1983, RIMS Kokyuroku series and Publ. RIMS; collected in DJKM 1983, Publ. RIMS 19 [DJKM 1981]) reinterpreted the Sato Grassmannian via the boson-fermion correspondence and the action of the affine Lie algebra . Tau-functions become matrix elements for ; vertex operators on Fock space generate the symmetries of the KP hierarchy; the Hirota bilinear identities are Plücker identities for the affine Lie group orbit. Reductions to subgroups produce all the classical soliton hierarchies: KdV via , modified KP via the universal central extension, -KdV via , , Toda chains via further affine types.

The Segal-Wilson reformulation (1985). Segal and Wilson (1985, Publ. Math. IHES 61 [Segal-Wilson 1985]) recast the Sato Grassmannian via the loop group : is identified with the homogeneous space where denotes loops that extend holomorphically to the disc. The KP hierarchy is the action of the "positive Birkhoff" subgroup, and Krichever's algebro-geometric construction is realised as the embedding of the spectral curve into via the line-bundle-determinantal sections. This frames soliton theory inside the geometry of loop groups, opening the door to applications in conformal field theory (Frenkel-Lepowsky-Meurman vertex algebra constructions), in the geometric Langlands programme (the Grassmannian as the local affine Grassmannian of ), and in the geometry of moduli spaces of vector bundles on Riemann surfaces (Beauville-Laszlo 1994).

The Mulase-Shiota algebro-geometric programme. Mulase 1984 [Mulase 1984] (Adv. Math. 54) characterised the algebro-geometric points of as those stabilised by a finitely-generated commutative subalgebra of multiplication operators, equivalently the affine coordinate ring of a compact algebraic curve . This recovers Krichever's 1977 construction as a special case of the Sato picture. Shiota 1986 [Shiota 1986] (Invent. Math. 83) proved the converse Novikov conjecture: a ppav is the Jacobian of a curve iff its theta function is a KP tau-function. This solves the Schottky problem by identifying the Jacobian locus as the locus on which the theta function satisfies the bilinear KP equation — an algebraic characterisation via a single PDE in three variables. Refinements: Welters 1984 (Acta Math. 157) and Krichever 2010 (Funkt. Anal. Pril. 44) reduced Shiota's condition to the existence of a single trisecant line of the Kummer variety , the deepest known geometric characterisation of .

The Witten-Kontsevich theorem. Witten 1991 (Surveys in Differential Geometry 1 [Witten 1991]) conjectured, and Kontsevich 1992 (Comm. Math. Phys. 147 [Kontsevich 1992]) proved, that the generating function of -class intersection numbers $$ F(t_0, t_1, \ldots) = \sum_{g, n} \frac{1}{n!} \sum_{a_1, \ldots, a_n} \big\langle \tau_{a_1} \cdots \tau_{a_n}\big\rangle_g \prod_i t_{a_i}, $$ on the Deligne-Mumford moduli space of stable curves, is the unique tau-function of the KdV (= 2-KP) hierarchy fixed by the string equation. Kontsevich's proof realises as a Hermitian-matrix integral over Hermitian matrices, with the Feynman diagrams of the cubic matrix model identified with Strebel-Jenkins differentials triangulating . The theorem identifies the moduli space of curves as a tau-function object — moduli-of-curves data is integrable-systems data. The picture extends to: Mirzakhani 2007 (J. Amer. Math. Soc. 20) Weil-Petersson volumes as a virtual KdV tau-function via the recursion relations; Eynard-Orantin 2007 (Comm. Number Theory Phys. 1) topological recursion as a vast generalisation of the Sato picture to spectral-curve data with arbitrary branch structure; Givental 2001 (Mosc. Math. J. 1) cohomological field theories as towers of tau-functions over cohomology of ; Okounkov-Pandharipande 2006 (Ann. Math. 163) Hurwitz numbers as KP tau-functions via the ELSV formula.

Reductions and generalisations. The Drinfeld-Sokolov 1984 reduction (J. Soviet Math. 30) extracts the -KdV (Gelfand-Dickey) hierarchies from KP via the constraint ; the case is KdV, is Boussinesq, and the general case carries -symmetry. The 2-component KP / Toda hierarchy of Ueno-Takasaki 1984 (Adv. Stud. Pure Math. 4) extends to two pseudo-differential operators and includes the Toda chain, Ablowitz-Ladik lattice, and the discrete KdV equation. The matrix KP hierarchy (Sato 1981 generalisation) replaces with , producing the AKNS, mNLS, multi-component NLS systems. The supersymmetric extension (Manin-Radul 1985) replaces by a Grassmann algebra, producing the super-KP hierarchy underlying superstring perturbation theory. Quantum KP (Frenkel-Reshetikhin 1996) is the -deformation, with the quantum tau-function a -character on the quantum affine algebra .

Synthesis. The Sato programme identifies the KP hierarchy as a universal integrable system: every classical soliton equation is a reduction of KP, every algebraic curve is a point of the Sato Grassmannian, every moduli problem with a natural integrable structure has a tau-function, and the algebraic geometry of the Schottky problem is equivalent to integrable-systems data. The Witten-Kontsevich theorem is the apex statement of the picture in the moduli-of-curves direction; the Shiota theorem is the apex statement in the Schottky / abelian-varieties direction. Together they identify the moduli space , the Jacobian locus , and the universal integrable hierarchy on as three faces of a single object.

Full proof set Master

Lemma (formal residue Plücker identity). Let have wave function and adjoint wave function $\psi^W(t, z) = (1 + O(z^{-1})) e^{-\sum_n t_n z^n}W\mathrm{Gr}\infty$ iff* $$ \mathrm{Res}_{z = \infty}\big(\psi_W(t, z) \cdot \psi^*_W(t', z)\big) = 0 \quad \text{for all } t, t'. $$

Proof. The residue at of a Laurent series is the coefficient . The wave functions are normalised admissible bases of and respectively (where is the natural pairing induced by , ). The residue vanishes iff the two subspaces are dual in this pairing, iff and are simultaneously commensurable with via the same index — iff .

Proposition (Lax-Sato compatibility). Let be a first-order pseudo-differential operator with leading term . The Lax-Sato flows for pairwise commute.

Proof. By the Jacobi identity for commutators of operators, . Expanding via the Jacobi identity and using that , the cross terms cancel by virtue of the identity in the algebra (a consequence of the splitting and the Jacobi identity), and the surviving term vanishes by the same splitting argument applied to . The full check is in Sato 1981 / Dickey 2003 Soliton Equations and Hamiltonian Systems §1. The commuting flows define an infinite-dimensional integrable system on the orbit of .

Theorem (Sato 1981 — main correspondence, full statement). There is a bijection between:

  1. Points of index zero, modulo the action of multiplication by elements of ;
  2. Tau-functions satisfying the Hirota bilinear identities , modulo overall scalar;
  3. Formal solutions of the KP hierarchy with as .

Proof. The chain (1) (2) (3) is constructed in Steps 1-3 of the Key Theorem above. For the reverse:

(3) (2). Given a KP solution , integrate to recover up to a function of ; the higher KP flows determine the remaining freedom uniquely. The Hirota equations are then automatic from the KP hierarchy.

(2) (1). Given , define the wave function and the adjoint similarly. The Plücker residue identity (Lemma above) is equivalent to the Hirota identity on . The subspace is then a well-defined index-zero point of .

The bijection respects the KP-flow action: in (1) the flow is , in (2) it is the time translation , in (3) it is the actual KP evolution.

Corollary (linear flow on ). The KP hierarchy is a linear dynamical system on in the precise sense that the orbit of under the flow group is the integral submanifold of a commuting family of vector fields induced by multiplication by on .

Proof. Immediate from Step 2 of the Key Theorem and the abelian-ness of the flow group .

Connections Master

  • Finite-gap integration 05.09.09. The finite-gap KdV / KP solutions are the points of stabilised by a finitely-generated commutative subalgebra of multiplication operators on (Krichever-Mulase 1985). The Krichever map sends each compact algebraic curve with marked data to a Sato Grassmannian point, and the Its-Matveev formula is the special case of the Sato tau-function for these algebro-geometric subspaces.

  • Theta function 06.06.05. The Riemann theta function of a compact Riemann surface enters the Sato tau-function via the Krichever embedding: for finite-gap subspaces . The Shiota 1986 characterisation of the Jacobian locus via the KP equation is the converse: the theta function is a KP tau-function iff is a Jacobian.

  • Integrable system 05.02.03. The Sato Grassmannian is the universal moduli space of integrable systems. Every classical Liouville-Arnold-integrable Hamiltonian system on a finite-dimensional symplectic manifold has an infinite-dimensional analogue as a flow on via the spectral-curve construction (Hitchin 1987, Beauville-Narasimhan-Ramanan 1989).

  • Action-angle coordinates 05.02.04. The Sato action is the infinite-dimensional analogue of the Liouville-Arnold angle coordinates: the times are linear coordinates on the orbit through , and the orbit is a tower of finite-dimensional invariant tori (one per finite-gap reduction). The Lax-Sato evolution is the universal linearisation.

  • Schottky problem 06.06.06. The Shiota characterisation of the Jacobian locus by the KP equation identifies the Schottky problem with the integrability of a single PDE. Trisecant-line refinements (Welters 1984, Krichever 2010) give the deepest known geometric characterisation of .

  • Moduli of curves and intersection theory. The Witten-Kontsevich theorem identifies the generating function of -class intersection numbers on as a KdV tau-function. This realises moduli-space cohomology as integrable-systems data and underpins the topological recursion programme (Eynard-Orantin 2007).

Historical & philosophical context Master

Kadomtsev and Petviashvili 1970 [Kadomtsev-Petviashvili 1970] (Soviet Physics Doklady 15) introduced the KP equation as the leading-order correction to KdV when weak transverse modulation is allowed, originally to study stability of shallow-water solitary waves. Hirota 1971 [Hirota 1971] (Phys. Rev. Lett. 27) discovered the bilinear method, writing exact -soliton solutions of KdV as polynomial-exponential combinations and exposing the underlying quadratic structure that would later be understood as Plücker identities.

Sato 1981 [Sato 1981] (RIMS Kokyuroku 439) introduced the Sato Grassmannian and the tau-function in a series of lectures at Kyoto. The Japanese-language paper was disseminated through the RIMS Kokyuroku and translated piecewise into English; the canonical English exposition followed in Sato-Sato 1983 Soliton equations as dynamical systems on infinite-dimensional Grassmann manifolds (Lecture Notes in Num. Appl. Anal. 5). Date, Jimbo, Kashiwara, and Miwa 1981-1983 [DJKM 1981] (RIMS Kokyuroku series, collected in Publ. RIMS 19) developed the transformation-group / vertex-operator interpretation, identifying the KP hierarchy with the orbit of a single highest-weight vector under the affine Lie group and giving the boson-fermion correspondence as the dictionary between bosonic tau-functions and fermionic Fock-space matrix elements.

Mulase 1984 [Mulase 1984] (Adv. Math. 54) gave the algebro-geometric characterisation of finite-gap subspaces inside . Segal and Wilson 1985 [Segal-Wilson 1985] (Publ. Math. IHES 61) reformulated the Sato Grassmannian via loop groups, opening the door to applications in CFT and the geometric Langlands programme. Shiota 1986 [Shiota 1986] (Invent. Math. 83) proved Novikov's 1979 conjecture, characterising the Jacobian locus by the KP equation and solving the Schottky problem.

Witten 1991 [Witten 1991] (Surveys in Differential Geometry 1) conjectured the KdV-tau-function structure of -intersection-number generating functions, motivated by 2D quantum gravity and matrix models. Kontsevich 1992 [Kontsevich 1992] (Comm. Math. Phys. 147) proved the conjecture via Hermitian matrix integrals and Strebel-Jenkins ribbon graphs. The Witten-Kontsevich theorem is the apex statement linking moduli of curves to integrable systems, and the seed of the modern topological-recursion (Eynard-Orantin), Mirzakhani-Weil-Petersson, Givental cohomological-field-theory, and Okounkov-Pandharipande Gromov-Witten / Hurwitz programmes.

Bibliography Master

@article{KadomtsevPetviashvili1970,
  author = {Kadomtsev, B. B. and Petviashvili, V. I.},
  title = {On the stability of solitary waves in weakly dispersing media},
  journal = {Soviet Physics Doklady},
  volume = {15},
  year = {1970},
  pages = {539--541},
}

@article{Hirota1971,
  author = {Hirota, Ryogo},
  title = {Exact solution of the {K}orteweg--de {V}ries equation for multiple collisions of solitons},
  journal = {Physical Review Letters},
  volume = {27},
  year = {1971},
  pages = {1192--1194},
}

@article{Sato1981,
  author = {Sato, Mikio},
  title = {Soliton equations as dynamical systems on infinite-dimensional {G}rassmann manifolds},
  journal = {RIMS Kokyuroku},
  volume = {439},
  year = {1981},
  pages = {30--46},
}

@article{DJKM1983,
  author = {Date, Etsuro and Jimbo, Michio and Kashiwara, Masaki and Miwa, Tetsuji},
  title = {Transformation groups for soliton equations},
  journal = {Publications of RIMS, Kyoto University},
  volume = {19},
  year = {1983},
  pages = {943--1001},
}

@article{Mulase1984,
  author = {Mulase, Motohico},
  title = {Cohomological structure in soliton equations and {J}acobian varieties},
  journal = {Advances in Mathematics},
  volume = {54},
  year = {1984},
  pages = {57--66},
}

@article{SegalWilson1985,
  author = {Segal, Graeme and Wilson, George},
  title = {Loop groups and equations of {K}d{V} type},
  journal = {Publications Math{\'e}matiques de l'IH{\'E}S},
  volume = {61},
  year = {1985},
  pages = {5--65},
}

@article{Shiota1986,
  author = {Shiota, Takahiro},
  title = {Characterization of {J}acobian varieties in terms of soliton equations},
  journal = {Inventiones Mathematicae},
  volume = {83},
  year = {1986},
  pages = {333--382},
}

@article{Krichever1977,
  author = {Krichever, I. M.},
  title = {Methods of algebraic geometry in the theory of non-linear equations},
  journal = {Russian Mathematical Surveys},
  volume = {32},
  year = {1977},
  pages = {185--213},
}

@article{Witten1991,
  author = {Witten, Edward},
  title = {Two-dimensional gravity and intersection theory on moduli space},
  journal = {Surveys in Differential Geometry},
  volume = {1},
  year = {1991},
  pages = {243--310},
}

@article{Kontsevich1992,
  author = {Kontsevich, Maxim},
  title = {Intersection theory on the moduli space of curves and the matrix {A}iry function},
  journal = {Communications in Mathematical Physics},
  volume = {147},
  year = {1992},
  pages = {1--23},
}

@article{Welters1984,
  author = {Welters, Gerald},
  title = {A criterion for {J}acobi varieties},
  journal = {Annals of Mathematics},
  volume = {120},
  year = {1984},
  pages = {497--504},
}

@book{MiwaJimboDate2000,
  author = {Miwa, Tetsuji and Jimbo, Michio and Date, Etsuro},
  title = {Solitons: Differential Equations, Symmetries and Infinite-Dimensional Algebras},
  publisher = {Cambridge University Press},
  year = {2000},
}

@book{Dickey2003,
  author = {Dickey, L. A.},
  title = {Soliton Equations and {H}amiltonian Systems},
  publisher = {World Scientific},
  edition = {2nd},
  year = {2003},
}

@article{DrinfeldSokolov1984,
  author = {Drinfeld, V. G. and Sokolov, V. V.},
  title = {{L}ie algebras and equations of {K}orteweg--de {V}ries type},
  journal = {Journal of Soviet Mathematics},
  volume = {30},
  year = {1985},
  pages = {1975--2036},
}

@article{Mirzakhani2007,
  author = {Mirzakhani, Maryam},
  title = {Weil-{P}etersson volumes and intersection theory on the moduli space of curves},
  journal = {Journal of the American Mathematical Society},
  volume = {20},
  year = {2007},
  pages = {1--23},
}