• Making use of infinite projected entangled pair states, we investigate the ground state phase diagram of the nearest-neighbor spin-1 bilinear-biquadratic Heisenberg model on the triangular lattice. In agreement with previous studies, we find the ferromagnetic, 120 degree magnetically ordered, ferroquadrupolar and antiferroquadrupolar phases, and confirm that all corresponding phase transitions are first order. Moreover, we provide an accurate estimate of the location of the ferroquadrupolar to 120 degree magnetically ordered phase transition, thereby fully establishing the phase diagram. Also, we do not encounter any signs of the existence of a quantum paramagnetic phase. In particular, contrary to the equivalent square lattice model, we demonstrate that on the triangular lattice the one-dimensional Haldane phase does not reach all the way up to the two-dimensional limit. Finally, we investigate the possibility of an intermediate partially-magnetic partially-quadrupolar phase close to $\theta = \pi/2$, and we show that, also contrary to the square lattice case, this phase is not present on the triangular lattice.
  • We show how to accurately study 2D quantum critical phenomena using infinite projected entangled-pair states (iPEPS). We identify the presence of a finite correlation length in the optimal iPEPS approximation to Lorentz-invariant critical states which we use to perform a finite correlation-length scaling (FCLS) analysis to determine critical exponents. This is analogous to the one-dimensional (1D) finite entanglement scaling with infinite matrix product states. We provide arguments why this approach is also valid in 2D by identifying a class of states that despite obeying the area law of entanglement seems hard to describe with iPEPS. We apply these ideas to interacting spinless fermions on a honeycomb lattice and obtain critical exponents which are in agreement with Quantum Monte Carlo results. Furthermore, we introduce a new scheme to locate the critical point without the need of computing higher order moments of the order parameter. Finally, we also show how to obtain an improved estimate of the order parameter in gapless systems, with the 2D Heisenberg model as an example.
  • Using infinite projected entangled pair states, we study the ground state phase diagram of the spin-1 bilinear-biquadratic Heisenberg model on the square lattice directly in the thermodynamic limit. We find an unexpected partially nematic partially magnetic phase in between the antiferroquadrupolar and ferromagnetic regions. Furthermore, we describe all observed phases and discuss the nature of the phase transitions involved.
  • Infinite projected entangled pair states simulations of the $S=1$ bilinear-biquadratic Heisenberg model on the square lattice reveal an emergent Haldane phase in between the previously predicted antiferromagnetic and 3-sublattice 120$^\circ$ magnetically ordered phases. This intermediate phase preserves SU(2) spin and translational symmetry but breaks lattice rotational symmetry, and it can be adiabatically connected to the Haldane phase of decoupled $S=1$ chains. Our results contradict previous studies which found a direct transition between the two magnetically ordered states.
  • In spite of their intrinsic one-dimensional nature matrix product states have been systematically used to obtain remarkably accurate results for two-dimensional systems. Motivated by basic entropic arguments favoring projected entangled-pair states as the method of choice, we assess the relative performance of infinite matrix product states and infinite projected entangled-pair states on cylindrical geometries. By considering the Heisenberg and half-filled Hubbard models on the square lattice as our benchmark cases, we evaluate their variational energies as a function of both bond dimension as well as cylinder width. In both examples we find crossovers at moderate cylinder widths, i.e. for the largest bond dimensions considered we find an improvement on the variational energies for the Heisenberg model by using projected entangled-pair states at a width of about 11 sites, whereas for the half-filled Hubbard model this crossover occurs at about 7 sites.
  • 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.
  • Competing inhomogeneous orders are a central feature of correlated electron materials including the high-temperature superconductors. The two- dimensional 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 a stripe order that has a highly compressible wavelength on an energy scale of a few Kelvin, with wavelength fluctuations coupled to pairing order. The favored filled stripe order is different from that seen in real materials. Our results demonstrate the power of modern numerical methods to solve microscopic models even in challenging settings.
  • We present a conjugate-gradient method for the ground-state optimization of projected entangled-pair states (PEPS) in the thermodynamic limit, as a direct implementation of the variational principle within the PEPS manifold. Our optimization is based on an efficient and accurate evaluation of the gradient of the global energy functional by using effective corner environments, and is robust with respect to the initial starting points. It has the additional advantage that physical and virtual symmetries can be straightforwardly implemented. We provide the tools to compute static structure factors directly in momentum space, as well as the variance of the Hamiltonian. We benchmark our method on Ising and Heisenberg models, and show a significant improvement on the energies and order parameters as compared to algorithms based on imaginary-time evolution.
  • An infinite projected entangled-pair state (iPEPS) is a variational tensor network ansatz for 2D wave functions in the thermodynamic limit where the accuracy can be systematically controlled by the bond dimension $D$. We show that for the doped Hubbard model in the strongly correlated regime ($U/t=8$, $n=0.875$) iPEPS yields lower variational energies than state-of-the-art variational methods in the large 2D limit, which demonstrates the competitiveness of the method. In order to obtain an accurate estimate of the energy in the exact infinite $D$ limit we introduce and test an extrapolation technique based on a truncation error computed in the iPEPS imaginary time evolution algorithm. The extrapolated energies are compared with accurate quantum Monte Carlo results at half filling and with various other methods in the doped, strongly correlated regime.
  • We revisit the SU(6) Heisenberg model on the honeycomb lattice, which has been predicted to be a chiral spin liquid by mean-field theory [G. Szirmai et al., Phys. Rev. A 84, 011611 (2011)]. Using exact diagonalizations of finite clusters, infinite projected entangled pair states simulations, and variational Monte Carlo simulations based on Gutzwiller projected wave functions, we provide strong evidence in favour of the competing plaquette state, which was reported to be higher but close by in energy according to mean-field theory. This is further confirmed by the investigation of the model with a ring exchange term, which shows that there is a transition between the plaquette state and the chiral state at a finite value of the ring exchange term.
  • 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 \emph{local gauge fixing}, successfully implemented for finite-size PEPS [M. Lubasch, J. Ignacio Cirac, M.-C. Ba\~{n}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.
  • We present the ground state extension of the efficient quantum Monte Carlo algorithm for lattice fermions of arXiv:1411.0683. Based on continuous-time expansion of imaginary-time projection operator, the algorithm is free of systematic error and scales \emph{linearly} with projection time and interaction strength. Compared to the conventional quantum Monte Carlo methods for lattice fermions, this approach has greater flexibility and is easier to combine with powerful machinery such as histogram reweighting and extended ensemble simulation techniques. We discuss the implementation of the continuous-time projection in detail using the spinless $t-V$ model as an example and compare the numerical results with exact diagonalization, density-matrix-renormalization-group and infinite projected entangled-pair states calculations. Finally we use the method to study the fermionic quantum critical point of spinless fermions on a honeycomb lattice and confirm previous results concerning its critical exponents.
  • We study the extent of the spin liquid phases in the Kitaev-Heisenberg model using infinite Projected Entangled-Pair States tensor network ansatz wave functions directly in the thermodynamic limit. To assess the accuracy of the ansatz wave functions we perform benchmarks against exact results for the Kitaev model and find very good agreement for various observables. In the case of the Kitaev-Heisenberg model we confirm the existence of 6 different phases: N\'eel, stripy, ferromagnetic, zigzag and two spin liquid phases. We find finite extents for both spin liquid phases and discontinuous phase transitions connecting them to symmetry-broken phases.
  • Variational studies of the t-J model on the square lattice based on infinite projected-entangled pair states (iPEPS) 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 anti-phase 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.
  • 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 Ising order parameter to the Dirac fermions at low energy drastically affects the quantum critical behavior. Encouraged by a recently discovery of 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 worm sampling technique. We estimate the transition point $V/t= 1.356(1)$ with the critical exponents $\nu =0.80(3)$ and $\eta =0.302(7)$. Compatible results for the transition point are also obtained with infinite projected entangled-pair states.
  • We construct a family of simple fermionic projected entangled pair states (fPEPS) on the square lattice with bond dimension $D=3$ which are exactly hole-doped resonating valence bond (RVB) wavefunctions with short-range singlet bonds. Under doping the insulating RVB spin liquid evolves immediately into a superconductor with mixed $d+is$ pairing symmetry whose pair amplitude grows as the square-root of the doping. The relative weight between $s$-wave and $d$-wave components can be controlled by a single variational parameter $c$. We optimize our ansatz w.r.t. $c$ for the frustrated $t-J_1-J_2$ model (including both nearest and next-nearest neighbor antiferromagnetic interactions $J_1$ and $J_2$, respectively) for $J_2\simeq J_1/2$ and obtain an energy very close to the infinite-PEPS state (using full update optimization and same bond dimension). The orbital symmetry of the optimized RVB superconductor has predominant d-wave character, although we argue a residual (complex s-wave) time reversal symmetry breaking component should always be present. Connections of the results to the physics of superconducting cuprates and pnictides are outlined.
  • Using infinite projected entangled-pair states (iPEPS), we show that the Shastry-Sutherland model in an external magnetic field has low-magnetization plateaus which, in contrast to previous predictions, correspond to crystals of bound states of triplets, and {\it not} to crystals of triplets. The first sizable plateaus appear at magnetization 1/8, 2/15 and 1/6, in agreement with experiments on the orthogonal-dimer antiferromagnet SrCu$_2$(BO$_3$)$_2$, and they can be naturally understood as regular patterns of bound states, including the intriguing 2/15 one. We also show that, even in a confined geometry, two triplets bind into a {\it localized} bound state with $S_z=2$. Finally, we discuss the role of competing domain-wall and supersolid phases as well as that of additional anisotropic interactions.
  • The low temperature phase diagram of $^4$He adsorbed on a single graphene sheet is studied by computer simulation of a system comprising nearly thousand helium atoms. In the first layer, two commensurate solid phases are observed, with fillings 1/3 and 7/16 respectively, separated by a domain wall phase, as well as an incommensurate crystal at higher coverage. No evidence of a thermodynamically stable superfliuid phase is found for the first adlayer. Second layer promotion occurs at a coverage of 0.111(4) $\AA^{-2}$. In the second layer two phases are observed, namely a superfluid and an incommensurate solid, with no commensurate solid intervening between these two phases. The computed phase diagram closely resembles that predicted for helium on graphite.
  • We study a quantum version of the three-state Potts model that includes as special cases the effective models of bosons and fermions on the square lattice in the Mott insulating limit. It can be viewed as a model of quantum permutations with amplitudes J_parallel and J_perp for identical and different colors, respectively. For J_parallel=J_perp>0, it is equivalent to the SU(3) Heisenberg model, which describes the Mott insulating phase of 3-color fermions, while the parameter range J_perp<min(0,-J_parallel) can be realized in the Mott insulating phase of 3-color bosonic atoms. Using linear flavor wave theory, infinite projected entangled-pair states (iPEPS), and continuous-time quantum Monte-Carlo simulations, we construct the full T=0 phase diagram, and we explore the T>0 properties for J_perp<0. For dominant antiferromagnetic J_parallel interactions, a three-sublattice long-range ordered stripe state is selected out of the ground state manifold of the antiferromagnetic Potts model by quantum fluctuations. Upon increasing |J_perp|, this state is replaced by a uniform superfluid for J_perp<0, and by an exotic three-sublattice superfluid followed by a two-sublattice superfluid for J_perp>0. The transition out of the uniform superfluid (that can be realized with bosons) is shown to be a peculiar type of Kosterlitz-Thouless transition with three types of elementary vortices.
  • Conflicting predictions have been made for the ground state of the SU(3) Heisenberg model on the honeycomb lattice: Tensor network simulations found a plaquette order [Zhao et al, Phys. Rev. B 85, 134416 (2012)], where singlets are formed on hexagons, while linear flavor-wave theory (LFWT) suggested a dimerized, color ordered state [Lee and Yang, Phys. Rev. B 85, 100402 (2012)]. In this work we show that the former state is the true ground state by a systematic study with infinite projected-entangled pair states (iPEPS), for which the accuracy can be systematically controlled by the so-called bond dimension $D$. Both competing states can be reproduced with iPEPS by using different unit cell sizes. For small $D$ the dimer state has a lower variational energy than the plaquette state, however, for large $D$ it is the latter which becomes energetically favorable. The plaquette formation is also confirmed by exact diagonalizations and variational Monte Carlo studies, according to which both the dimerized and plaquette states are non-chiral flux states.
  • We simulate the Shastry-Sutherland model in two dimensions by means of infinite projected entangled-pair states (iPEPS) - a variational tensor network method where the accuracy can be systematically controlled by the so-called bond dimension. Besides the well established dimer and N\'eel phase iPEPS confirms the presence of an intermediate phase with plaquette long-range order, and we determine its phase boundaries with high accuracy. The first order phase transition for $J=0.675(2)$ between dimer and plaquette phase is compatible with previous series expansion results. iPEPS predicts a weak first-order phase transition between plaquette and N\'eel phase occurring for $J=0.765(15)$. We do not find a stable intermediate columnar-dimer phase, even when we bias the state towards this order.
  • We present a numerical study of the SU(N) Heisenberg model with the fundamental representation at each site for the kagome lattice (for N=3) and the checkerboard lattice (for N=4), which are the line graphs of the honeycomb and square lattices and thus belong to the class of bisimplex lattices. Using infinite projected entangled-pair states (iPEPS) and exact diagonalizations, we show that in both cases the ground state is a simplex solid state with a two-fold ground state degeneracy, in which the N spins belonging to a simplex (i.e. a complete graph) form a singlet. Theses states can be seen as generalizations of valence bond solid states known to be stabilized in certain SU(2) spin models.
  • In a recent letter [D. Poletti et al., EPL 93, 37008 (2011)] a model of attractive spinless fermions on the honeycomb lattice at half filling has been studied by mean-field theory, where distinct homogenous phases at rather large attraction strength $V>3.36$, separated by (topological) phase transitions, have been predicted. In this comment we argue that without additional interactions the ground states in these phases are not stable against phase separation. We determine the onset of phase separation at half filling $V_{ps}\approx 1.7$ by means of infinite projected entangled-pair states (iPEPS) and exact diagonalization.
  • We present a numerical study of the SU(3) Heisenberg model of three-flavor fermions on the triangular and square lattice by means of the density-matrix renormalization group (DMRG) and infinite projected entangled-pair states (iPEPS). For the triangular lattice we confirm that the ground state has a three-sublattice order with a finite ordered moment which is compatible with the result from linear flavor wave theory (LFWT). The same type of order has recently been predicted also for the square lattice [PRL 105, 265301 (2010)] from LFWT and exact diagonalization. However, for this case the ordered moment cannot be computed based on LFWT due to divergent fluctuations. Our numerical study clearly supports this three-sublattice order, with an ordered moment of m=0.2-0.4 in the thermodynamic limit.
  • Using infinite projected entangled pair states (iPEPS), exact diagonalization, and flavor-wave theory, we show that the SU(4) Heisenberg model undergoes a spontaneous dimerization on the square lattice, in contrast to its SU(2) and SU(3) counterparts, which develop Neel and three-sublattice stripe-like long-range order. Since the ground state of a dimer is not a singlet for SU(4) but a 6-dimensional irrep, this leaves the door open for further symmetry breaking. We provide evidence that, unlike in SU(4) ladders, where dimers pair up to form singlet plaquettes, here the SU(4) symmetry is additionally broken, leading to a gapless spectrum in spite of the broken translational symmetry.