Exact Phase-Space Rotation in the Trapped Quantum Calogero Model

Technical Summary: Exact Phase-Space Rotation in the Trapped Quantum Calogero Model

Problem Statement
The construction of microscopically exact transport equations for interacting quantum many-body systems remains a significant open challenge, closely linked to Hilbert's Sixth Problem for quantum systems. While Generalized Hydrodynamics (GHD) provides an empirically exact description of Euler-scale dynamics for integrable systems, a fully first-principles microscopic derivation of the underlying transport equation for interacting quantum models has historically been lacking. Previous work established such a derivation for the untrapped (free-streaming) quantum Calogero model using a Hermitian Wigner operator constructed from quantum Lax matrices. The central problem addressed in this paper is whether this microscopic strategy survives the introduction of an external harmonic confining potential, a setting that is generic for experimental realizations of one-dimensional quantum gases.

Methodology
The author extends the microscopic framework developed for the untrapped Calogero model to the harmonically trapped case. The methodology proceeds as follows:

1. Model Definition: The study considers the trapped Calogero Hamiltonian, which includes kinetic energy, an inverse-square interaction term, and a harmonic confinement term with frequency $\Omega$.
2. Quantum Lax Pair: The construction utilizes the known quantum Lax pair for the trapped model. This involves $N \times N$ operator-valued matrices $\hat{X}$ (position) and $\hat{L}$ (momentum-like), along with a symmetric matrix $\hat{A}$. The time evolution of these matrices is governed by the Lax relations:
   • $\frac{d}{dt}\hat{X} + i[\hat{X}, \hat{A}] = \hat{L}$
   • $\frac{d}{dt}\hat{L} + i[\hat{L}, \hat{A}] = -\Omega^2 \hat{X}$
     The presence of the harmonic trap modifies the second relation compared to the untrapped case, introducing the restoring force term $-\Omega^2 \hat{X}$.
3. Wigner Operator Construction: A Hermitian Wigner operator, $\hat{W}(x, t, \lambda)$, is defined using the quantum Lax matrices. It is constructed via a two-dimensional delta function of the operators $\hat{X}$ and $\hat{L}$, symmetrized (Weyl ordered) to ensure hermiticity. The operator is expanded in a multivariate Taylor series involving symmetrized products of $\hat{X}$ and $\hat{L}$.
4. Moment Recursion: The time evolution of the operator is analyzed by differentiating the symmetrized moments $O_{m,n}(t) = v^T S[\hat{X}^m \hat{L}^n] v$. By applying the product rule and utilizing the Lax equations, the commutator terms involving $\hat{A}$ vanish due to the properties of the vector $v$ (where $\hat{A}v = 0$). This leaves a recursion relation for the moments driven solely by the source terms from the Lax equations.
5. Derivation of the Transport Equation: The derived recursion relation for the moments is substituted back into the time derivative of the Wigner operator. Through distributional identities involving derivatives of delta functions, the infinite series of moment derivatives is re-summed to yield a closed partial differential equation.

Key Contributions and Results
The paper establishes the following exact results:

• Exact Microscopic Transport Equation: The Wigner operator $\hat{W}(x, t, \lambda)$ and its expectation value $\rho(x, t, \lambda)$ satisfy the exact operator identity:
  $$\partial_t \hat{W} + \lambda \partial_x \hat{W} - \Omega^2 x \partial_\lambda \hat{W} = 0$$
  This equation holds for arbitrary particle number $N$, arbitrary interaction strength $\ell$, and arbitrary initial density matrices, without invoking any scaling limits (such as large $N$ or low energy).
• Rigid Phase-Space Rotation: The resulting dynamics describe a rigid rotation in the phase space $(x, \lambda)$ with angular frequency $\Omega$. The solution for the phase-space density is given by:
  $$\rho(x, t, \lambda) = \rho\left(x \cos \Omega t - \frac{\lambda}{\Omega} \sin \Omega t, 0, \lambda \cos \Omega t + \Omega x \sin \Omega t\right)$$
  This confirms that the system exhibits isochronous dynamics with a period of $2\pi/\Omega$ for any initial state.
• Behavior of Conserved Charges: In the untrapped model, moments of rapidity are individually conserved. In the trapped model, individual moments $O_{0,n}$ are no longer conserved. Instead, they form rotating multiplets. The paper demonstrates that within the quadratic sector, the unique conserved combination is proportional to the trapped Hamiltonian $\hat{H}_\Omega$, specifically $O_{0,2} + \Omega^2 O_{2,0}$. This serves as a nontrivial consistency check, as energy conservation was not assumed in the derivation but emerged naturally.
• Limiting Cases: The derived equation reduces exactly to the free-streaming equation of the untrapped model in the limit $\Omega \to 0$.

Significance
The paper claims to provide a microscopic realization of the isochronous dynamics previously known for trapped Calogero-type systems, extending this understanding to arbitrary initial states and arbitrary quantum entropy densities. By constructing a Hermitian Wigner operator that obeys an exact phase-space rotation equation, the work closes the gap between generalized hydrodynamics and first-principles microscopic derivations for interacting quantum systems in confined geometries. The result generalizes the periodic "breathing" solutions of classical trapped Calogero fluids to the full quantum regime. Furthermore, the framework suggests a potential generalization to arbitrary external potentials $U(x)$, where the phase-space equation would take the form $\partial_t \rho + \lambda \partial_x \rho - U'(x) \partial_\lambda \rho = 0$, though the exact solvability via rotation is noted as a special feature of the harmonic case.