Exact Quench Dynamics from Thermal Pure Quantum States

This paper presents an exact solution for the entanglement entropy dynamics following a quench from a thermal pure quantum state in a free fermion system, revealing a characteristic double-plateau structure that contrasts with typical linear growth and saturation, as confirmed through 2D conformal field theory on the Klein bottle, matrix Riccati equation numerics, and a quasiparticle picture.

The Gist

Imagine you have a giant, crowded dance floor representing a quantum system. Usually, when physicists study how this dance floor changes over time, they start with a quiet, orderly crowd where everyone is standing still or moving in a simple, predictable line. When you suddenly change the music (a "quench"), the crowd starts to get messy, and the "entanglement" (a measure of how connected the dancers are) grows steadily until it hits a limit. It's like a slow, steady wave of chaos spreading across the room.

This paper is about a very different kind of starting line.

Instead of a quiet crowd, the authors start with a dance floor that is already in a state of maximum chaos and connection. They call this a "Thermal Pure Quantum" (TPQ) state. Imagine a room where every single dancer is already holding hands with a partner on the exact opposite side of the room, even though they are far apart. The room is already "hot" and fully entangled.

The question the authors asked is: What happens if we suddenly change the music in this already chaotic room?

The Big Surprise: The "Double-Plateau" Dance

In normal scenarios, the "messiness" (entanglement) just grows and then stops. But in this specific setup, the authors discovered a bizarre, unique pattern they call a "double-plateau" structure.

Think of it like a rollercoaster ride for the room's connection:

1. The First Flat Spot: For a while, the connection level stays high and steady. Nothing seems to change.
2. The Dip: Suddenly, the connection level drops. It's as if the dancers are briefly letting go of their partners.
3. The Second Flat Spot: The connection level stabilizes again at a new, lower level.

This "up-down-up" (or rather, "flat-down-flat") pattern is completely new and unexpected. It's like the room takes a deep breath, exhales some of its tension, and then settles into a new rhythm.

How Did They Solve This?

Solving this mathematically is usually impossible because the system is so complex. It's like trying to predict the path of every single grain of sand in a hurricane. However, the authors used three different "superpowers" to crack the code, and they all agreed on the answer:

1. The Magic Mirror (Conformal Field Theory):
   They used a branch of math called Conformal Field Theory (CFT), but with a twist. Instead of looking at the system on a normal cylinder (like a soda can), they looked at it on a Klein Bottle.

   • Analogy: Imagine a soda can where the top and bottom are connected, but the left and right sides are twisted and connected in a way that creates a loop with no inside or outside. By doing the math on this weird, twisted shape, they could predict the exact "dance moves" of the entanglement.

2. The Super-Computer (Numerical Simulation):
   They built a massive digital model of the system (a chain of spins, like a row of tiny magnets) and ran it on a computer. Instead of guessing, they solved the equations exactly.

   • Analogy: It's like simulating a billion billiard balls hitting each other on a computer to see exactly how they scatter, confirming that the "double-plateau" pattern is real and not just a math trick.

3. The Messenger Pairs (Quasiparticle Picture):
   They explained why this happens using a simple story about messengers.

   • The Setup: In their starting state, every dancer has a partner standing on the exact opposite side of the room.
   • The Action: When the music changes, these partners start walking toward each other.
   • The Dip: As the partners walk into the specific section of the room we are watching, they "cancel out" the messiness. It's like two waves meeting and smoothing each other out. This causes the entanglement to drop.
   • The Second Plateau: Once the partners have walked past each other and left the section, the messiness settles at a new, lower level.

Why Does This Matter?

This isn't just about abstract math. It helps us understand how information is stored and scrambled in complex systems.

• Black Holes: The "crosscap" state they used is mathematically similar to the inside of a black hole. Understanding how these states evolve might help us understand how black holes process information.
• Quantum Computers: As we build quantum computers, we need to know how these highly connected states behave when things go wrong or change. This paper gives us a precise map of that behavior.

The Takeaway

The authors took a system that was already "maxed out" on chaos and showed that when you poke it, it doesn't just get messier. Instead, it performs a specific, predictable dance: it holds steady, lets go of some tension, and settles into a new state. They proved this using three different methods that all told the same story, revealing a hidden rhythm in the quantum world that we never knew existed.

Technical Summary

1. Problem Statement

The paper addresses the non-equilibrium dynamics of isolated quantum systems, specifically focusing on entanglement entropy evolution following a quantum quench from a Thermal Pure Quantum (TPQ) state.

• Context: Standard quench dynamics typically start from low-entanglement states (e.g., ground states or product states), leading to a characteristic linear growth of entanglement entropy followed by saturation.
• The Challenge: TPQ states are high-energy, pure states that are locally indistinguishable from thermal Gibbs ensembles but possess volume-law entanglement (maximal entanglement) from the outset. Simulating the dynamics of such highly entangled states is computationally prohibitive (exponential complexity) for standard numerical methods.
• Specific Goal: The authors investigate the real-time evolution of bipartite entanglement entropy in a free fermion system (mapped to an XX spin chain) starting from a structured TPQ state $|\Psi_\beta\rangle$, defined by imaginary-time evolution of a crosscap state $|C\rangle$.

2. Methodology

The authors employ three complementary, exact approaches to solve the problem, ensuring rigorous validation of their results:

A. 2D Conformal Field Theory (CFT) Approach

• Framework: The system is modeled as a $c=1$ free compact boson CFT (dual to a massless Dirac fermion).
• Geometry: The initial state preparation (crosscap boundary conditions) maps the problem to a Klein bottle geometry rather than the standard cylinder or torus.
• Calculation:
  • The TPQ state is defined as $|\Psi_\beta\rangle \propto e^{-\beta H/4}|C\rangle$.
  • Using the replica trick, the entanglement entropy $S_A(t)$ is derived from the $n$-th moment of the reduced density matrix, $\text{Tr}[\rho_A^n]$.
  • This is computed as a path integral involving vertex operator correlators on the Klein bottle.
  • The result is expressed analytically using Dedekind eta functions ($\eta$) and Jacobi theta functions ($\theta_{1,2}$).

B. Exact Numerical Simulation (Gaussian States)

• Model: The system is realized as a spin-1/2 XX chain at half-filling, which maps to non-interacting fermions via the Jordan-Wigner transformation.
• Technique: Since the initial state and the Hamiltonian are quadratic, the system remains a Gaussian state throughout the evolution.
• Implementation:
  • The state is fully characterized by a fermionic covariance matrix $\Gamma$.
  • Imaginary-time evolution (state preparation) is solved exactly using the Matrix Riccati equation: $\frac{d\Gamma}{d\tau} = -H - \Gamma H \Gamma$.
  • Real-time evolution (quench) is solved via unitary transformation: $\Gamma(t) = e^{Ht}\Gamma_\beta e^{-Ht}$.
  • Entanglement entropy is calculated directly from the eigenvalues of the restricted covariance matrix.

C. Asymptotically Exact Quasiparticle Picture

• Concept: The authors adapt the standard quasiparticle picture (usually for low-entanglement quenches) to the TPQ scenario.
• Mechanism: Instead of creating entanglement, the dynamics involve the transport and destruction of pre-existing non-local entanglement.
• Model: The initial state is viewed as a sea of antipodally entangled quasiparticle pairs $(k, -k)$ located at sites $x$ and $x+L/2$.
• Dynamics: As pairs propagate ballistically, the entanglement entropy decreases when both members of a pair enter the subsystem $A$, and increases when they leave.

3. Key Results

The Double-Plateau Structure

Contrary to the standard linear growth and saturation observed in ground-state quenches, the entanglement entropy $S_A(t)$ exhibits a unique double-plateau structure:

1. Initial Plateau: The entropy remains constant at its initial volume-law value ($S_A(0) \propto l$) for a short duration.
2. Decrease: The entropy drops as antipodal pairs enter the subsystem, effectively "removing" entanglement from the boundary.
3. Second Plateau: The entropy rises again and saturates at a lower value as the pairs traverse the subsystem and exit.

Analytical Formula

The CFT approach yields an exact analytical formula for the time-dependent entanglement entropy:
$$S_A(t, \sigma) = \frac{1}{6} \log \left| \frac{\eta(i\beta/2\pi)^{-6} \theta_1(\frac{\sigma}{2\pi} | \frac{i\beta}{2\pi}) \theta_2(\frac{\beta+4it}{4\pi i} | \frac{i\beta}{2\pi})}{\theta_2(\frac{\beta+4it+2i\sigma}{4\pi i} | \frac{i\beta}{2\pi}) \theta_2(\frac{\beta+4it-2i\sigma}{4\pi i} | \frac{i\beta}{2\pi})} \right|^2$$
where $\sigma$ is the subsystem length, $\beta$ is the inverse temperature parameter, and $t$ is time.

Validation

• Numerical Agreement: The exact numerical simulation of the XX chain matches the CFT prediction perfectly (after appropriate finite-size scaling).
• Quasiparticle Agreement: The quasiparticle picture, utilizing the specific entanglement density $s(k)$ derived from the BCS-like structure of the TPQ state, quantitatively reproduces the double-plateau profile.

4. Significance and Contributions

1. Exact Solution for High-Entanglement Dynamics: The paper provides one of the rare exact solutions for the real-time dynamics of a system starting from a volume-law entangled state, overcoming the usual computational barriers.
2. New Dynamical Phenomenology: It identifies a distinct "double-plateau" behavior in entanglement entropy, fundamentally different from the linear growth seen in standard quenches. This highlights that the transport of pre-existing non-local correlations dominates the dynamics.
3. Geometric Insight (Klein Bottle): It successfully extends the CFT quench framework to non-orientable manifolds (Klein bottle), demonstrating how crosscap boundary conditions (representing antipodal entanglement) alter the topology of the spacetime used for calculation.
4. Bridge Between Formalisms: The work rigorously connects three distinct frameworks—CFT on non-trivial geometries, exact Gaussian numerics (Riccati equations), and quasiparticle transport—showing their consistency in describing complex non-equilibrium phenomena.
5. Implications for Holography and Black Holes: Since TPQ states constructed from crosscaps are used to model microstates of single-sided black holes with smooth interiors, these results offer insights into how quantum information is scrambled and redistributed in such holographic systems.

5. Conclusion

The authors demonstrate that quenching from a Thermal Pure Quantum state leads to a non-monotonic entanglement evolution characterized by a double-plateau structure. This behavior arises from the ballistic transport of pre-existing antipodal entanglement. The study establishes a robust theoretical and numerical framework for exploring non-equilibrium dynamics in highly entangled regimes, with potential applications in quantum simulation and holographic duality.