Hamiltonian formulation of the supersymmetric KdV equation

This paper presents a constrained Hamiltonian formulation of a specific supersymmetric KdV equation (with parameter $a=2$) using the Dirac-Bergmann algorithm, revealing a unique nonlocal contribution to the Hamiltonian density arising from fermionic consistency conditions and demonstrating the equivalence between the derived Hamiltonian dynamics and the component-form supersymmetric KdV system.

The Gist

Imagine the universe is filled with waves. Some waves are simple, like ripples in a pond. Others are complex, like the massive, crashing waves of a storm that seem to keep their shape and energy forever. In physics, there is a famous mathematical rule called the KdV equation that describes these special, shape-shifting waves (called solitons). It's like a perfect recipe for how these waves move and interact.

Now, physicists love to add "flavor" to their recipes. One popular flavor is supersymmetry. Think of this as adding a new type of ingredient to the mix: "fermionic" fields. If the original wave ingredients are like solid, predictable bricks (bosons), these new fermionic ingredients are like ghostly, jittery particles that behave very differently. When you mix them together, you get a "Supersymmetric KdV" system. It's a more complex, ghostly version of the original wave equation.

The Problem: A Broken Scale
The authors of this paper wanted to study this ghostly wave system using a specific tool called Hamiltonian mechanics. You can think of Hamiltonian mechanics as a precise scale or a set of rules that tells you exactly how a system will move forward in time based on its current energy and position.

Usually, to use this scale, you need a "Lagrangian" (a master formula describing the system's energy). However, for this specific supersymmetric wave system (where a parameter called $a$ equals 2), the Lagrangian is degenerate.

Here is an analogy: Imagine trying to weigh a cloud on a bathroom scale. The scale is designed for solid objects. If you put a cloud on it, the needle doesn't just point to a number; it spins wildly or gets stuck because the cloud doesn't have a single, solid "weight" in the traditional sense. The math breaks down because the relationship between the "speed" of the fields and their "momentum" is tangled and unclear. You can't just flip a switch (a standard mathematical operation called a Legendre transformation) to get the answer.

The Solution: The Detective's Algorithm
Since the standard scale didn't work, the authors used a special detective method called the Dirac-Bergmann algorithm.

1. Finding the Clues (Constraints): Instead of forcing the system to fit the standard rules, they looked for "clues" or constraints. These are hidden rules that the system must obey. For example, they found that two of the ghostly fields (let's call them $\psi$ and $\xi$) aren't actually independent; they are locked together in a specific way ($\psi = \xi_x$). It's like realizing that in your ghostly wave recipe, the amount of "jitter" is strictly determined by the slope of the wave itself. You can't change one without changing the other.
2. Building the New Scale: Once they identified these locked relationships (constraints), they built a new, custom "Total Hamiltonian" (the master rulebook for the system's motion). This new rulebook included special "penalties" (Lagrange multipliers) to ensure the system never broke those locked rules.

The Surprise: The "Non-Local" Ghost
The most interesting discovery in their new rulebook was a non-local term.

In everyday life, things are usually "local." If you push the left end of a rope, the right end doesn't move instantly; the wave travels down the rope. But in this mathematical system, the authors found that the energy of the wave at one point depends on the entire history or the whole shape of the wave elsewhere.

Think of it like a telepathic wave. If you tweak the wave at point A, the energy calculation at point B instantly changes, not because a signal traveled, but because the system's internal rules link the two points together in a way that requires "looking at the whole picture" to understand the energy. This happens because the fundamental building block of the system isn't the field itself, but its derivative (its slope), which forces the math to reach out and grab information from far away.

The Result: Two Ways to See the Same Thing
The authors showed that their new, complex rulebook (the component-level Hamiltonian) works perfectly. When they ran the equations, the waves moved exactly as the original supersymmetric KdV system predicted.

They also translated this complex, component-by-component rulebook into a Superspace version. Think of "Superspace" as a higher-dimensional map where the "ghostly" and "solid" ingredients are blended into a single, smooth object. In this high-level view, the weird "telepathic" (non-local) connection disappears, and the rulebook looks clean and compact again. This proved that their messy, detailed math and the elegant, high-level math were describing the exact same reality.

In Summary
The paper is a guide on how to weigh a "cloud" (a complex supersymmetric wave system) when your standard scale is broken. By using a detective's algorithm to find hidden rules (constraints) and building a custom scale, the authors successfully described how these ghostly waves move. They discovered that to understand the energy of these waves, you sometimes have to look at the whole system at once (non-locality), but if you view the system from a higher, more abstract angle (superspace), the mystery simplifies into a beautiful, compact formula.

Technical Summary

Technical Summary: Hamiltonian Formulation of the Supersymmetric KdV Equation

Problem Statement
The paper addresses the constrained Hamiltonian formulation of the supersymmetric Korteweg–de Vries (sKdV) equation, specifically the case characterized by the parameter $a=2$ (sKdV-2). While the classical KdV equation is a well-understood completely integrable system with a rich multi-Hamiltonian structure, its supersymmetric extensions present unique challenges. The authors note that for the sKdV-$a$ family, a Lagrangian formulation exists only for specific values of $a$. The case $a=0$ is trivial, reducing to the standard KdV equation. However, the case $a=2$ is nontrivial and possesses a consistent Lagrangian. The central problem is that this specific Lagrangian is degenerate (singular) because it is linear in the time derivatives of all fields. Consequently, the standard Legendre transformation fails to yield a Hamiltonian, necessitating a systematic constrained analysis to determine the system's dynamics and constraints.

Methodology
The authors employ the Dirac–Bergmann algorithm (DBA) to construct the Hamiltonian formulation for the sKdV-2 system. The methodology proceeds as follows:

1. Lagrangian Construction: Starting from a nontrivial Lagrangian density ($\mathcal{L}_{sKdV-2}$) involving bosonic ($u$) and fermionic ($\psi, \xi$) component fields, the authors identify the system as degenerate due to the vanishing determinant of the Hessian matrix with respect to time derivatives.
2. Primary Constraints: Using the left derivative convention for the Legendre transformation, the canonical momenta are calculated. Because the velocities cannot be uniquely solved for in terms of momenta, primary constraints are identified:
   • $c_1 = \Pi_u + \frac{1}{2}u_x$
   • $c_2 = \Pi_\psi + \frac{1}{2}\psi$
   • $c_3 = \Pi_\xi$
3. Constraint Consistency and Secondary Constraints: The consistency of these constraints over time (requiring their Poisson brackets with the total Hamiltonian to vanish) is analyzed. This process reveals a secondary constraint, $\tilde{c}_1 = \psi - \xi_x = 0$, which establishes a correlation between the fermionic degrees of freedom.
4. Determination of Multipliers: The consistency conditions are used to solve for the Lagrange multipliers ($\lambda_i$) associated with the constraints. Notably, the determination of the multiplier $\lambda_3$ involves an inverse derivative operator ($\partial_x^{-1}$), introducing nonlocality into the system.
5. Hamiltonian Construction: The total Hamiltonian is constructed by summing the canonical Hamiltonian density and the contributions from all primary and secondary constraints weighted by their respective multipliers.

Key Results

• Total Hamiltonian Density: The authors explicitly derive the total Hamiltonian density, which includes a nonlocal term arising from the inverse derivative operator in the Lagrange multiplier $\lambda_3$. This nonlocality is a direct consequence of the constrained dynamics and the fact that the fundamental physical variable is identified as $\xi_x$ rather than $\xi$.
• Equations of Motion: The derived Hamilton equations of motion are shown to reproduce the component form of the sKdV-2 system (equations 7 and 8 in the text), confirming the validity of the formulation.
• Superspace Representation: The authors demonstrate that the component-level Hamiltonian can be expressed in a compact superspace form:
  $$\bar{H} = \frac{1}{2} \int \left[ -2D^2\Phi (D^3\Phi)^2 + D^4\Phi D^5\Phi \right] dx d\theta$$
  In this superspace representation, the nonlocality observed in the component formulation disappears, as the relations between the fermionic component fields ($\psi = \xi_x$) effectively remove the inverse derivative operator.
• Consistency Checks: The superspace Hamiltonian is shown to be consistent with previous studies on the sKdV-2 system (specifically referencing Mathieu [42]), up to the redefinition of the superfield necessitated by the introduction of the velocity potential ($u \equiv u_x$). Furthermore, in the limit where the fermionic component vanishes ($\xi \to 0$), the Hamiltonian correctly reduces to the standard KdV Hamiltonian.

Significance and Claims
The paper claims to provide a complete constrained Hamiltonian formulation for the sKdV-2 system, a case that had not been fully detailed in this manner previously. The significance of the work lies in:

1. Handling Degeneracy: Successfully applying the Dirac–Bergmann algorithm to a supersymmetric system with a degenerate Lagrangian, explicitly identifying the full set of constraints.
2. Revealing Nonlocality: Highlighting a distinctive feature of this supersymmetric extension: the emergence of nonlocal contributions to the Hamiltonian density at the component level, which are intrinsic to the constrained nature of the system.
3. Unifying Formulations: Establishing a bridge between component-level and superspace formulations, demonstrating how the nonlocality in the component description is resolved in the superspace representation.
4. Methodological Extension: The authors suggest that the approach presented can serve as a template for investigating the constrained Hamiltonian structures of other supersymmetric integrable systems, such as the supersymmetric nonlinear Schrödinger equation (NLSE), as well as nonlocal nonlinear integrable PDEs and equations with delay.

The work is dedicated to the memory of Yavuz Nutku, acknowledging his foundational contributions to the Hamiltonian formulation of integrable systems.