Quench dynamics of the quantum XXZ chain with staggered interactions: Exact results and simulations on digital quantum computers

This paper investigates the exact quench dynamics of the quantum $S=1/2$ XXZ antiferromagnetic chain with staggered interactions in the flat-band limit, deriving closed-form expressions for entanglement entropies and Loschmidt echoes while validating these theoretical predictions through numerical simulations on IBM-Q quantum devices.

The Gist

Imagine a long line of dancers holding hands, representing a chain of tiny quantum magnets (spins). In this specific dance, the partners are arranged in pairs (dimers), and they are holding hands very tightly in an alternating pattern: tight, loose, tight, loose.

This paper is about what happens when you suddenly swap the strength of their grip. You take the "tight" pairs and make them "loose," and the "loose" pairs become "tight." This sudden change is called a quantum quench.

The researchers wanted to know: How does this line of dancers react over time? Do they eventually calm down and settle into a new rhythm, or do they keep dancing wildly forever?

Here is a breakdown of their findings using simple analogies:

1. The "Flat-Band" Dance Floor

Usually, when you push a line of dancers, a wave of movement travels down the line, and eventually, everyone settles down. But in this specific setup (called the "flat-band limit"), the dancers are stuck in their local pairs. They can't easily pass their energy to their neighbors.

• The Result: The system never relaxes. Instead of calming down, the dancers keep oscillating back and forth in a perfect, repeating pattern forever. It's like a pendulum that never stops swinging because there is no friction.

2. The "Bell Basis" Secret Code

To predict exactly how the dancers move, the researchers used a special mathematical language called the Bell basis.

• The Analogy: Imagine trying to describe a complex dance routine by breaking it down into four basic "moves" (like a spin, a jump, a slide, or a pose). The researchers found that no matter how long the line of dancers is, the entire complex dance can be described as a mix of just these four basic moves.
• The Discovery: They figured out a precise rulebook (selection rules) that tells them exactly which combinations of these moves are allowed and which are forbidden. This allowed them to write down the exact solution for the dance at any point in time, without needing a supercomputer.

3. Measuring "Entanglement" (The Invisible String)

In quantum physics, "entanglement" is like an invisible string connecting two dancers. If one moves, the other knows instantly, even if they are far apart.

• The Finding: The researchers measured how strong these invisible strings get over time. They found that the "string tension" (entanglement entropy) goes up and down in a perfect, repeating cycle. It never settles at a steady level.
• The Surprise: They discovered that the length of the chain changes the speed of the cycle. If the chain has an even number of dancers, the rhythm is fast. If it has an odd number, the rhythm is exactly half as fast. It's like a drumbeat that changes tempo depending on whether the band has an even or odd number of members.

4. The "Return Probability" (Did we come home?)

The researchers also tracked something called the Loschmidt Echo.

• The Analogy: Imagine the dancers start in a specific formation. As time passes, they dance around. The Loschmidt Echo asks: "At this exact moment, how much do the dancers look like they did at the very start?"
• The Result: Sometimes, the dancers return to their exact starting pose (the echo is 100%). Other times, they look completely different (the echo is 0%).
• The "Zeros": They found specific moments in time where the dancers never look like the start, no matter how long the chain is (for certain lengths). These are called "Loschmidt zeros." It's like a clock that hits a specific time where the hands stop moving entirely.

5. Testing on a "Quantum Robot" (IBM Quantum Computer)

Since they couldn't build a real quantum chain of atoms in a lab easily, they used a digital quantum computer (IBM's quantum processor) to simulate the dance.

• Method A (The Hadamard Test): They used a clever trick to "peek" at the mathematical coefficients of the dance moves. This worked well for small groups (4 to 6 dancers) but got too noisy and complicated for larger groups.
• Method B (Randomized Measurements): For larger groups (up to 12 dancers), they used a different strategy. They let the quantum computer run the dance, then took "snapshots" of the dancers in random poses. By taking thousands of these random snapshots and using a statistical trick called "Classical Shadows," they could reconstruct the whole dance pattern.
• The Outcome: The digital simulation matched their mathematical predictions almost perfectly, proving that their theory is correct and that current quantum computers can simulate these complex quantum dances.

Summary

This paper is a success story of combining pure math with modern technology.

1. Math: They solved a complex quantum puzzle exactly, showing that in this specific "flat-band" system, the system never forgets its past and keeps dancing in a loop forever.
2. Tech: They proved that today's noisy quantum computers are powerful enough to simulate these complex dynamics, provided you use the right measurement tricks (like taking random snapshots) to filter out the noise.

It's a bit like proving you can predict the exact path of a bouncing ball in a frictionless room, and then successfully filming that prediction on a shaky, low-quality camera.

Technical Summary

1. Problem Statement

The paper investigates the non-equilibrium dynamics of a quantum many-body system following a quantum quench. Specifically, it focuses on the $S=1/2$ XXZ antiferromagnetic chain in the flat-band limit with staggered (dimerized) and anisotropic interactions.

• The System: A chain of $N=2n$ spins with alternating bond strengths ($J$ and $J'$) and an anisotropy parameter $\Delta$.
• The Protocol: A global quench is performed by initializing the system in the ground state of a fully dimerized Hamiltonian (strong bonds on odd links, $J>0, J'=0$) and suddenly interchanging the bond strengths (strong bonds on even links, $J=0, J'>0$).
• The Challenge: While integrable systems typically relax to a Generalized Gibbs Ensemble, this specific system exhibits confinement due to vanishing group velocity of excitations (flat band). Consequently, the system fails to thermalize, leading to persistent oscillations in entanglement and observables. The authors aim to derive exact analytical solutions for these dynamics and validate them using noisy intermediate-scale quantum (NISQ) devices.

2. Methodology

The study employs a dual approach: exact analytical derivation and digital quantum simulation.

A. Analytical Approach: The Bell Basis

The authors exploit the fact that both the pre-quench and post-quench Hamiltonians can be written as sums of non-interacting dimers.

• Basis Transformation: They work in the Bell basis (valence bond basis), where each dimer's eigenstates are the four Bell states ($|\psi_0\rangle$ to $|\psi_3\rangle$).
• Expansion Coefficients: The initial state (a product of singlets on odd links) is expanded in terms of the post-quench eigenstates (products of Bell states on even links).
  • They derive selection rules for non-zero expansion coefficients ($a_k$):
    1. The sum of dimer-type indices must be even.
    2. The combined count of specific triplet dimers ($\psi_2, \psi_3$) must be even.
  • All non-zero coefficients have equal magnitude: $|a_k| = 2^{-(n-1)}$.
• Exact Solutions: Using these coefficients and the known energy eigenvalues, they derive closed-form, size-independent expressions for:
  • Time-dependent states $|\Psi_0(t)\rangle$.
  • Entanglement Entropies: Both von Neumann ($S_1$) and second-order Rényi ($S_2$).
  • Loschmidt Echo ($L(t)$): The return probability $|\langle \Psi_0 | \Psi_0(t) \rangle|^2$ and the return rate function $r(t) = -\frac{1}{N} \ln L(t)$.

B. Digital Quantum Simulation

The authors validate their theory using IBM-Q quantum processors (specifically the ibm_fez device) via two distinct experimental protocols:

1. Hadamard Test for Amplitude Estimation:

   • Used to estimate the expansion coefficients $a_k = \langle \Phi_k | \Psi_0 \rangle$ directly.
   • Circuits prepare the initial state and the Bell basis states, then use an ancilla qubit to measure the real part of the overlap.
   • Limitation: Effective only for small systems ($N=4$) due to exponential growth in circuit depth and noise for larger $N$.

2. Trotter-Free Time Evolution + Randomized Measurements:

   • Time Evolution: Constructs the time-evolved state using a circuit that is Trotter-error-free for this specific model. The evolution operator is decomposed into exact exponentials of Pauli terms ($XX, YY, ZZ$) without approximation errors, regardless of time step size.
   • Measurement: Employs randomized Pauli measurements (measuring qubits in random $X, Y, Z$ bases).
   • Post-Processing: Uses two classical shadow techniques to extract observables:
     • Hamming Distance Method: Estimates purity via statistical correlations of bitstrings.
     • Classical Shadows: Reconstructs the density matrix to estimate the second-order Rényi entropy and Loschmidt echo.

3. Key Contributions and Results

A. Exact Analytical Results

• Entanglement Scaling:
  • For odd $n$, the entanglement entropy follows $S_{n=odd}(t) = 1 + S_{n=2}(t/2)$.
  • For even $n > 2$, it follows $S_{n=even}(t) = 2 S_{n=2}(t/2)$.
  • This reveals that the oscillation period doubles for odd $n > 2$ compared to the $n=2$ case due to boundary effects.
• Periodicity: The dynamics are periodic with period $2q\pi$ if $J\Delta = p/q$ (rational). If $\Delta$ is irrational, the dynamics are non-periodic.
• Loschmidt Zeros and DPTs:
  • Identified exact zeros of the Loschmidt echo in finite chains for specific anisotropies (e.g., $\Delta=0, 1/2, 7/4$) at times $t = \pi, 3\pi, \dots$.
  • These zeros correspond to Dynamical Quantum Phase Transitions (DQPTs) where the return rate function $r(t)$ diverges.
  • Finite-Size Scaling:
    • For odd $n$, the Loschmidt echo at critical times $t^* = (2m+1)\pi$ scales universally as $L^* = (1/2)^{2(n-1)}$, independent of $\Delta$.
    • For even $n$, the behavior depends on $n \pmod 4$:
      • $n = 4\nu + 2$: Exact zeros occur ($L^*=0$).
      • $n = 4\nu + 4$: $L^* = (1/2)^{2(n-2)}$.

B. Simulation Results

• Hadamard Test: Successfully reconstructed the state and calculated entropies for $N=4$ on real hardware, matching exact results. Accuracy degraded for $N=6$ due to circuit depth and noise.
• Randomized Measurements:
  • Successfully simulated systems up to $N=12$ qubits.
  • Results for both Rényi entropy and Loschmidt echo showed satisfactory agreement with exact analytical solutions, even without advanced error mitigation.
  • Demonstrated that the classical shadow technique and Hamming distance method yield consistent results when applied to the same measurement data.

4. Significance

• Theoretical Insight: The paper provides a rare example of a non-trivial interacting quantum system where exact, closed-form solutions for dynamical observables (entanglement and Loschmidt echo) can be derived for arbitrary system sizes. It clarifies the role of flat bands and confinement in preventing thermalization.
• DQPT Characterization: It establishes universal finite-size scaling laws for Loschmidt zeros in the flat-band limit, showing that DQPTs occur at critical times independent of the anisotropy parameter for odd system sizes.
• Quantum Hardware Benchmarking: The work serves as a benchmark for NISQ devices.
  • It demonstrates that Trotter-free circuits are crucial for simulating specific models on noisy hardware.
  • It validates randomized measurement protocols (classical shadows) as a scalable and robust method for extracting entanglement and echo properties on real quantum processors, overcoming the limitations of direct state tomography.
• Future Directions: The study opens avenues for exploring finite-rate quenches, shallow-circuit variants of the Hadamard test for intermediate system sizes, and further investigation into the system-size dependence of Loschmidt zeros.

In summary, this paper bridges rigorous analytical many-body physics with practical quantum simulation, offering a comprehensive understanding of quench dynamics in a flat-band XXZ chain and demonstrating the capability of current quantum hardware to probe complex non-equilibrium phenomena.