Solvable Random Unitary Dynamics in a Disordered Tomonaga-Luttinger Liquid

This paper presents the first analytical derivation of the frame potential for a disordered Tomonaga-Luttinger liquid, revealing power-law decay and late-time saturation governed by a single coupling parameter, with specific applications to random field XXZ spin chains that offer direct insights for quantum algorithm design.

The Gist

Imagine you have a long, one-dimensional string of beads (a quantum system) that is constantly being shaken by a random, jittery hand (disorder). In physics, we usually study how this shaking affects specific things, like how a wave moves along the string. But this paper asks a different question: How "random" does the entire system become over time?

To answer this, the authors use a tool called the Frame Potential. Think of this as a "chaos meter."

• If the meter reads 1, the system is perfectly ordered and predictable (like a metronome).
• If the meter drops toward 0, the system has become maximally random, like a shuffled deck of cards where every outcome is equally likely.

Here is the story of what they found, broken down into simple concepts:

1. The Setup: A Noisy Quantum String

The scientists looked at a specific type of quantum system called a Tomonaga-Luttinger Liquid (TLL). You can imagine this as a very special, one-dimensional highway where particles (like electrons or atoms) move together in a coordinated dance.

• The Disorder: They added "quenched Gaussian forward-scattering disorder." In plain English, this means they sprinkled the highway with random, static bumps that only push the particles forward or backward slightly, but don't knock them off the road entirely.
• The Goal: They wanted to calculate exactly how fast the "chaos meter" (Frame Potential) drops as the system evolves.

2. The Big Breakthrough: A Perfectly Solvable Puzzle

Usually, calculating randomness in these messy, interacting systems is a nightmare. It's like trying to predict the exact path of every leaf in a storm while they are all bumping into each other.

• The Trick: The authors found a special case where the math works out perfectly. Because the disorder only pushes the particles in a specific way (forward scattering), the messy equations simplify into a neat, solvable shape (a "quadratic" structure).
• The Result: They derived a closed-form formula. This is a "recipe" that tells you exactly how the chaos meter drops at any given time, without needing to run a supercomputer simulation.

3. The Two Stages of Chaos

Their formula reveals two distinct phases of randomness:

• Stage 1: The Early Drop (Power Law)
  At the beginning, the chaos meter drops steadily, like a ball rolling down a hill. The speed of this drop depends on how "squishy" the system is and how strong the random bumps are.
• Stage 2: The Late Plateau (The Limit)
  Eventually, the meter stops dropping and levels off at a specific low value. This is the "maximum randomness" the system can achieve.
  • The Sweet Spot: They found that the system becomes most random (the meter drops the lowest) when the particles are on the verge of becoming a ferromagnet (where they all want to line up in the same direction). It's counter-intuitive: the system is most chaotic right before it tries to organize itself.

4. The "Multiple Quench" Trick

The paper also tested a strategy to make the system even more random. Imagine you are shaking the string.

• Single Shake: You shake it once for a long time.
• Multiple Quenches: Instead of one long shake, you shake it, stop, shake it again with a different random pattern, stop, and repeat.
• The Finding: This "stop-and-start" method works like a turbocharger. The paper shows that doing this multiple times exponentially increases the randomness. It's like shuffling a deck of cards, then shuffling it again with a different technique, then again—the deck becomes perfectly randomized much faster than just shuffling it once for a long time.

5. Checking the Work

To make sure their fancy math wasn't just a theoretical fantasy, they compared their formulas against:

• Exact Diagonalization: Crunching the numbers for small systems where the answer is known to be 100% correct.
• Simulations: Using powerful computer algorithms (TEBD) to simulate larger systems.
• The Verdict: The math matched the computer simulations perfectly across the entire range of conditions they tested.

Summary

In short, this paper provides a perfectly accurate map for how randomness builds up in a specific type of disordered quantum string. They discovered that:

1. You can calculate this randomness exactly using a new formula.
2. The system gets most chaotic near a specific magnetic point.
3. You can supercharge this chaos by shaking the system in multiple short bursts rather than one long burst.

This is a "blueprint" for understanding how quantum systems scramble information, which is crucial for designing better quantum algorithms and simulations.

Technical Summary

1. Problem Statement

The paper addresses a significant gap in the understanding of random unitary dynamics in experimentally accessible, interacting one-dimensional (1D) quantum systems.

• Context: While disorder in 1D interacting systems (Tomonaga-Luttinger Liquids or TLLs) is well-studied via traditional correlation functions (using replica or RG techniques), a complementary perspective from quantum information has emerged. This perspective quantifies the "randomness" of a unitary ensemble using the Frame Potential (FP), a diagnostic for how close a system is to a Haar-random unitary ensemble (essential for quantum algorithms like randomized measurements).
• The Challenge: Analytical treatments of the FP exist for quantum circuits, all-to-all interacting models (e.g., SYK), and stochastic Brownian models. However, no analytical results existed for local, interacting 1D systems with quenched disorder. The primary obstacle is that disorder averaging typically generates non-local replica interactions in field theory, making analytical solutions intractable.
• Goal: To derive a closed-form, non-perturbative expression for the Frame Potential of a disordered TLL and validate it against microscopic lattice models.

2. Methodology

The authors employ a combination of bosonization, Schwinger-Keldysh field theory, and exact diagonalization/TEBD numerics.

• Model Definition:

  • The system is a TLL of length $L$ with Hamiltonian $H_X = H_{\text{TLL}} + H_{\text{dis}, X}$, where $X$ denotes different disorder realizations ($U, V$).
  • Disorder: Quenched Gaussian forward-scattering disorder, coupling linearly to the bosonic field gradient: $H_{\text{dis}} = \int dx \, \xi(x) \partial_x \phi(x)$.
  • Microscopic Realization: The model is mapped to a random-field XXZ spin chain with a specific filtering of the random field to suppress backward scattering (enforcing the forward-scattering condition).

• Theoretical Framework (Keldysh Formalism):

  • The Frame Potential $F^{(k)}(T)$ is defined as the disorder-averaged $2k$-th moment of the trace overlap of unitary evolutions: $F^{(k)}(T) = \mathbb{E}_{U,V} |\text{Tr}(\rho W_U W_V^\dagger)|^{2k}$.
  • The authors map the trace overlap to a path integral on a closed Keldysh contour. For the $k$-th moment, this involves $2k$ Keldysh loops.
  • Key Insight: Because the disorder couples linearly to the bosonic fields, the disorder-averaged action remains exactly quadratic (Gaussian). This allows for an exact analytical solution without perturbation theory.
  • The calculation reduces to evaluating functional determinants of the kernel $\Omega(q; \omega, \omega')$, which combines the clean TLL Green's function and the disorder-induced self-energy.

• Numerical Validation:

  • Exact Diagonalization (ED): Used for the free-fermion point ($\Delta=0$) to access long times.
  • Time-Evolving Block Decimation (TEBD): Used for interacting regimes ($\Delta \neq 0$) with bond dimension $\chi=256$.
  • Filtering: To match the continuum theory, the authors apply higher-order filters to the random field in the lattice model to suppress $2k_F$ (backward) scattering components.

3. Key Contributions and Results

A. Closed-Form Analytical Expression

The authors derive an exact, non-perturbative expression for the normalized Frame Potential ratio $R^{(k)}(T) = F^{(k)}(T)/F^{(k)}(0)$:
$$\ln R^{(k)}(T) = -\frac{1}{2} \sum_q \ln \left[ 1 + k A_q(T) \right] e^{-\alpha|q|}$$
where $A_q(T)$ depends on the dimensionless coupling $g = 8\pi K \gamma / u^2$, the sound velocity $u$, the Luttinger parameter $K$, and the disorder strength $\gamma$.

B. Dynamical Regimes

1. Short-Time Behavior: The FP decays as a power law in time, analogous to standard TLL correlation functions. The decay rate depends on both the Luttinger parameter $K$ and velocity $u$.
2. Late-Time Behavior: The FP does not decay to zero but saturates to a plateau. This plateau value is controlled by a single parameter $g$.
   • The plateau value decreases monotonically with increasing disorder strength $\gamma$, Luttinger parameter $K$, and decreasing sound velocity $u$.
   • Physical Interpretation: Stronger randomness is achieved when the system is highly compressible ($K$ large) and evolves slowly ($u$ small). The optimal regime for randomness is near the Heisenberg ferromagnetic point ($\Delta \to -1$).

C. Multiple-Quench Protocol

The authors generalize the result to a protocol involving $m$ independent quenches (concatenated evolutions with different disorder realizations).

• Result: The logarithm of the Frame Potential scales linearly with the number of quenches $m$: $\ln R^{(k)}_m \approx m \ln R^{(k)}$.
• Implication: This leads to an exponential suppression of the Frame Potential with the number of quenches, significantly enhancing the randomness of the unitary ensemble compared to a single quench.

D. Numerical Verification

• Free-Fermion Benchmark: At $\Delta=0$, the analytical predictions match ED and TEBD results perfectly.
• Interacting Regime: For $\Delta \in [-0.7, 0.5]$, the theory matches TEBD data quantitatively when backward scattering is suppressed (via filtering).
• Deviation: At $\Delta=0.5$ (where $K < 3/2$), backward scattering becomes relevant, causing deviations between the pure forward-scattering theory and numerics. Higher-order filters restore agreement.
• Scaling Collapse: The late-time plateau data for various parameters collapses onto a single universal curve predicted by the analytical formula involving the modified Bessel function $K_0$.

4. Significance and Outlook

• Theoretical Breakthrough: This is the first analytical derivation of the Frame Potential for an experimentally realizable, interacting 1D system with local disorder. It overcomes the "non-local replica interaction" barrier by exploiting the specific quadratic structure of forward-scattering disorder in bosonized theories.
• Algorithm Design: The results provide a concrete guide for optimizing analog quantum simulation platforms (e.g., cold atoms, superconducting circuits) to generate random unitary ensembles (k-designs).
  • Optimization Strategy: Tune system parameters toward the ferromagnetic limit ($K \to \infty, u \to 0$) and utilize multiple-quench protocols to exponentially enhance randomness.
• Experimental Relevance: The study connects abstract quantum information metrics (Frame Potential) to physical parameters (compressibility, sound velocity) in materials like organic conductors, quantum wires, and spin chains.
• Future Directions: The authors suggest extending the theory to finite temperatures, incorporating backward scattering perturbatively to map phase boundaries, and applying these findings to optimize randomized measurement protocols for entropy estimation.

In summary, the paper establishes a rigorous, solvable framework for understanding quantum randomness in disordered 1D matter, bridging the gap between condensed matter physics and quantum information theory.