  We describe an iterative method to optimize the multi-scale entanglementrenormalization ansatz (MERA) for the low-energy subspace of local Hamiltonianson a D-dimensional lattice. For translation invariant systems the cost of thisoptimization is logarithmic in the linear system size. Specialized algorithmsfor the treatment of infinite systems are also described. Benchmark simulationresults are presented for a variety of 1D systems, namely Ising, Potts, XX andHeisenberg models. The potential to compute expected values of localobservables, energy gaps and correlators is investigated.
We describe an iterative method to optimize the multi-scale entanglement renormalization ansatz (MERA) for the low-energy subspace of local Hamiltonians on a D-dimensional lattice. For translation invariant systems the cost of this optimization is logarithmic in the linear system size. Specialized algorithms for the treatment of infinite systems are also described.