
Making use of infinite projected entangled pair states, we investigate the
ground state phase diagram of the nearestneighbor spin1 bilinearbiquadratic
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 onedimensional Haldane phase does not reach
all the way up to the twodimensional limit. Finally, we investigate the
possibility of an intermediate partiallymagnetic partiallyquadrupolar 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 entangledpair states (iPEPS). We identify the presence of a finite
correlation length in the optimal iPEPS approximation to Lorentzinvariant
critical states which we use to perform a finite correlationlength scaling
(FCLS) analysis to determine critical exponents. This is analogous to the
onedimensional (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 spin1 bilinearbiquadratic 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$
bilinearbiquadratic Heisenberg model on the square lattice reveal an emergent
Haldane phase in between the previously predicted antiferromagnetic and
3sublattice 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 onedimensional nature matrix product states have
been systematically used to obtain remarkably accurate results for
twodimensional systems. Motivated by basic entropic arguments favoring
projected entangledpair states as the method of choice, we assess the relative
performance of infinite matrix product states and infinite projected
entangledpair states on cylindrical geometries. By considering the Heisenberg
and halffilled 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 entangledpair states at a width of about 11 sites, whereas for the
halffilled Hubbard model this crossover occurs at about 7 sites.

We present a scheme to perform an iterative variational optimization with
infinite projected entangledpair states (iPEPS), a tensor network ansatz for a
twodimensional 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 transfermatrix method. Benchmark
results for challenging problems are presented, including the 2D Heisenberg
model, the ShastrySutherland model, and the tJ 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 hightemperature 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 conjugategradient method for the groundstate optimization of
projected entangledpair 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
imaginarytime evolution.

An infinite projected entangledpair 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 stateoftheart
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 meanfield 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 meanfield 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
highlyentangled 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 finitesize 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 transversefield 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 continuoustime
expansion of imaginarytime 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
continuoustime projection in detail using the spinless $tV$ model as an
example and compare the numerical results with exact diagonalization,
densitymatrixrenormalizationgroup and infinite projected entangledpair
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 KitaevHeisenberg model
using infinite Projected EntangledPair 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 KitaevHeisenberg 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 symmetrybroken phases.

Variational studies of the tJ model on the square lattice based on infinite
projectedentangled pair states (iPEPS) confirm an extremely close competition
between a uniform dwave superconducting state and different stripe states. The
sitecentered stripe with an inphase dwave order has an equal or only
slightly lower energy than the stripe with antiphase dwave 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 chargedensitywave
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 entangledpair
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
holedoped resonating valence bond (RVB) wavefunctions with shortrange 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 squareroot 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 $tJ_1J_2$ model (including
both nearest and nextnearest neighbor antiferromagnetic interactions $J_1$ and
$J_2$, respectively) for $J_2\simeq J_1/2$ and obtain an energy very close to
the infinitePEPS state (using full update optimization and same bond
dimension). The orbital symmetry of the optimized RVB superconductor has
predominant dwave character, although we argue a residual (complex swave)
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 entangledpair states (iPEPS), we show that the
ShastrySutherland model in an external magnetic field has lowmagnetization
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 orthogonaldimer 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 domainwall 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 threestate 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 3color
fermions, while the parameter range J_perp<min(0,J_parallel) can be realized
in the Mott insulating phase of 3color bosonic atoms. Using linear flavor wave
theory, infinite projected entangledpair states (iPEPS), and continuoustime
quantum MonteCarlo 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 threesublattice longrange 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 threesublattice superfluid
followed by a twosublattice 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 KosterlitzThouless 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 flavorwave 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 projectedentangled pair states (iPEPS), for
which the accuracy can be systematically controlled by the socalled 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 nonchiral flux
states.

We simulate the ShastrySutherland model in two dimensions by means of
infinite projected entangledpair states (iPEPS)  a variational tensor network
method where the accuracy can be systematically controlled by the socalled
bond dimension. Besides the well established dimer and N\'eel phase iPEPS
confirms the presence of an intermediate phase with plaquette longrange 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 firstorder phase
transition between plaquette and N\'eel phase occurring for $J=0.765(15)$. We
do not find a stable intermediate columnardimer 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 entangledpair states (iPEPS) and exact diagonalizations, we
show that in both cases the ground state is a simplex solid state with a
twofold 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 meanfield 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 entangledpair states
(iPEPS) and exact diagonalization.

We present a numerical study of the SU(3) Heisenberg model of threeflavor
fermions on the triangular and square lattice by means of the densitymatrix
renormalization group (DMRG) and infinite projected entangledpair states
(iPEPS). For the triangular lattice we confirm that the ground state has a
threesublattice 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 threesublattice order, with an ordered moment of
m=0.20.4 in the thermodynamic limit.

Using infinite projected entangled pair states (iPEPS), exact
diagonalization, and flavorwave 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 threesublattice
stripelike longrange order. Since the ground state of a dimer is not a
singlet for SU(4) but a 6dimensional 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.