Continuum Fractons: Quantization and the Few Body Problem

Technical Summary: Continuum Fractons: Quantization and the Few Body Problem

Problem Statement
This work addresses the canonical quantization of non-relativistic, dipole-conserving fractons in the continuum. While classical continuum fractons exhibit unique dynamical features—such as "Machian" dynamics (immobility of isolated particles), configuration space attractors, and strong ergodicity breaking—their quantum mechanical counterparts were previously unexplored. The central questions are: What are the energy spectra and eigenfunctions of the quantized Hamiltonian? How does dynamical evolution proceed, and do the quantum systems retain the classical attractors and broken ergodicity? The authors focus on systems of $N$ distinguishable particles, with detailed analysis for $N=2$ and $N=3$, treating strict fractons that are non-interacting beyond a finite separation.

Methodology
The authors formulate a continuum quantum mechanics starting from a classical Hamiltonian invariant under global dipole conservation ($D = \sum x_i$) and translation. The kinetic term involves a pair inertia function $K(|x_i - x_j|)$, which has compact support and vanishes at the boundaries of its support.

1. Canonical Quantization: The system is quantized by promoting coordinates and momenta to operators. The authors adopt a specific operator ordering, $\hat{H} = -\sum_{i<j} (\partial_{x_i} - \partial_{x_j}) K(|x_i - x_j|) (\partial_{x_i} - \partial_{x_j}) / 2$, to define the Hamiltonian.
2. Reduction of Degrees of Freedom: A unitary canonical transformation is applied to separate the conserved total dipole moment and total momentum, reducing the problem to $N-1$ degrees of freedom.
3. Sturm-Liouville (SL) Theory: For the two-body problem, the Hamiltonian is mapped to a Sturm-Liouville operator on a finite interval. The spectral properties are analyzed based on the behavior of $K(x)$ near the edges of its support, specifically $K(x) \sim |x - x_{edge}|^\theta$.
4. Liouville Transform: To understand the dynamics and spectral transition, the authors employ a Liouville transformation to map the SL operator to an equivalent Schrödinger operator with an effective potential $V(\xi)$. This allows the use of standard scattering theory and WKB approximations.
5. Lattice Regularization for $N=3$: For the three-body problem, the differential operator becomes piecewise defined with internal matching conditions that lack a standard SL classification. To resolve this, the authors introduce a lattice regularization that preserves dipole conservation. This construction implicitly selects a specific set of matching conditions for the continuum limit, which are then analyzed numerically.

Key Contributions and Results

• Zero Modes and Degeneracy: The authors identify two sources of infinite zero modes and degeneracy. First, the conserved total dipole and momentum generate an infinite-dimensional zero-energy space. Second, for compact support $K(x)$, any square-integrable wavefunction outside the support of $K$ is an exact zero-energy eigenstate. The analysis focuses on the finite-energy spectrum within the support of $K$.

• Two-Body Spectral Transition ($\theta = 2$):

  • The spectrum undergoes a sharp transition controlled by the exponent $\theta$ describing how $K(x)$ vanishes at the edges.
  • $\theta < 2$ (Discrete Spectrum): The operator is in the limit-circle case. It requires self-adjoint boundary conditions (e.g., vanishing current at boundaries). The spectrum is discrete, and wavepackets reflect off the edges, behaving like bound states.
  • $\theta > 2$ (Continuous Spectrum): The operator is in the limit-point case, possessing a unique self-adjoint extension without imposed boundary conditions. The spectrum is continuous. Crucially, wavepackets do not reflect; instead, they slow down and "pile up" (squeeze) into asymptotically narrow regions at the edges where $K \to 0$.
  • Semiclassical Connection: This transition mirrors the classical "freezing time" $\tau_f$. For $\theta < 2$, $\tau_f$ is finite (bounded motion); for $\theta > 2$, $\tau_f$ diverges (unbounded motion in time), corresponding to the continuous spectrum.

• Three-Body Dynamics and Numerical Evidence:

  • The three-body Hilbert space is block-diagonal, separating into sectors corresponding to isolated pairs and a central "genuinely three-particle" hexagonal region.
  • Localization: Low-energy eigenstates are found to localize along the classical attractor lines (regions where particles are clustered), consistent with the classical limit.
  • Tunneling: Unlike classical trajectories which settle into a single attractor region, quantum wavefunctions tunnel between the six symmetry-related attractor regions (permutations of particles), organizing into singlets and doublets of the Dihedral symmetry group.
  • Spectral Transition: Numerical scaling of energy gaps suggests a spectral transition in the three-body system similar to the two-body case, with a critical value $\theta_c$ flowing to 2 in the continuum limit.
  • Dynamics: For $\theta > \theta_c$, wavepacket evolution is significantly slower and shows partial permeation into attractor regions, contrasting with the rapid settling seen for $\theta < \theta_c$.

• Order-by-Disorder: The authors propose a "quantum order-by-disorder" mechanism where quantum fluctuations select the subset of the ground state manifold with the softest orthogonal directions. This explains the preference for low-energy states to reside in regions where the number of active kinetic terms (non-zero $K$) is minimized (Region 1), maximizing flat directions.

Significance and Claims
The paper claims to establish that the lack of ergodicity observed in classical continuum fractons survives quantization for large systems. The spectral transition at $\theta=2$ is identified as a fundamental feature distinguishing regimes of quantum fracton dynamics:

1. Dynamical Behavior: The transition separates regimes of wavepacket reflection (discrete spectrum) from regimes of edge accumulation (continuous spectrum), the latter being a quantum analog of classical attractors.
2. Validity of Continuum Limits: The work argues against the naive ultra-local limit ($K(x) \to \delta(x)$) often used in field theory formulations of fractons. The authors demonstrate that the spectral properties and clustering behavior are critically sensitive to the edge behavior of the inertia function $K(x)$, which is lost in the ultra-local limit.
3. Mathematical Structure: The paper highlights that while the two-body problem is solvable via standard Sturm-Liouville theory, the $N \ge 3$ problem introduces complex matching conditions and Hilbert space fragmentation (Krylov sectors) that require non-perturbative approaches, such as the lattice regularization employed here.

The authors conclude that while the two-body problem is essentially solved, the many-body regime at finite density remains an open question regarding whether quantum effects restore ergodicity, though their "order-by-disorder" argument suggests ergodicity breaking persists at low temperatures.