Conformal data and renormalization group flow in critical quantum spin chains using periodic uniform matrix product states

TL;DR

It is established that a Bloch-state ansatz based on periodic uniform matrix product states (PUMPS), originally designed to capture single-quasiparticle excitations in gapped systems, is in fact capable of accurately approximating all low-energy eigenstates of critical quantum spin chains on the circle.

Abstract

We establish that a Bloch-state ansatz based on periodic uniform Matrix Product States (puMPS), originally designed to capture single-quasiparticle excitations in gapped systems, is in fact capable of accurately approximating all low-energy eigenstates of critical quantum spin chains on the circle. When combined with the methods of [Milsted, Vidal, Phys. Rev. B 96 245105] based on the Koo-Saleur formula, puMPS Bloch states can then be used to identify each low-energy eigenstate of a chain made of up to hundreds of spins with its corresponding scaling operator in the emergent conformal field theory (CFT). This enables the following two tasks, that we demonstrate using the quantum Ising model and a recently proposed generalization thereof due to O'Brien and Fendley [Phys. Rev. Lett. 120, 206403]. (i) From the spectrum of low energies and momenta we extract conformal data (specifying the emergent CFT) with unprecedented numerical accuracy. (ii) By changing the lattice size, we investigate nonperturbatively the RG flow of the low-energy spectrum between two CFTs. In our example, where the flow is from the Tri-Critical Ising CFT to the Ising CFT, we obtain excellent agreement with an analytical result [Klassen and Melzer, Nucl. Phys. B 370 511] conjectured to describe the flow of the first spectral gap directly in the continuum.

Figures (15)

• Figure 1: Top: Scaling operator spectra of (a) the Ising and (b) the TCI CFTs (with a selection of operators labeled). Bottom: Approximate scaling dimensions and conformal-tower identification for the OF model at $\lambda=0.4$ with (c) $N=256, D=52$ and (d) $N=32, D=32$, corresponding to points of Fig. \ref{['fig:flow_N']}. We label a selection of states according to a numerical identification of the corresponding CFT operators milsted_extraction_2017supplemental. Note: We displace data points slightly along the x-axis to show degeneracies.
• Figure 2: Spectral RG flow (crosses) of the first 5 energy levels (as apparent scaling dimensions $\Delta$) at momentum zero, excluding $\Delta=0$, extracted from the OF model with $\lambda=0.4$, using puMPS with $D \le 52$. For comparison, we also plot the exact scaling dimensions of the Ising and TCI CFTs (dots, diamonds). The crossover between the two highest levels plotted, which we confirm by tracking conformal tower membership using $H_n$ matrix elements, is consistent with these states belonging to different Kramers-Wannier self-duality sectors.
• Figure 3: Connection of the spectral RG flow of Fig. \ref{['fig:flow_N']} (left) to the "flow" of OF model energy levels as a function of $\lambda$ at fixed system size $N=32$, computed using puMPS with $D=28$. Note how the apparent scaling dimensions agree with the TCI CFT values at the TCI point $\lambda_\textsl{\tiny TCI} \approx 0.428$.
• Figure 4: Flow of the first spectral gap from Fig. \ref{['fig:flow_N']} compared supplemental with the integrable field theory result of klassen_spectral_1992, conjectured to describe the equivalent flow in the continuum.
• Figure 5: The tensor networks for the derivative of the auxiliary energy function with respect to the central tensor in \ref{['energy_derivative']} for an infinite translation invariant MPS. Top: the denominator, i.e. the square of the norm of locally deformed puMPS with central tensor $A_C$. It equals the vector norm of the tensor $A_C$. Bottom: the numerator, where red tensors form a matrix product operator representation of the shifted Hamiltonian $\tilde{H}$.