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.

It is shown that, for a Hamiltonian with a band structure, the half width of
local spectral density of states, or strength function, is closely related to
the width of the nonperturbative (NPT) parts of energy eigenfunctions. In the
Wigner-band random-matrix model, making use of a generalized Brillouin-Wigner
perturbation theory, we derive analytical expressions for the width of the NPT
parts under weak and strong perturbation. An iterative algorithm is given, by
which the NPT widths can be computed efficiently, and is used in numerical test
of the analytical predictions.