We study the 3D Kitaev and Kitaev-Heisenberg models respectively on the hyperhoneycomb and hyperoctagon lattices, both at zero and finite-temperature, in the thermodynamic limit. Our analysis relies on advanced tensor network (TN) simulations based o
n graph Projected Entangled-Pair States (gPEPS). We map out the TN phase diagrams of the models and characterize their underlying gapped and gapless phases both at zero and finite temperature. In particular, we demonstrate how cooling down the hyperhoneycomb system from high-temperature leads to fractionalization of spins to itinerant Majorana fermions and gauge fields that occurs in two separate temperature regimes, leaving their fingerprint on specific heat as a double-peak feature as well as on other quantities such as the thermal entropy, spin-spin correlations and bond entropy. Using the Majorana representation of the Kitaev model, we further show that the low-temperature thermal transition to the Kitaev quantum spin liquid (QSL) phase is associated with the non-trivial Majorana band topology and the presence of Weyl nodes, which manifests itself via non-vanishing Chern number and finite thermal Hall conductivity. Beyond the pure Kitaev limit, we study the 3D Kitaev-Heisenberg (KH) model on the hyperoctagon lattice and extract the full phase diagram for different Heisenberg couplings. We further explore the thermodynamic properties of the magnetically-ordered regions in the KH model and show that, in contrast to the QSL phase, here the thermal phase transition follows the standard Landau symmetry-breaking theory.
Ultracold atoms in optical lattices are one of the most promising experimental setups to simulate strongly correlated systems. However, efficient numerical algorithms able to benchmark experiments at low-temperatures in interesting 3d lattices are la
cking. To this aim, here we introduce an efficient tensor network algorithm to accurately simulate thermal states of local Hamiltonians in any infinite lattice, and in any dimension. We apply the method to simulate thermal bosons in optical lattices. In particular, we study the physics of the (soft-core and hard-core) Bose-Hubbard model on the infinite pyrochlore and cubic lattices with unprecedented accuracy. Our technique is therefore an ideal tool to benchmark realistic and interesting optical-lattice experiments.
We study the zero-temperature phase diagram of the spin-$frac{1}{2}$ Heisenberg model with breathing anisotropy (i.e., with different coupling strength on the upward and downward triangles) on the kagome lattice. Our study relies on large scale tenso
r network simulations based on infinite projected entangled-pair state and infinite projected entangled-simplex state methods adapted to the kagome lattice. Our energy analysis suggests that the U(1) algebraic quantum spin-liquid (QSL) ground-state of the isotropic Heisenberg model is stable up to very large breathing anisotropy until it breaks down to a critical lattice-nematic phase that breaks rotational symmetry in real space through a first-order quantum phase transition. Our results also provide further insight into the recent experiment on vanadium oxyfluoride compounds which has been shown to be relevant platforms for realizing QSL in the presence of breathing anisotropy.
We construct a short-range resonating valence-bond state (RVB) on the ruby lattice, using projected entangled-pair states (PEPS) with bond dimension $D=3$. By introducing non-local moves to the dimer patterns on the torus, we distinguish four distinc
t sectors in the space of dimer coverings, which is a signature of the topological nature of the RVB wave function. Furthermore, by calculating the reduced density matrix of a bipartition of the RVB state on an infinite cylinder and exploring its entanglement entropy, we confirm the topological nature of the RVB wave function by obtaining non-zero topological contribution, $gamma=-rm{ln} 2$, consistent with that of a $mathbb{Z}_2$ topological quantum spin liquid. We also calculate the ground-state energy of the spin-$frac{1}{2}$ antiferromagnetic Heisenberg model on the ruby lattice and compare it with the RVB energy. Finally, we construct a quantum-dimer model for the ruby lattice and discuss it as a possible parent Hamiltonian for the RVB wave function.
We develop a strategy for tensor network algorithms that allows to deal very efficiently with lattices of high connectivity. The basic idea is to fine-grain the physical degrees of freedom, i.e., decompose them into more fundamental units which, afte
r a suitable coarse-graining, provide the original ones. Thanks to this procedure, the original lattice with high connectivity is transformed by an isometry into a simpler structure, which is easier to simulate via usual tensor network methods. In particular this enables the use of standard schemes to contract infinite 2d tensor networks - such as Corner Transfer Matrix Renormalization schemes - which are more involved on complex lattice structures. We prove the validity of our approach by numerically computing the ground-state properties of the ferromagnetic spin-1 transverse-field Ising model on the 2d triangular and 3d stacked triangular lattice, as well as of the hard-core and soft-core Bose-Hubbard models on the triangular lattice. Our results are benchmarked against those obtained with other techniques, such as perturbative continuous unitary transformations and graph projected entangled pair states, showing excellent agreement and also improved performance in several regimes.
We present a general graph-based Projected Entangled-Pair State (gPEPS) algorithm to approximate ground states of nearest-neighbor local Hamiltonians on any lattice or graph of infinite size. By introducing the structural-matrix which codifies the de
tails of tensor networks on any graphs in any dimension $d$, we are able to produce a code that can be essentially launched to simulate any lattice. We further introduce an optimized algorithm to compute simple tensor updates as well as expectation values and correlators with a mean-field-like effective environments. Though not being variational, this strategy allows to cope with PEPS of very large bond dimension (e.g., $D=100$), and produces remarkably accurate results in the thermodynamic limit in many situations, and specially when the correlation length is small and the connectivity of the lattice is large. We prove the validity of our approach by benchmarking the algorithm against known results for several models, i.e., the antiferromagnetic Heisenberg model on a chain, star and cubic lattices, the hardcore Bose-Hubbard model on square lattice, the ferromagnetic Heisenberg model in a field on the pyrochlore lattice, as well as the $3$-state quantum Potts model in field on the kagome lattice and the spin-$1$ bilinear-biquadratic Heisenberg model on the triangular lattice. We further demonstrate the performance of gPEPS by studying the quantum phase transition of the $2d$ quantum Ising model in transverse magnetic field on the square lattice, and the phase diagram of the Kitaev-Heisenberg model on the hyperhoneycomb lattice. Our results are in excellent agreement with previous studies.
The infinite Projected Entangled-Pair State (iPEPS) algorithm is one of the most efficient techniques for studying the ground-state properties of two-dimensional quantum lattice Hamiltonians in the thermodynamic limit. Here, we show how the algorithm
can be adapted to explore nearest-neighbor local Hamiltonians on the ruby and triangle-honeycomb lattices, using the Corner Transfer Matrix (CTM) renormalization group for 2D tensor network contraction. Additionally, we show how the CTM method can be used to calculate the ground state fidelity per lattice site and the boundary density operator and entanglement entropy (EE) on an infinite cylinder. As a benchmark, we apply the iPEPS method to the ruby model with anisotropic interactions and explore the ground-state properties of the system. We further extract the phase diagram of the model in different regimes of the couplings by measuring two-point correlators, ground state fidelity and EE on an infinite cylinder. Our phase diagram is in agreement with previous studies of the model by exact diagonalization.
The ruby lattice is a four-valent lattice interpolating between honeycomb and triangular lattices. In this work we investigate the topological spin-liquid phases of a spin Hamiltonian with Kitaev interactions on the ruby lattice using exact diagonali
zation and perturbative methods. The latter interactions combined with the structure of the lattice yield a model with $mathbb{Z}_2 times mathbb{Z}_2$ gauge symmetry. We mapped out the phase digram of the model and found gapped and gapless spin-liquid phases. While the low energy sector of the gapped phase corresponds to the well-known topological color code model on a honeycomb lattice, the low-energy sector of the gapless phases is described by an effective spin model with three-body interactions on a triangular lattice. A gap is opened in the spectrum in a small magnetic field. We argue that the latter phases could be possibly described by exotic excitations, whose their spectrum is richer than the Ising phase of the Kitaev model.
We use the topological entanglement entropy (TEE) as an efficient tool to fully characterize the Abelian phase of a $mathbb{Z}_2 times mathbb{Z}_2$ spin liquid emerging as the ground state of topological color code (TCC), which is a class of stabiliz
er states on the honeycomb lattice. We provide the fusion rules of the quasiparticle (QP) excitations of the model by introducing single- or two-body operators on physical spins for each fusion process which justify the corresponding fusion outcome. Beside, we extract the TEE from Renyi entanglement entropy (EE) of the TCC, analytically and numerically by finite size exact diagonalization on the disk shape regions with contractible boundaries. We obtain that the EE has a local contribution, which scales linearly with the boundary length in addition to a topological term, i.e. the TEE, arising from the condensation of closed strings in the ground state. We further investigate the ground state dependence of the TEE on regions with non-contractible boundaries, i.e. by cutting the torus to half cylinders, from which we further identify multiple independent minimum entropy states (MES) of the TCC and then extract the U and S modular matrices of the system, which contain the self and mutual statistics of the anyonic QPs and fully characterize the topological phase of the TCC. Eventually, we show that, in spite of the lack of a local order parameter, TEE and other physical quantities obtained from ground state wave function such as entanglement spectrum (ES) and ground state fidelity are sensitive probes to study the robustness of a topological phase. We find that the topological order in the presence of a magnetic field persists until the vicinity of the transition point, where the TEE and fidelity drops to zero and the ES splits severely, signaling breakdown of the topological phase of the TCC.
