It is an open question how well tensor network states in the form of an
infinite projected entangled pair states (iPEPS) tensor network can approximate
gapless quantum states of matter. Here we address this issue for two different
physical scenarios: i) a conformally invariant (2+1)d quantum critical point
in the incarnation of the transverse field Ising model on the square lattice
and ii) spontaneously broken continuous symmetries with gapless Goldstone modes
exemplified by the S=1/2 antiferromagnetic Heisenberg and XY models on the
square lattice. We find that the energetically best wave functions display {\em
finite} correlation lengths and we introduce a powerful finite correlation
length scaling framework for the analysis of such finite-D iPEPS states. The
framework is important i) to understand the mild limitations of the finite-D
iPEPS manifold in representing Lorentz-invariant, gapless many body quantum
states and ii) to put forward a practical scheme in which the finite
correlation length ξ(D) combined with field theory inspired formulae can be
used to extrapolate the data to infinite correlation length, i.e. to the
thermodynamic limit. The finite correlation length scaling framework opens the
way for further exploration of quantum matter with an (expected)
Lorentz-invariant, massless low-energy description, with many applications
ranging from condensed matter to high-energy physics.Comment: 16 pages, 11 figure