The classical Heisenberg model in two spatial dimensions constitutes one of the most paradigmatic spin models, taking an important role in statistical and condensed matter physics to understand magnetism. Still, despite its paradigmatic character and
the widely accepted ban of a (continuous) spontaneous symmetry breaking, controversies remain whether the model exhibits a phase transition at finite temperature. Importantly, the model can be interpreted as a lattice discretization of the $O(3)$ non-linear sigma model in $1+1$ dimensions, one of the simplest quantum field theories encompassing crucial features of celebrated higher-dimensional ones (like quantum chromodynamics in $3+1$ dimensions), namely the phenomenon of asymptotic freedom. This should also exclude finite-temperature transitions, but lattice effects might play a significant role in correcting the mainstream picture. In this work, we make use of state-of-the-art tensor network approaches, representing the classical partition function in the thermodynamic limit over a large range of temperatures, to comprehensively explore the correlation structure for Gibbs states. By implementing an $SU(2)$ symmetry in our two-dimensional tensor network contraction scheme, we are able to handle very large effective bond dimensions of the environment up to $chi_E^text{eff} sim 1500$, a feature that is crucial in detecting phase transitions. With decreasing temperatures, we find a rapidly diverging correlation length, whose behaviour is apparently compatible with the two main contradictory hypotheses known in the literature, namely a finite-$T$ transition and asymptotic freedom, though with a slight preference for the second.
We study the system of multi-body interacting bosons on a two dimensional optical lattice and analyze the formation of bound bosonic pairs in the context of the Bose-Hubbard model. Assuming a repulsive two-body interaction we obtain the signatures of
pair formation in the regions between the Mott insulator lobes of the phase diagram for different choices of higher order local interactions. Considering the most general Bose-Hubbard model involving local multi-body interactions we investigate the ground state properties utilizing the cluster mean-field theory approach and further confirm the results by means of sophisticated infinite Projected Entangled Pair States calculations. By using various order parameters, we show that the choice of higher-order interaction can lead to pair superfluid phase in the system between two different Mott lobes. We also analyze the effect of temperature and density-dependent tunneling to establish the stability of the PSF phase.
Many-body localization is a striking mechanism that prevents interacting quantum systems from thermalizing. The absence of thermalization behaviour manifests itself, for example, in a remanence of local particle number configurations, a quantity that
is robust over a parameter range. Local particle numbers are directly accessible in programmable quantum simulators, in systems of cold atoms even in two spatial dimensions. Yet, the classical simulation aimed at building trust in quantum simulations is highly challenging. In this work, we present a comprehensive tensor network simulation of a many-body localized systems in two spatial dimensions using a variant of an iPEPS algorithm. The required translational invariance can be restored by implementing the disorder into an auxiliary spin system, providing an exact disorder average under dynamics. We can quantitatively assess signatures of many-body localization for the infinite system: Our methods are powerful enough to provide crude dynamical estimates for the transition between localized and ergodic phases. Interestingly, in this setting of finitely many disorder values, which we also compare with simulations involving non-interacting fermions and for which we discuss the emergent physics, localization emerges in the interacting regime, for which we provide an intuitive argument, while Anderson localization is absent.
Coupling a quantum many-body system to an external environment dramatically changes its dynamics and offers novel possibilities not found in closed systems. Of special interest are the properties of the steady state of such open quantum many-body sys
tems, as well as the relaxation dynamics towards the steady state. However, new computational tools are required to simulate open quantum many-body systems, as methods developed for closed systems cannot be readily applied. We review several approaches to simulate open many-body systems and point out the advances made in recent years towards the simulation of large system sizes.
Understanding dissipation in 2D quantum many-body systems is a remarkably difficult open challenge. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady-states of 2D quan
tum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin-1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.
We study numerically the spin-1/2 XXZ model in a field on an infinite Kagome lattice. We use different algorithms based on infinite Projected Entangled Pair States (iPEPS) for this, namely: (i) with simplex tensors and 9-site unit cell, and (ii) coar
se-graining three spins in the Kagome lattice and mapping it to a square-lattice model with nearest-neighbor interactions, with usual PEPS tensors, 6- and 12-site unit cells. Similarly to our previous calculation at the SU(2)-symmetric point (Heisenberg Hamiltonian), for any anisotropy from the Ising limit to the XY limit, we also observe the emergence of magnetization plateaus as a function of the magnetic field, at $m_z = frac{1}{3}$ using 6- 9- and 12-site PEPS unit cells, and at $m_z = frac{1}{9}, frac{5}{9}$ and $frac{7}{9}$ using a 9-site PEPS unit cell, the later set-up being able to accommodate $sqrt{3} times sqrt{3}$ solid order. We also find that, at $m_z = frac{1}{3}$, (lattice) nematic and $sqrt{3} times sqrt{3}$ VBC-order states are degenerate within the accuracy of the 9-site simplex-method, for all anisotropy. The 6- and 12-site coarse-grained PEPS methods produce almost-degenerate nematic and $1 times 2$ VBC-Solid orders. Within our accuracy, the 6-site coarse-grained PEPS method gives slightly lower energies, which can be explained by the larger amount of entanglement this approach can handle, even when the PEPS unit-cell is not commensurate with the expected ground state. Furthermore, we do not observe chiral spin liquid behaviors at and close to the XY point, as has been recently proposed. Our results are the first tensor network investigations of the XXZ spin chain in a field, and reveal the subtle competition between nearby magnetic orders in numerical simulations of frustrated quantum antiferromagnets, as well as the delicate interplay between energy optimization and symmetry in tensor networks.
Symmetry-protected trivial (SPt) phases of matter are the product-state analogue of symmetry-protected topological (SPT) phases. This means, SPt phases can be adiabatically connected to a product state by some path that preserves the protecting symme
try. Moreover, SPt and SPT phases can be adiabatically connected to each other when interaction terms that break the symmetries protecting the SPT order are added in the Hamiltonian. It is also known that spin-1 SPT phases in quantum spin chains can emerge as effective intermediate phases of spin-2 Hamiltonians. In this paper we show that a similar scenario is also valid for SPt phases. More precisely, we show that for a given spin-2 quantum chain, effective intermediate spin-1 SPt phases emerge in some regions of the phase diagram, these also being adiabatically connected to non-trivial intermediate SPT phases. We characterize the phase diagram of our model by studying quantities such as the entanglement entropy, symmetry-related order parameters, and 1-site fidelities. Our numerical analysis uses Matrix Product States (MPS) and the infinite Time-Evolving Block Decimation (iTEBD) method to approximate ground states of the system in the thermodynamic limit. Moreover, we provide a field theory description of the possible quantum phase transitions between the SPt phases. Together with the numerical results, such a description shows that the transitions may be described by Conformal Field Theories (CFT) with central charge c=1. Our results are in agreement, and further generalize, those in [Y. Fuji, F. Pollmann, M. Oshikawa, Phys. Rev. Lett. 114, 177204 (2015)].
Here we study the emergence of different Symmetry-Protected Topological (SPT) phases in a spin-2 quantum chain. We consider a Heisenberg-like model with bilinear, biquadratic, bicubic, and biquartic nearest-neighbor interactions, as well as uniaxial
anisotropy. We show that this model contains four different effective spin-1 SPT phases, corresponding to different representations of the $(mathbb{Z}_2 times mathbb{Z}_2) + T$ symmetry group, where $mathbb{Z}_2$ is some $pi$-rotation in the spin internal space and $T$ is time-reversal. One of these phases is equivalent to the usual spin-1 Haldane phase, while the other three are different but also typical of spin-1 systems. The model also exhibits an $SO(5)$-Haldane phase. Moreover, we also find that the transitions between the different effective spin-1 SPT phases are continuous, and can be described by a $c=2$ conformal field theory. At such transitions, indirect evidence suggests a possible effective field theory of four massless Majorana fermions. The results are obtained by approximating the ground state of the system in the thermodynamic limit using Matrix Product States via the infinite Time Evolving Block Decimation method, as well as by effective field theory considerations. Our findings show, for the first time, that different large effective spin-1 SPT phases separated by continuous quantum phase transitions can be stabilized in a simple quantum spin chain.
