Expectation values after an integrable boundary quantum quench

This paper develops a general framework based on form factors to analyze the real-time dynamics of integrable boundary quantum quenches, specifically investigating the time evolution of the pre-quench vacuum in the Lee-Yang model and validating these analytical results with numerical calculations using the truncated conformal space approach.

The Gist

Imagine a long, narrow hallway (a "strip") where the rules of physics are perfectly predictable and orderly. This is the world of the Lee–Yang model, a specific type of quantum system the authors are studying.

At one end of this hallway, the walls are painted a specific color (let's call it "Boundary A"). At the other end, the wall is also "Boundary A." The system is sitting quietly in its most relaxed state, like a calm lake.

The "Quench": A Sudden Paint Job
Suddenly, at time zero, someone rushes to the right end of the hallway and instantly repaints the wall from "Boundary A" to "Boundary B." In physics terms, this is called a boundary quantum quench. It's not a slow change; it's an abrupt switch.

The paper asks a simple but deep question: What happens next?

When you change the wall color, it doesn't just stay changed at that one spot. The "news" of the new color ripples through the hallway. The authors want to track exactly how this ripple moves, how it changes the energy of the system, and how the "calm lake" settles into a new state.

The Two Tools Used to Solve the Puzzle

To figure this out, the authors used two very different but complementary methods, like using both a telescope and a microscope to study a star.

1. The "Perfect Map" (Form Factors)
First, they used a mathematical technique called Form Factors. Think of this as having a perfect, pre-drawn map of how particles behave in this specific hallway.

• Because the system is "integrable" (meaning it follows strict, solvable rules), the authors could calculate exactly how the "ripple" of the new boundary condition travels.
• They found that the ripple moves at the speed of light (in this quantum world).
• They discovered a fascinating "echo" effect. When the ripple hits the opposite wall (the left side), it bounces back. It keeps bouncing back and forth between the two walls, creating a rhythmic pattern.
• The Surprise: Usually, when a ripple hits a wall, it might just fade away or bounce back loudly. But here, the authors found that the "direct" ripple and the "reflected" ripple cancel each other out in a very specific way. Instead of fading slowly, the system settles down in a way that follows a specific mathematical rhythm (oscillating and slowing down like $1/t$). It's like two waves crashing together and creating a perfectly calm spot for a moment before the next wave arrives.

2. The "Digital Simulator" (TCSA)
To make sure their "Perfect Map" wasn't just a pretty theory, they built a Digital Simulator (called Truncated Conformal Space Approach, or TCSA).

• Imagine trying to simulate a storm on a computer. You can't calculate every single drop of water, so you only calculate the biggest, most important drops. This is what "truncation" means: simplifying the math by ignoring the tiniest details to make the computer run.
• The authors ran their simulation to see if the digital "ripple" matched the "Perfect Map."
• The Problem: At first, the simulation looked messy. It had "static" or "noise" (oscillations) that the perfect map didn't predict.
• The Fix: The authors realized this noise wasn't a mistake in the physics; it was an artifact of the simulation's limits (ignoring the tiny drops). They developed a clever "noise-canceling" technique. By mathematically subtracting the known errors of the simulation, they cleaned up the data.
• The Result: Once the noise was removed, the simulation matched the "Perfect Map" perfectly. The digital ripple behaved exactly as the theory predicted.

The Big Picture
The paper is essentially a success story of cross-verification.

• Theory said: "If you change the wall, the ripple will bounce back and forth, and the system will settle down in this specific, rhythmic way."
• Simulation said: "We tried to build it, and it looked messy at first, but once we fixed our tools, it matched the theory exactly."

Why Does This Matter?
The authors used this specific "Lee–Yang" hallway as a test case. It's a simple, non-unitary (a bit weird mathematically) model, but it's the perfect training ground. By proving that their "Form Factor" map and their "Digital Simulator" agree on this simple model, they have built a reliable toolkit.

They are essentially saying: "We have a new, reliable way to predict what happens when you suddenly change the rules at the edge of a quantum system. We tested it, and it works." This gives them confidence to apply these same tools to more complex, real-world quantum systems in the future.

Technical Summary

Technical Summary: Expectation values after an integrable boundary quantum quench

Problem Statement
This work investigates a specific class of non-equilibrium dynamics in quantum field theory: an integrable boundary quantum quench. Unlike global quenches where a Hamiltonian parameter is changed uniformly, or local quenches involving bond modifications, this protocol involves the instantaneous switching of a boundary condition. Specifically, the system is initialized in the ground state of a Hamiltonian with fixed boundary conditions ($\gamma$ and $\alpha$). At time $t=0$, a boundary-changing operator $\psi_{\beta\alpha}$ is inserted at one boundary, switching the condition from $\alpha$ to $\beta$. The authors aim to compute the real-time evolution of vacuum-to-vacuum matrix elements of local bulk operators $\Phi(x,t)$ inserted after the quench, defined as:
$$\gamma,\beta \langle 0 | \Phi(x, t) \psi_{\beta\alpha}(0, 0) | 0 \rangle_{\gamma,\alpha}$$
The study focuses on the Lee–Yang model, a non-unitary minimal model ($c = -22/5$), first in its conformal field theory (CFT) limit and subsequently in its massive integrable deformation.

Methodology
The paper employs a dual approach, combining analytical form-factor expansions with numerical verification via the Truncated Conformal Space Approach (TCSA).

1. Conformal Field Theory (CFT) Analysis:

   • The problem is mapped from a strip geometry to the upper half-plane (UHP) using an exponential map.
   • The correlation functions are reduced to chiral four-point functions involving the bulk field and boundary-changing operators.
   • Using the degenerate nature of the Lee–Yang primary field $\phi$ (level 2), the authors solve the resulting second-order differential equations to determine the space-time dependence of the correlators.
   • Two specific scenarios are analyzed: switching from identical boundaries ($\alpha=\gamma$) to different ones, and switching from different boundaries to identical ones.

2. Massive Integrable Theory Analysis:

   • The authors develop a general framework based on a "double form-factor expansion."
   • Boundary State Formalism: The pre-quench state is treated via a boundary state $|B\rangle_\alpha$ constructed from the reflection factors of the massive theory.
   • Double Expansion: The time evolution is expanded by inserting a complete set of post-quench multiparticle states (boundary channel) and subsequently resolving the identity in the bulk channel. This leads to an expression involving:
     • Boundary-changing form factors $F^{\psi_{\beta\alpha}}_n$.
     • Bulk form factors $F^\Phi_n$.
     • Boundary reflection factors $R_\beta(\theta)$.
     • Boundary-state amplitudes.
   • Asymptotic Analysis: The large-time behavior is analyzed using the stationary phase approximation. A key finding is that for integrable boundary conditions, the reflection factor satisfies $R_\beta(0) = -1$, leading to a destructive interference between direct and reflected particle contributions. This cancels the leading $t^{-1/2}$ decay term, resulting in a $t^{-1}$ algebraic decay modulated by oscillations at the particle mass frequency ($e^{imt}$).

3. Numerical Verification (TCSA):

   • The authors adapt the TCSA to boundary-changing situations. The Hamiltonian is truncated at a specific energy level $\Lambda$ in a finite volume $L$.
   • To address the discrepancy between finite-volume TCSA data and continuum form-factor predictions (specifically regarding light-cone singularities), the authors implement a renormalization procedure. This involves subtracting the known CFT cutoff dependence (derived analytically) from the massive TCSA results and adding back the exact CFT result, effectively isolating the massive corrections.
   • The method validates the form-factor expansion by comparing the renormalized TCSA time evolution of the bulk one-point function against the analytical one-particle approximation.

Key Contributions and Results

• General Framework: The paper establishes a systematic method for analyzing boundary quenches using a double form-factor expansion, linking boundary-changing operators to bulk dynamics via boundary states and reflection factors.
• Analytical Results in Lee–Yang Model:
  • Exact CFT correlators are derived for both $I \to \phi$ and $\phi \to I$ boundary switches.
  • In the massive theory, the leading asymptotic behavior of the bulk expectation value is derived. The authors demonstrate that the boundary reflection causes a cancellation of the leading stationary phase term, changing the relaxation decay from $t^{-1/2}$ to $t^{-1}$.
  • The light-cone structure is identified as a purely Minkowski effect, characterized by the propagation of the boundary change at the speed of light and subsequent reflections between boundaries.
• Numerical Validation:
  • The TCSA results, after renormalization, show remarkable agreement with the analytical form-factor predictions for the time evolution of the bulk field expectation value.
  • The study clarifies that the "smearing" of the light-cone peak in raw TCSA data is an artifact of the ultraviolet cutoff (truncation of high-energy rapidity modes) rather than a failure of the form-factor approach. By imposing a matching rapidity cutoff on the analytical integral, the TCSA signal is perfectly reproduced.
  • The matrix elements of bulk fields between boundary multiparticle states are computed and verified against TCSA data for one-particle states.

Significance and Claims
The authors claim that their work provides a viable computational approach to boundary-changing quenches in integrable systems. The primary significance lies in the successful cross-verification of two distinct methods:

1. Exact Solvability: The form-factor bootstrap approach, which relies on the integrability of the model.
2. Hamiltonian Truncation: The TCSA, which is applicable to generic quantum field theories but suffers from cutoff artifacts.

The paper argues that the comparison serves as a non-trivial benchmark for integrable calculations and offers a guide for improving renormalized truncation methods in non-equilibrium settings. The authors explicitly state that the Lee–Yang results demonstrate the viability of their framework, while acknowledging that extending this to unitary models (like Ising or Potts) or spin-chain realizations (via tensor networks) remains a future direction. They do not claim immediate experimental realization but suggest that such protocols could be relevant for quantum simulation platforms.