Simulation of fermionic lattice models in two dimensions with Projected Entangled-Pair States: Next-nearest neighbor Hamiltonians

In a recent contribution [Phys. Rev. B 81, 165104 (2010)] fermionic Projected Entangled-Pair States (PEPS) were used to approximate the ground state of free and interacting spinless fermion models, as well as the $t$-$J$ model. This paper revisits these three models in the presence of an additional next-nearest hopping amplitude in the Hamiltonian. First we explain how to account for next-nearest neighbor Hamiltonian terms in the context of fermionic PEPS algorithms based on simulating time evolution. Then we present benchmark calculations for the three models of fermions, and compare our results against analytical, mean-field, and variational Monte Carlo results, respectively. Consistent with previous computations restricted to nearest-neighbor Hamiltonians, we systematically obtain more accurate (or better converged) results for gapped phases than for gapless ones.Comment: 10 pages, 11 figures, minor change

Simulation of strongly correlated fermions in two spatial dimensions with fermionic Projected Entangled-Pair States

We explain how to implement, in the context of projected entangled-pair states (PEPS), the general procedure of fermionization of a tensor network introduced in [P. Corboz, 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.Comment: 25 pages, 35 figures (revised version

Simulation of anyons with tensor network algorithms

Interacting systems of anyons pose a unique challenge to condensed matter simulations due to their non-trivial exchange statistics. These systems are of great interest as they have the potential for robust universal quantum computation, but numerical tools for studying them are as yet limited. We show how existing tensor network algorithms may be adapted for use with systems of anyons, and demonstrate this process for the 1-D Multi-scale Entanglement Renormalisation Ansatz (MERA). We apply the MERA to infinite chains of interacting Fibonacci anyons, computing their scaling dimensions and local scaling operators. The scaling dimensions obtained are seen to be in agreement with conformal field theory. The techniques developed are applicable to any tensor network algorithm, and the ability to adapt these ansaetze for use on anyonic systems opens the door for numerical simulation of large systems of free and interacting anyons in one and two dimensions.Comment: Fixed typos, matches published version. 16 pages, 21 figures, 4 tables, RevTeX 4-1. For a related work, see arXiv:1006.247

Magnetization of SrCu2(BO3)2 in ultrahigh magnetic fields up to 118 T

The magnetization process of the orthogonal-dimer antiferromagnet SrCu2(BO3)2 is investigated in high magnetic fields of up to 118 T. A 1/2 plateau is clearly observed in the field range 84 to 108 T in addition to 1/8, 1/4 and 1/3 plateaux at lower fields. Using a combination of state-of-the-art numerical simulations, the main features of the high-field magnetization, a 1/2 plateau of width 24 T, a 1/3 plateau of width 34 T, and no 2/5 plateau, are shown to agree quantitatively with the Shastry-Sutherland model if the ratio of inter- to intra-dimer exchange interactions J'/J=0.63. It is further predicted that the intermediate phase between the 1/3 and 1/2 plateau is not uniform but consists of a 1/3 supersolid followed by a 2/5 supersolid and possibly a domain-wall phase, with a reentrance into the 1/3 supersolid above the 1/2 plateau.Comment: 5 pages + 10 pages supplemental materia

Stripe order in the underdoped region of the two-dimensional Hubbard model

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

Comment on "Topological quantum phase transitions of attractive spinless fermions in a honeycomb lattice" by Poletti D. et al

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.Comment: 2 pages, 1 figur

Entanglement Entropy of Random Fractional Quantum Hall Systems

The entanglement entropy of the $\nu = 1/3$ and $\nu = 5/2$ quantum Hall states in the presence of short range random disorder has been calculated by direct diagonalization. A microscopic model of electron-electron interaction is used, electrons are confined to a single Landau level and interact with long range Coulomb interaction. For very weak disorder, the values of the topological entanglement entropy are roughly consistent with expected theoretical results. By considering a broader range of disorder strengths, the fluctuation in the entanglement entropy was studied in an effort to detect quantum phase transitions. In particular, there is a clear signature of a transition as a function of the disorder strength for the $\nu = 5/2$ state. Prospects for using the density matrix renormalization group to compute the entanglement entropy for larger system sizes are discussed.Comment: 29 pages, 16 figures; fixed figures and figure captions; revised fluctuation calculation

Entanglement renormalization and boundary critical phenomena

The multiscale entanglement renormalization ansatz is applied to the study of boundary critical phenomena. We compute averages of local operators as a function of the distance from the boundary and the surface contribution to the ground state energy. Furthermore, assuming a uniform tensor structure, we show that the multiscale entanglement renormalization ansatz implies an exact relation between bulk and boundary critical exponents known to exist for boundary critical systems.Comment: 6 pages, 4 figures; for a related work see arXiv:0912.164