General rightsIt is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), other than for strictly personal, individual use, unless the work is under an open content license (like Creative Commons). Disclaimer/Complaints regulationsIf you believe that digital publication of certain material infringes any of your rights or (privacy) interests, please let the Library know, stating your reasons. In case of a legitimate complaint, the Library will make the material inaccessible and/or remove it from the website. Please Ask the Library: http://uba.uva.nl/en/contact, or a letter to: Library of the University of Amsterdam, Secretariat, Singel 425, 1012 WP Amsterdam, The Netherlands. You will be contacted as soon as possible. Variational studies of the t-J model on the square lattice based on infinite projected-entangled pair states confirm an extremely close competition between a uniform d-wave superconducting state and different stripe states. The site-centered stripe with an in-phase d-wave order has an equal or only slightly lower energy than the stripe with antiphase d-wave order. The optimal stripe filling is not constant but increases with J=t. A nematic anisotropy reduces the pairing amplitude and the energies of stripe phases are lowered relative to the uniform state with increasing nematicity.

† These authors contributed equally to the calculations in this work.Competing inhomogeneous orders are a central feature of correlated electron materials including the high-temperature superconductors. The twodimensional Hubbard model serves as the canonical microscopic physical model for such systems. Multiple orders have been proposed in the underdoped part of the phase diagram, which corresponds to a regime of maximum numerical difficulty. By combining the latest numerical methods in exhaustive simulations, we uncover the ordering in the underdoped ground state. We find 1 arXiv:1701.00054v3 [cond-mat.str-el]

We present release 1.3 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an international open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. Development is centered on common XML and binary data formats, on libraries to simplify and speed up code development, and on full-featured simulation programs. The programs enable non-experts to start carrying out numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), as well as the density matrix renormalization group (DMRG). Changes in the new release include a DMRG program for interacting models, support for translation symmetries in the diagonalization programs, the ability to define custom measurement operators, and support for inhomogeneous systems, such as lattice models with traps. The software is available from our web server at http://alps.comp-phys.org/.

We explain how to implement, in the context of projected entangled-pair states ͑PEPSs͒, the general procedure of fermionization of a tensor network introduced in P. Corboz and G. Vidal, Phys. Rev. B 80, 165129 ͑2009͒. The resulting fermionic PEPS, similar to previous proposals, can be used to study the ground state of interacting fermions on a two-dimensional lattice. As in the bosonic case, the cost of simulations depends on the amount of entanglement in the ground state and not directly on the strength of interactions. The present formulation of fermionic PEPS leads to a straightforward numerical implementation that allowed us to recycle much of the code for bosonic PEPS. We demonstrate that fermionic PEPS are a useful variational ansatz for interacting fermion systems by computing approximations to the ground state of several models on an infinite lattice. For a model of interacting spinless fermions, ground state energies lower than Hartree-Fock results are obtained, shifting the boundary between the metal and charge-density wave phases. For the t-J model, energies comparable with those of a specialized Gutzwiller-projected ansatz are also obtained.

The main characteristic of Mott insulators, as compared to band insulators, is to host low-energy spin fluctuations. In addition, Mott insulators often possess orbital degrees of freedom when crystal-field levels are partially filled. While in the majority of Mott insulators, spins and orbitals develop long-range order, the possibility for the ground state to be a quantum liquid opens new perspectives. In this paper, we provide clear evidence that the spin-orbital SUð4Þ symmetric Kugel-Khomskii model of Mott insulators on the honeycomb lattice is a quantum spin-orbital liquid. The absence of any form of symmetry breakinglattice or SUðNÞ-is supported by a combination of semiclassical and numerical approaches: flavor-wave theory, tensor network algorithm, and exact diagonalizations. In addition, all properties revealed by these methods are very accurately accounted for by a projected variational wave function based on the -flux state of fermions on the honeycomb lattice at 1=4 filling. In that state, correlations are algebraic because of the presence of a Dirac point at the Fermi level, suggesting that the symmetric Kugel-Khomskii model on the honeycomb lattice is an algebraic quantum spin-orbital liquid. This model provides an interesting starting point to understanding the recently discovered spin-orbital-liquid behavior of Ba 3 CuSb 2 O 9 . The present results also suggest the choice of optical lattices with honeycomb geometry in the search for quantum liquids in ultracold four-color fermionic atoms.

Spinless fermions on a honeycomb lattice provide a minimal realization of lattice Dirac fermions. Repulsive interactions between nearest neighbors drive a quantum phase transition from a Dirac semimetal to a charge-density-wave state through a fermionic quantum critical point, where the coupling of the Ising order parameter to the Dirac fermions at low energy drastically affects the quantum critical behavior. Encouraged by a recent discovery (Huffman and Chandrasekharan 2014 Phys. Rev. B 89 111101) of the absence of the fermion sign problem in this model, we study the fermionic quantum critical point using the continuous-time quantum Monte Carlo method with a worm-sampling technique. We estimate the transition point = V t 1.356(1) with the critical exponents ν = 0.80(3) and η = 0.302 (7). Compatible results for the transition point are also obtained with infinite projected entangled-pair states.

We present a scheme to perform an iterative variational optimization with infinite projected entangled-pair states (iPEPS), a tensor network ansatz for a two-dimensional wave function in the thermodynamic limit, to compute the ground state of a local Hamiltonian. The method is based on a systematic summation of Hamiltonian contributions using the corner transfer-matrix method. Benchmark results for challenging problems are presented, including the 2D Heisenberg model, the Shastry-Sutherland model, and the t-J model, which show that the variational scheme yields considerably more accurate results than the previously best imaginary time evolution algorithm, with a similar computational cost and with a faster convergence towards the ground state.Comment: 11 pages, 9 figures, revised (published) version, correction in Fig.

The infinite Projected Entangled Pair States (iPEPS) algorithm [J. Jordan et al, PRL 101, 250602 (2008)] has become a useful tool in the calculation of ground state properties of 2d quantum lattice systems in the thermodynamic limit. Despite its many successful implementations, the method has some limitations in its present formulation which hinder its application to some highly-entangled systems. The purpose of this paper is to unravel some of these issues, in turn enhancing the stability and efficiency of iPEPS methods. For this, we first introduce the fast full update scheme, where effective environment and iPEPS tensors are both simultaneously updated (or evolved) throughout time. As we shall show, this implies two crucial advantages: (i) dramatic computational savings, and (ii) improved overall stability. Besides, we extend the application of the local gauge fixing, successfully implemented for finite-size PEPS [M. Lubasch, J. Ignacio Cirac, M.-C. Bañuls, PRB 90, 064425 (2014)], to the iPEPS algorithm. We see that the gauge fixing not only further improves the stability of the method, but also accelerates the convergence of the alternating least squares sweeping in the (either "full" or "fast full") tensor update scheme. The improvement in terms of computational cost and stability of the resulting "improved" iPEPS algorithm is benchmarked by studying the ground state properties of the quantum Heisenberg and transverse-field Ising models on an infinite square lattice.