The first sentence of the abstract should read as "We describe an iterative method to optimize the multiscale entanglement renormalization ansatz for the...".
We show how to build a multiscale entanglement renormalization ansatz (MERA) representation of the ground state of a many-body Hamiltonian H by applying the recently proposed tensor network renormalization [G. Evenbly and G. Vidal, Phys. Rev. Lett. 115, 180405 (2015)] to the Euclidean time evolution operator e(-βH) for infinite β. This approach bypasses the costly energy minimization of previous MERA algorithms and, when applied to finite inverse temperature β, produces a MERA representation of a thermal Gibbs state. Our construction endows tensor network renormalization with a renormalization group flow in the space of wave functions and Hamiltonians (and not merely in the more abstract space of tensors) and extends the MERA formalism to classical statistical systems.
The continuous multiscale entanglement renormalization ansatz (cMERA) [J. Haegeman et al., Phys. Rev. Lett. 110, 100402 (2013)] gives a variational wave functional for ground states of quantum field-theoretic Hamiltonians. A cMERA is defined as the result of applying to a reference unentangled state a unitary evolution generated by a quasilocal operator, the entangler. This makes the extension of the formalism to the case where boundaries and defects are present nontrivial. Here we show how this generalization works, using the (1+1)-dimensional free boson cMERA as a proof-of-principle example, and restricting ourselves to conformal boundaries and defects. In our prescription, the presence of a boundary or defect induces a modification of the entangler localized only to its vicinity, in analogy with the so-called principle of minimal updates for the lattice tensor network MERA.