«When I was finishing this article, one question wouldn't leave me: is the universality of the KP hierarchy truly a fundamental property of nature, or just a lucky mathematical abstraction? We showed that AKNS is embedded in KP through wave functions, but it remains unclear why such distinct physical systems obey the same structure. Perhaps we are seeing merely a projection of a deeper geometry that we do not yet know how to formulate.» – Dr. Daniel Stern
What are we trying to understand?
Imagine you are looking at a subway map with tangled transfers between lines. Some stations are connected directly; others require a long journey with multiple changes. Now, imagine someone found a hidden tunnel that instantly connects two seemingly distant stations. This is roughly what mathematical physics does when it studies the connections between various integrable systems.
The Ablowitz-Kaup-Newell-Segur hierarchy (or, more accurately, Ablowitz-Kaup-Newell-Segur) \u2013 sounds intimidating, doesn't it? Let's simply call it AKNS. This is a family of equations that describe waves in nonlinear media: from light propagation in an optical fiber to the behavior of waves on a water surface. These equations are not just an abstraction; they describe real physical processes that can be measured and observed.
But here is what is interesting: there exists another, more general system \u2013 the Kadomtsev-Petviashvili, or KP, hierarchy. It plays the role of a sort of \u00abmother tongue\u00bb for many integrable systems. It has long been known that AKNS can be obtained from KP through certain restrictions, as if you took a three-dimensional shape and looked at its two-dimensional shadow. But the question remained: can this be shown directly, through the fundamental objects of the theory \u2013 the \u03c4-functions?
This is exactly the question explored in the work I want to tell you about. And the answer it gives turns out to be more elegant than one might have expected.
What is a τ-function and why is it important?
What is a \u03c4-function and why is it important?
If a hierarchy of equations is a subway map, then the \u03c4-function is the schedule for all trains simultaneously. It contains complete information about the solution of the entire system of equations. If you know this \u03c4-function \u2013 you know everything about the system's behavior.
The logarithmic derivatives of the \u03c4-function give us the so-called correlation functions \u2013 quantities that show how field values are connected at different points in space and time. It is as if you threw a stone into a pond and tried to understand how a ripple at one point affects waves at another.
To calculate these derivatives, mathematicians developed the method of matrix resolvent analysis. The resolvent is, roughly speaking, a way to \u00abinvert\u00bb a differential operator, turning it into something one can work with algebraically. Imagine you have a function that transforms an input signal into an output, and the resolvent allows you to recover the input signal from the output.
But this method had one drawback: it was quite abstract and did not always yield intuitively clear results. The new approach brings wave functions into play \u2013 objects that physicists have been accustomed to using since the days of quantum mechanics.
The AKNS Hierarchy: what is it all about?
Let's get specific. The AKNS hierarchy is built around a so-called Lax pair \u2013 two operators connected in a special way. The first operator describes the spatial structure of the system, the second \u2013 its evolution in time.
The key object here is a two-by-two matrix containing two fields, usually denoted by the letters q and r. These fields can describe, for instance, the amplitudes of counter-propagating waves in a nonlinear medium. The entire hierarchy of equations arises from the requirement that the spatial and temporal operators \u00abcommute\u00bb \u2013 that is, the order of their application does not matter.
This commutation condition generates an infinite chain of related equations. The first terms of this chain give us the classical equations of mathematical physics: the Korteweg-de Vries equation describing solitary waves on shallow water, the nonlinear Schrödinger equation used in optics and condensed matter physics, and others.
What makes these equations special? They are exactly solvable. Unlike most nonlinear equations, which can only be solved approximately or numerically, the equations of the AKNS hierarchy admit analytical solutions. Moreover, these solutions possess amazing stability properties: solitons \u2013 localized wave packets \u2013 can collide with each other and separate while maintaining their shape.
Wave functions: the connecting link
The new approach introduces a pair of wave functions that depend on the spatial coordinate, time variables, and a special parameter called the spectral parameter. This parameter can be understood as \u00abfrequency\u00bb or \u00abenergy\u00bb \u2013 in different contexts it has a different physical sense.
The wave function in this case is not a quantum-mechanical probability amplitude, although mathematically the objects are similar. It is a solution to a linear system of equations connected to the original nonlinear hierarchy. The key word here is \u00ablinear\u00bb. One of the main ideas of integrability theory is that complex nonlinear systems can be reduced to simpler linear ones, but at the cost of introducing additional parameters and spaces.
Two wave functions, let's denote them \u03a8 and \u03a8\u0304, are linked by the system's symmetry. One grows in the direction of positive spatial coordinate values, the other \u2013 negative. Together, they form a complete basis of solutions for the linear problem.
These functions satisfy a system of equations that includes both spatial and temporal evolution. Asymptotically, at large distances, they behave like exponents with the spectral parameter in the exponent. This asymptotic behavior is critical for building the theory \u2013 it sets the normalization and determines exactly how the wave functions \u00absee\u00bb the structure of the solution.
From wave functions to the resolvent
The matrix resolvent is the Green's function for a differential operator. If the operator describes how the system responds to an external influence, then the Green's function shows the response at a point in space to an impulse influence at another point. In electrostatics, this would be the field of a point charge; in elasticity theory \u2013 the displacement from a point force.
Usually, the resolvent is calculated using abstract methods of functional analysis. But if we have explicit wave functions, we can write down the resolvent constructively: it is built from products of wave functions at two different points in space.
The formula looks like this: if the first point is to the right of the second, the resolvent is given by the product of the first wave function at the first point and the conjugate second wave function at the second point. If the order of points is reversed \u2013 the wave functions swap places. At the boundary, where the points coincide, a jump occurs, which defines the delta function in the equation for the resolvent.
This is not just a technical trick. Such a representation shows that the resolvent is not an abstract object, but something built from physically meaningful wave functions. Each wave function carries information about how the system propagates disturbances in a certain direction.
Correlation functions: what do they measure?
Now let's move on to correlation functions. Imagine you are studying the surface of the ocean. You can measure the wave height at one point \u2013 this gives you local information. But if you want to understand the structure of the swell, you need to know how the height at one point relates to the height at a second, third, or fourth point. These connections are described by correlation functions.
In mathematical theory, a k-point correlation function is the derivative of the logarithm of the \u03c4-function with respect to k different parameters. These parameters can be time variables or spectral parameters, depending on the problem statement.
A classical result obtained by Okounkov and Reshetikhin in the early 2000s expressed these functions through traces of products of resolvents. The trace is the sum of the diagonal elements of a matrix, an invariant under many transformations. The formula involved summation over permutations \u2013 all possible ways to connect k points into a cyclic chain.
The new approach rewrites this formula, replacing abstract resolvents with products of wave functions. Instead of traces of resolvents, we get integrals of products of wave functions over spectral parameters. This makes calculations more explicit and often simpler, especially when the wave functions have good analytical properties.
Details of the new formula
Let's go into a bit more detail. The new formula for the k-point correlation function represents a sum over all permutations of k elements. Each term of the sum includes the sign of the permutation \u2013 plus for even permutations, minus for odd ones. This ensures the antisymmetry characteristic of fermionic systems, although the AKNS hierarchy itself does not necessarily describe fermions.
Each permutation contributes a term in the form of a product of traces. Each trace is taken from the product of a wave function at one point with a conjugate wave function at another point, with the spectral parameters linked via the permutation. Then integrals are taken over all spectral parameters along certain contours in the complex plane.
The choice of integration contours is critical. They must encompass singularities related to the operator's spectrum \u2013 poles, cuts, or other singularities of the resolvent. Depending on the specific solution of the AKNS hierarchy, these singularities can be discrete (like eigenvalues) or continuous (like spectral branches).
The advantage of such a formulation is that it clearly demonstrates the role of wave functions as fundamental dynamic objects. If for a specific solution the wave functions are known explicitly \u2013 and for many classical solutions they are \u2013 the correlation functions can be calculated directly, without resorting to abstract operator methods.
Connection with the Kadomtsev-Petviashvili hierarchy
Now for the most interesting part: how does all this relate to the KP hierarchy? The Kadomtsev-Petviashvili hierarchy is, in a sense, the \u00abmother of all integrable hierarchies\u00bb. Many famous systems are obtained from it by imposing constraints \u2013 reductions.
The \u03c4-function of the KP hierarchy is defined through Hirota bilinear equations or, identically, through infinite-dimensional determinants. It depends on an infinite set of time variables, and its logarithmic derivatives satisfy strict relations called Plücker equations or Hirota relations.
The AKNS hierarchy is obtained from KP by imposing certain conditions on the time variables. These conditions isolate a subset of solutions possessing additional symmetry. Physically, this corresponds to the transition from waves propagating in two spatial dimensions to waves in one dimension.
The key question: if we take an arbitrary solution of the AKNS hierarchy with its \u03c4-function, can we show that this \u03c4-function is a KP \u03c4-function? That is, does it satisfy all the structural requirements characteristic of KP?
Using the new formula for correlation functions via wave functions and the connection of these wave functions with Baker-Akhiezer wave functions, one can show that the answer is positive. The idea is to establish a link between the AKNS wave functions and the Baker-Akhiezer wave functions, which play a central role in KP theory.
Baker-Akhiezer wave functions
Baker-Akhiezer wave functions are special meromorphic functions on an algebraic curve (or its generalizations) satisfying certain analytical conditions. For the KP hierarchy, they depend on a point on the curve (which parameterizes the spectral parameter) and on an infinite set of time variables.
These functions have a characteristic asymptotic behavior: they behave like an exponent of a linear combination of time variables multiplied by a slowly varying function. The coefficients in the exponent are related to the local parameter on the curve near a special point \u2013 usually the point at infinity.
It turns out that the AKNS wave functions introduced in the discussed work can be interpreted as restrictions of the Baker-Akhiezer wave functions to a specific subspace of time variables. This restriction corresponds to the reduction from KP to AKNS.
Technically, this means the following: we take a KP wave function, which depends on the full set of times, and impose conditions that some combinations of these times are fixed or linked by certain relations. The resulting function satisfies the equations of the AKNS hierarchy and coincides with the wave functions defined through the Lax pair of this hierarchy.
Proof through correlation functions
To show that the AKNS \u03c4-function is a KP \u03c4-function, one needs to verify that the AKNS correlation functions satisfy the defining relations for KP. These relations include the Plücker equations \u2013 polylinear identities that the derivatives of the logarithm of the \u03c4-function must satisfy.
Using the expression of AKNS correlation functions via wave functions and the link between these wave functions and Baker-Akhiezer wave functions, one can show that the necessary identities are satisfied automatically. This is a non-trivial check \u2013 it requires careful handling of wave function asymptotics and integral representations.
But in the end, it turns out that the structure embedded in the AKNS wave functions already contains all the necessary information about the KP structure. The AKNS \u03c4-function is not just \u00absimilar\u00bb to the KP \u03c4-function \u2013 it literally is a KP \u03c4-function under a special choice of dependence on time variables.
Why is this important?
This result is important for several reasons. First, it establishes a direct and constructive link between two fundamental hierarchies of integrable systems. We knew such a link existed on a formal level, but now we have an explicit construction via wave functions.
Second, it provides a new tool for calculating correlation functions and other characteristics of AKNS hierarchy solutions. If we know the wave functions \u2013 and for many physically interesting solutions we do \u2013 we can efficiently calculate correlations of any order.
Third, it deepens our understanding of the universality of the KP hierarchy. It plays the role of a sort of \u00absource code\u00bb from which more specialized systems are derived through various reductions. Each such reduction inherits the rich structure of KP but manifests it in a simpler and more physically transparent form.
Finally, this is an example of how a purely mathematical construction \u2013 wave functions introduced almost formally \u2013 turns out to be the key to understanding deep connections between different physical theories. Mathematics here is not just a language of description; it is a tool of discovery.
Open questions and perspectives
As always in science, every answer spawns new questions. Can this approach be generalized to other reductions of the KP hierarchy? Do similar constructions exist for discrete integrable systems, where spatial and temporal variables take discrete values?
How does this formalism relate to quantum deformations of integrable systems? In quantum field theory, correlation functions acquire an additional meaning related to probability amplitudes. Is the structure discovered for classical hierarchies preserved when moving to the quantum case?
Are there applications of these results to specific physical systems? The nonlinear Schrödinger equation, for example, describes light propagation in optical fibers and the dynamics of Bose-Einstein condensates. Can the new formulas for correlation functions be used to predict observable effects in these systems?
And a more philosophical question: why are such different physical systems described by the same mathematical structure? What do waves on water, light pulses in a fiber, and quantum fluctuations of a condensate have in common such that they all obey the laws of a single hierarchy of equations? Is the universality of integrable systems a coincidence or a reflection of deep principles we do not yet fully understand?
Final thoughts
The work we discussed behaves like a typical example of modern mathematical physics. There are no loud claims of revolutionary discoveries, no promises to change the world. This is painstaking research into the links between mathematical structures, testing hypotheses, and deriving new formulas.
But it is precisely such work that forms the foundation of our understanding of nature. Every established link, every new formula is a brick in the building of theory. Today we learned that AKNS wave functions are connected to the resolvent in a certain way. Tomorrow this knowledge will help someone solve a specific problem in nonlinear optics or hydrodynamics. The day after tomorrow, it will become part of a textbook studied by future generations of physicists.
The beauty of mathematical physics lies in the fact that abstract constructions turn out to be applicable to the real world. \u03c4-functions and wave functions, correlation functions and resolvents \u2013 all of this is not a game of symbols on paper. It is a description of how real waves behave, how real light propagates, and how real quantum systems evolve.
And every time we find a new connection between mathematical objects, we learn something new about the world. Or, more precisely, about the language in which the laws of nature are written. Here is what we know now. Here is what remains a mystery. And here is why this is important.