Source author record

Frank Pollmann

Frank Pollmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

84works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

84 published item(s)

preprint2023arXiv

Hilbert Space Fragmentation in Open Quantum Systems

We investigate the phenomenon of Hilbert space fragmentation (HSF) in open quantum systems and find that it can stabilize highly entangled steady states. For concreteness, we consider the Temperley-Lieb model, which exhibits quantum HSF in an entangled basis, and investigate the Lindblad dynamics under two different couplings. First, we couple the system to a dephasing bath that reduces quantum fragmentation to a classical one with the resulting stationary state being separable. We observe that despite vanishing quantum correlations, classical correlations develop due to fluctuations of the remaining conserved quantities, which we show can be captured by a classical stochastic circuit evolution. Second, we use a coupling that preserves the quantum fragmentation structure. We derive a general expression for the steady state, which has a strong coherent memory of the initial state due to the extensive number of non-commuting conserved quantities. We show that it is highly entangled as quantified by the logarithmic negativity.

preprint2022arXiv

Crossing a topological phase transition with a quantum computer

Quantum computers promise to perform computations beyond the reach of modern computers with profound implications for scientific research. Due to remarkable technological advances, small scale devices are now becoming available for use. One of the most apparent applications for such a device is the study of complex many-body quantum systems, where classical computers are unable to deal with the generic exponential complexity of quantum states. Even zero-temperature equilibrium phases of matter and the transitions between them have yet to be fully classified, with topologically protected phases presenting major difficulties. We construct and measure a continuously parametrized family of states crossing a symmetry protected topological phase transition on the IBM Q quantum computers. We present two complementary methods for measuring string order parameters that reveal the transition, and additionally analyse the effects of noise in the device using simple error models. The simulation that we perform is easily scalable and is a practical demonstration of the utility of near-term quantum computers for the study of quantum phases of matter and their transitions.

preprint2022arXiv

Entanglement entropy and negativity in the Ising model with defects

Defects in two-dimensional conformal field theories (CFTs) contain signatures of their characteristics. In this work, we compute the entanglement entropy (EE) and the entanglement negativity (EN) of subsystems in the presence of energy and duality defects in the Ising CFT using the density matrix renormalization group (DMRG) technique. We show that the EE for the duality defect exhibits fundamentally different characteristics compared to the energy defect due to the existence of localized and delocalized zero energy modes. Of special interest is the nontrivial `finite-size correction' in the EE obtained recently using free fermion computations. These corrections arise when the subsystem size is appreciable compared to the total system size and lead to a deviation from the usual logarithmic scaling characteristic of one-dimensional quantum-critical systems. Using matrix product states with open and infinite boundary conditions, we numerically demonstrate the disappearance of the zero mode contribution for finite subsystem sizes in the thermodynamic limit. Our results provide further support to the recent free fermion computations, but clearly contradict earlier analytical field theory calculations based on twisted torus partition functions. Subsequently, we compute the logarithm of the EN (log-EN) between two disjoint subsystems separated by a defect. We show that the log-EN scales logarithmically with the separation of the subsystems. However, the coefficient of this logarithmic scaling yields a continuously-varying effective central charge that is different from that obtained from analogous computations of the EE. The defects leave their fingerprints in the subleading term of the scaling of the log-EN. Furthermore, the log-EN receives similar `finite size corrections' like the EE which leads to deviations from its characteristic logarithmic scaling.

preprint2022arXiv

Finite-depth scaling of infinite quantum circuits for quantum critical points

The scaling of the entanglement entropy at a quantum critical point allows us to extract universal properties of the state, e.g., the central charge of a conformal field theory. With the rapid improvement of noisy intermediate-scale quantum (NISQ) devices, these quantum computers present themselves as a powerful tool to study critical many-body systems. We use finite-depth quantum circuits suitable for NISQ devices as a variational ansatz to represent ground states of critical, infinite systems. We find universal finite-depth scaling relations for these circuits and verify them numerically at two different critical points, i.e., the critical Ising model with an additional symmetry-preserving term and the critical XXZ model.

preprint2022arXiv

Interaction-driven plateau transition between integer and fractional Chern Insulators

We present numerical evidence of an interaction-driven quantum Hall plateau transition between a $|C|>1$ Chern Insulator (CI) and a $ν= 1/3$ Laughlin state in the Harper-Hofstadter model. We study the model at flux densities $p/q$, where the lowest Landau level (LLL) manifold comprises $p$ magnetic sub-bands. For weak interactions, the model realises integer CIs corresponding to filled sub-bands, while strongly interacting candidate states include fractional quantum Hall (FQH) states at LLL filling fractions $ν=r/t$. These phases may compete at the same particle density when $p=t$. As a concrete example, we numerically explore the physics at flux density $n_ϕ = 3/11$, where we show evidence that a direct transition occurs between a CI and a $ν= 1/3$ Laughlin state, which we characterise in terms of its critical, topological and entanglement properties. We also show that strong interactions generically stabilise a $ν= 1/3$ Laughlin state even when the LLL is split into multiple bands, and introduce a powerful methodology to extract its topological entanglement entropy by exploiting the scaling of magnetic length with $n_ϕ$.

preprint2022arXiv

Topological phases in the dynamics of the simple exclusion process

We study the dynamical large deviations of the classical stochastic symmetric simple exclusion process (SSEP) by means of numerical matrix product states. We show that for half-filling, long-time trajectories with a large enough imbalance between the number hops in even and odd bonds of the lattice belong to distinct symmetry protected topological (SPT) phases. Using tensor network techniques, we obtain the large deviation (LD) phase diagram in terms of counting fields conjugate to the dynamical activity and the total hop imbalance. We show the existence of high activity trivial and non-trivial SPT phases (classified according to string-order parameters) separated by either a critical phase or a critical point. Using the leading eigenstate of the tilted generator, obtained from infinite-system density matrix renormalisation group (DMRG) simulations, we construct a near-optimal dynamics for sampling the LDs, and show that the SPT phases manifest at the level of rare stochastic trajectories. We also show how to extend these results to other filling fractions, and discuss generalizations to asymmetric SEPs.

preprint2021arXiv

Distinguishing an Anderson Insulator from a Many-Body Localized phase through space-time snapshots with Neural Networks

Distinguishing the dynamics of an Anderson insulator from a Many-Body Localized (MBL) phase is an experimentally challenging task. In this work, we propose a method based on machine learning techniques to analyze experimental snapshot data to separate the two phases. We show how to train $3D$ convolutional neural networks (CNNs) using space-time Fock-state snapshots, allowing us to obtain dynamic information about the system. We benchmark our method on a paradigmatic model showing MBL ($t-V$ model with quenched disorder), where we obtain a classification accuracy of $\approx 80 \%$ between an Anderson insulator and an MBL phase. We underline the importance of providing temporal information to the CNNs and we show that CNNs learn the crucial difference between an Anderson localized and an MBL phase, namely the difference in the propagation of quantum correlations. Particularly, we show that the misclassified MBL samples are characterized by an unusually slow propagation of quantum correlations, and thus the CNNs label them wrongly as Anderson localized. Finally, we apply our method to the case with quasi-periodic potential, known as the Aubry-André model (AA model). We find that the CNNs have more difficulties in separating the two phases. We show that these difficulties are due to the fact that the MBL phase of the AA model is characterized by a slower information propagation for numerically accessible system sizes.

preprint2021arXiv

Fractonic critical point proximate to a higher-order topological insulator: How does UV blend with IR?

We propose an unconventional topological quantum phase transition connecting a higher-order topological insulator (HOTI) and a featureless Mott insulator sharing the same symmetry patterns. We construct an effective theory description of the quantum critical point (QCP) by combining a bosonization approach and the coupled-stripe construction of 1D critical spin ladders. The phase transition theory is characterized by a critical dipole liquid theory with subsystem $U(1)$ symmetry whose low energy modes contain a Bose surface along the $k_x,k_y$ axis. Such a quantum critical point manifests fracton dynamics and the breakdown of the area law entanglement entropy due to the existence of a Bose surface. We numerically confirm our findings by measuring the entanglement entropy, topological rank-2 Berry phase, and the static structure factor throughout the topological transition and compare it with our previous approach obtained from the percolation picture. A significant new element of our phase transition theory is that the infrared~(IR) effective theory is controlled by short wave-length fluctuations with peculiar UV-IR mixing.

preprint2021arXiv

Gapless topological phases and symmetry-enriched quantum criticality

We introduce topological invariants for gapless systems and study the associated boundary phenomena. More generally, the symmetry properties of the low-energy conformal field theory (CFT) provide discrete invariants, establishing the notion of symmetry-enriched quantum criticality. The charges of nonlocal scaling operators, or more generally of symmetry defects, are topological and imply the presence of localized edge modes. We primarily focus on the $1+1d$ case where the edge has a topological degeneracy, whose finite-size splitting can be exponential or algebraic in system size depending on the involvement of additional gapped sectors. An example of the former is given by tuning the spin-1 Heisenberg chain to a symmetry-breaking Ising phase. An example of the latter arises between the gapped Ising and cluster phases: this symmetry-enriched Ising CFT has an edge mode with finite-size splitting $\sim 1/L^{14}$. In addition to such new cases, our formalism unifies various examples previously studied in the literature. Similar to gapped symmetry-protected topological phases, a given CFT can split into several distinct symmetry-enriched CFTs. This raises the question of classification, to which we give a partial answer -- including a complete characterization of symmetry-enriched $1+1d$ Ising CFTs. Non-trivial topological invariants can also be constructed in higher dimensions, which we illustrate for a symmetry-enriched $2+1d$ CFT without gapped sectors.

preprint2021arXiv

Identifying Correlation Clusters in Many-Body Localized Systems

We introduce techniques for analysing the structure of quantum states of many-body localized (MBL) spin chains by identifying correlation clusters from pairwise correlations. These techniques proceed by interpreting pairwise correlations in the state as a weighted graph, which we analyse using an established graph theoretic clustering algorithm. We validate our approach by studying the eigenstates of a disordered XXZ spin chain across the MBL to ergodic transition, as well as the non-equilibrium dyanmics in the MBL phase following a global quantum quench. We successfully reproduce theoretical predictions about the MBL transition obtained from renormalization group schemes. Furthermore, we identify a clear signature of many-body dynamics analogous to the logarithmic growth of entanglement. The techniques that we introduce are computationally inexpensive and in combination with matrix product state methods allow for the study of large scale localized systems. Moreover, the correlation functions we use are directly accessible in a range of experimental settings including cold atoms.

preprint2021arXiv

Operator backflow and the classical simulation of quantum transport

Tensor product states have proved extremely powerful for simulating the area-law entangled states of many-body systems, such as the ground states of gapped Hamiltonians in one dimension. The applicability of such methods to the \emph{dynamics} of many-body systems is less clear: the memory required grows exponentially in time in most cases, quickly becoming unmanageable. New methods reduce the memory required by selectively discarding/dissipating parts of the many-body wavefunction which are expected to have little effect on the hydrodynamic observables typically of interest: for example, some methods discard fine-grained correlations associated with $n$-point functions, with $n$ exceeding some cutoff $\ell_*$. In this work, we present a theory for the sizes of `backflow corrections', i.e., systematic errors due to discarding this fine-grained information. In particular, we focus on their effect on transport coefficients. Our results suggest that backflow corrections are exponentially suppressed in the size of the cutoff $\ell_*$. Moreover, the backflow errors themselves have a hydrodynamical expansion, which we elucidate. We test our predictions against numerical simulations run on random unitary circuits and ergodic spin-chains. These results lead to the conjecture that transport coefficients in ergodic diffusive systems can be captured to a given precision $ε$ with an amount of memory scaling as $\exp[\mathcal{O}(\log(ε)^2)]$, significantly better than the naive estimate of memory $\exp[\mathcal{O}(\mathrm{poly}(ε^{-1}))]$ required by more brute-force methods.

preprint2021arXiv

Quantum phases of two-dimensional $\mathbb{Z}_2$ gauge theory coupled to single-component fermion matter

We investigate the rich quantum phase diagram of Wegner's theory of discrete Ising gauge fields interacting with $U(1)$ symmetric single-component fermion matter hopping on a two-dimensional square lattice. In particular limits the model reduces to (i) pure $\mathbb{Z}_2$ even and odd gauge theories, (ii) free fermions in a static background of deconfined $\mathbb{Z}_2$ gauge fields, (iii) the kinetic Rokhsar-Kivelson quantum dimer model at a generic dimer filling. We develop a local transformation that maps the lattice gauge theory onto a model of $\mathbb{Z}_2$ gauge-invariant spin $1/2$ degrees of freedom. Using the mapping, we perform numerical density matrix renormalization group calculations that corroborate our understanding of the limits identified above. Moreover, in the absence of the magnetic plaquette term, we reveal signatures of topologically ordered Dirac semimetal and staggered Mott insulator phases at half-filling. At strong coupling, the lattice gauge theory displays fracton phenomenology with isolated fermions being completely frozen and dimers exhibiting restricted mobility. In that limit, we predict that in the ground state dimers form compact clusters, whose hopping is suppressed exponentially in their size. We determine the band structure of the smallest clusters numerically using exact diagonalization. The rich phenomenology discussed in this paper can be probed in analog and digital quantum simulators of discrete gauge theories and in Kitaev spin-orbital liquids.

preprint2021arXiv

Quotient symmetry protected topological phenomena

Topological phenomena are commonly studied in phases of matter which are separated from a trivial phase by an unavoidable quantum phase transition. This can be overly restrictive, leaving out scenarios of practical relevance -- similar to the distinction between liquid water and vapor. Indeed, we show that topological phenomena can be stable over a large part of parameter space even when the bulk is strictly speaking in a trivial phase of matter. In particular, we focus on symmetry-protected topological phases which can be trivialized by extending the symmetry group. The topological Haldane phase in spin chains serves as a paradigmatic example where the $SO(3)$ symmetry is extended to $SU(2)$ by tuning away from the Mott limit. Although the Haldane phase is then adiabatically connected to a product state, we show that characteristic phenomena -- edge modes, entanglement degeneracies and bulk phase transitions -- remain parametrically stable. This stability is due to a separation of energy scales, characterized by quantized invariants which are well-defined when a subgroup of the symmetry only acts on high-energy degrees of freedom. The low-energy symmetry group is a quotient group whose emergent anomalies stabilize edge modes and unnecessary criticality, which can occur in any dimension.

preprint2021arXiv

Realising the Symmetry-Protected Haldane Phase in Fermi-Hubbard Ladders

Topology in quantum many-body systems has profoundly changed our understanding of quantum phases of matter. The paradigmatic model that has played an instrumental role in elucidating these effects is the antiferromagnetic spin-1 Haldane chain. Its ground state is a disordered state, with symmetry-protected fourfold-degenerate edge states due to fractional spin excitations. In the bulk, it is characterised by vanishing two-point spin correlations, gapped excitations, and a characteristic non-local order parameter. More recently it was understood that the Haldane chain forms a specific example of a more general classification scheme of symmetry protected topological (SPT) phases of matter that is based on ideas connecting to quantum information and entanglement. Here, we realise such a topological Haldane phase with Fermi-Hubbard ladders in an ultracold-atom quantum simulator. We directly reveal both edge and bulk properties of the system through the use of single-site and particle-resolved measurements as well as non-local correlation functions. Continuously changing the Hubbard interaction strength of the system allows us to investigate the robustness of the phase to charge (density) fluctuations far from the regime of the Heisenberg model employing a novel correlator.

preprint2021arXiv

Scaling of neural-network quantum states for time evolution

Simulating quantum many-body dynamics on classical computers is a challenging problem due to the exponential growth of the Hilbert space. Artificial neural networks have recently been introduced as a new tool to approximate quantum-many body states. We benchmark the variational power of the restricted Boltzmann machine quantum states and different shallow and deep neural autoregressive quantum states to simulate global quench dynamics of a non-integrable quantum Ising chain. We find that the number of parameters required to represent the quantum state at a given accuracy increases exponentially in time. The growth rate is only slightly affected by the network architecture over a wide range of different design choices: shallow and deep networks, small and large filter sizes, dilated and normal convolutions, with and without shortcut connections.

preprint2020arXiv

Dissipation-assisted operator evolution method for capturing hydrodynamic transport

We introduce the dissipation-assisted operator evolution (DAOE) method for calculating transport properties of strongly interacting lattice systems in the high temperature regime. DAOE is based on evolving observables in the Heisenberg picture, and applying an artificial dissipation that reduces the weight on non-local operators. We represent the observable as a matrix product operator, and show that the dissipation leads to a decay of operator entanglement, allowing us to capture the dynamics to long times. We test this scheme by calculating spin and energy diffusion constants in a variety of physical models. By gradually weakening the dissipation, we are able to consistently extrapolate our results to the case of zero dissipation, thus estimating the physical diffusion constant with high precision.

preprint2020arXiv

Entanglement dynamics of a many-body localized system coupled to a bath

The combination of strong disorder and interactions in closed quantum systems can lead to many-body localization (MBL). However this quantum phase is not stable when the system is coupled to a thermal environment. We investigate how MBL is destroyed in systems that are weakly coupled to a dephasive Markovian environment by focusing on their entanglement dynamics. We numerically study the third Rényi negativity $R_3$, a recently proposed entanglement proxy based on the negativity that captures the unbounded logarithmic growth in the closed case and that can be computed efficiently with tensor networks. We also show that the decay of $R_3$ follows a stretched exponential law, similarly to the imbalance, with however a smaller stretching exponent.

preprint2020arXiv

Entanglement Hamiltonian of the 1+1-dimensional free, compactified boson conformal field theory

Entanglement or modular Hamiltonians play a crucial role in the investigation of correlations in quantum field theories. In particular, in 1+1 space-time dimensions, the spectra of entanglement Hamiltonians of conformal field theories (CFTs) for certain geometries are related to the spectra of the physical Hamiltonians of corresponding boundary CFTs. As a result, conformal invariance allows exact computation of the spectra of the entanglement Hamiltonians for these models. In this work, we perform this computation of the spectrum of the entanglement Hamiltonian for the free compactified boson CFT over a finite spatial interval. We compare the analytical results obtained for the continuum theory with numerical simulations of a lattice-regularized model for the CFT using density matrix renormalization group technique. To that end, we use a lattice regularization provided by superconducting quantum electronic circuits, built out of Josephson junctions and capacitors. Up to non-universal effects arising due to the lattice regularization, the numerical results are compatible with the predictions of the exact computations.

preprint2020arXiv

Ergodicity-breaking arising from Hilbert space fragmentation in dipole-conserving Hamiltonians

We show that the combination of charge and dipole conservation---characteristic of fracton systems---leads to an extensive fragmentation of the Hilbert space, which in turn can lead to a breakdown of thermalization. As a concrete example, we investigate the out-of-equilibrium dynamics of one-dimensional spin-1 models that conserve charge (total $S^z$) and its associated dipole moment. First, we consider a minimal model including only three-site terms and find that the infinite temperature auto-correlation saturates to a finite value---showcasing non-thermal behavior. The absence of thermalization is identified as a consequence of the strong fragmentation of the Hilbert space into exponentially many invariant subspaces in the local $S^z$ basis, arising from the interplay of dipole conservation and local interactions. Second, we extend the model by including four-site terms and find that this perturbation leads to a weak fragmentation: the system still has exponentially many invariant subspaces, but they are no longer sufficient to avoid thermalization for typical initial states. More generally, for any finite range of interactions, the system still exhibits non-thermal eigenstates appearing throughout the entire spectrum. We compare our results to charge and dipole moment conserving random unitary circuit models for which we reach identical conclusions.

preprint2020arXiv

Fracton Critical Point in Higher-Order Topological Phase Transition

The theory of quantum phase transitions separating different phases with distinct symmetry patterns at zero temperature is one of the foundations of modern quantum many-body physics. In this paper we demonstrate that the existence of a 2D topological phase transition between a higher-order topological insulator (HOTI) and a trivial Mott insulator with the same symmetry eludes this paradigm. We present a theory of this quantum critical point (QCP) driven by the fluctuations and percolation of the domain walls between a HOTI and a trivial Mott insulator region. Due to the spinon zero modes that decorate the rough corners of the domain walls, the fluctuations of the phase boundaries trigger a spinon-dipole hopping term with fracton dynamics. Hence we find the QCP is characterized by a critical dipole liquid theory with subsystem $U(1)$ symmetry and the breakdown of the area law entanglement entropy which exhibits a logarithmic enhancement: $L \ln(L)$. Using the density matrix renormalization group (DMRG) method, we analyze the dipole stiffness together with structure factor at the QCP which provide strong evidence of a critical dipole liquid with a Bose surface. These numerical signatures further support the fracton dynamics of the QCP, and suggest a new paradigm for 2D quantum criticality proximate to a topological phase.

preprint2020arXiv

Higher-Order Entanglement and Many-Body Invariants for Higher-Order Topological Phases

We discuss how strongly interacting higher-order symmetry protected topological (HOSPT) phases can be characterized from the entanglement perspective: First, we introduce a topological many-body invariant which reveals the non-commutative algebra between flux operator and $C_n$ rotations. We argue that this invariant denotes the angular momentum carried by the instanton which is closely related to the discrete Wen-Zee response and fractional corner charge. Second, we define a new entanglement property, dubbed `higher-order entanglement', to scrutinize and differentiate various higher-order topological phases from a hierarchical sequence of the entanglement structure. We support our claims by numerically studying a super-lattice Bose-Hubbard model that exhibits different HOSPT phases.

preprint2020arXiv

Isometric Tensor Network States in Two Dimensions

Tensor network states (TNS) are a promising but numerically challenging tool for simulating two-dimensional (2D) quantum many-body problems. We introduce an isometric restriction of the TNS ansatz that allows for highly efficient contraction of the network. We consider two concrete applications using this ansatz. First, we show that a matrix-product state representation of a 2D quantum state can be iteratively transformed into an isometric 2D TNS. Second, we introduce a 2D version of the time-evolving block decimation algorithm (TEBD$^2$) for approximating the ground state of a Hamiltonian as an isometric TNS, which we demonstrate for the 2D transverse field Ising model.

preprint2020arXiv

Observing Quasiparticles through the Entanglement Lens

The low energy physics of interacting quantum systems is typically understood through the identification of the relevant quasiparticles or low energy excitations and their quantum numbers. We present a quantum information framework that goes beyond this to examine the nature of the entanglement in the corresponding quantum states. We argue that the salient features of the quasiparticles, including their quantum numbers, locality and fractionalization are reflected in the entanglement spectrum and in the mutual information. We illustrate these ideas in the specific context of the $d=1$ transverse field Ising model with an integrability breaking perturbation.

preprint2020arXiv

Parton theory of ARPES spectra in anti-ferromagnetic Mott insulators

Angle-resolved photoemission spectroscopy (ARPES) has revealed peculiar properties of mobile dopants in correlated anti-ferromagnets (AFMs). But describing them theoretically, even in simplified toy models, remains a challenge. Here we study ARPES spectra of a single mobile hole in the $t-J$ model. Recent progress in the microscopic description of mobile dopants allows us to use a geometric decoupling of spin and charge fluctuations at strong couplings, from which we conjecture a one-to-one relation of the one-dopant spectral function and the spectrum of a constituting spinon in the \emph{undoped} parent AFM. We thoroughly test this hypothesis for a single hole doped into a 2D Heisenberg AFM by comparing our semi-analytical predictions to previous quantum Monte Carlo results and our large-scale time-dependent matrix product state (td-MPS) calculations of the spectral function. Our conclusion is supported by a microscopic trial wavefuntion describing spinon-chargon bound states, which captures the momentum and $t/J$ dependence of the quasiparticle residue. Our conjecture suggests that ARPES measurements in the pseudogap phase of cuprates can directly reveal the Dirac-fermion nature of the constituting spinons. Specifically, we demonstrate that our trial wavefunction provides a microscopic explanation for the sudden drop of spectral weight around the nodal point associated with the formation of Fermi arcs, assuming that additional frustration suppresses long-range AFM ordering. We benchmark our results by studying the cross-over from two to one dimension, where spinons and chargons are confined and deconfined respectively.

preprint2020arXiv

Quantum phases of a one-dimensional Majorana-Bose-Hubbard model

Majorana zero modes (MZM-s) occurring at the edges of a 1D, p-wave, spinless superconductor, in absence of fluctuations of the phase of the superconducting order parameter, are quintessential examples of topologically-protected zero-energy modes occurring at the edges of 1D symmetry-protected topological phases. In this work, we numerically investigate the fate of the topological phase in the presence of phase-fluctuations using the density matrix renormalization group (DMRG) technique. To that end, we consider a one-dimensional array of MZM-s on mesoscopic superconducting islands at zero temperature. Cooper-pair and MZM-assisted single-electron tunneling, together with finite charging energy of the mesoscopic islands, give rise to a rich phase-diagram of this model. We show that the system can be in either a Mott-insulating phase, a Luttinger liquid (LL) phase of Cooper-pairs or a second gapless phase. In contrast to the LL of Cooper-pairs, this second phase is characterized by nonlocal string correlation functions which decay algebraically due to gapless charge-$e$ excitations. The three phases are separated from each other by phase-transitions of either Kosterlitz-Thouless or Ising type. Using a Jordan-Wigner transformation, we map the system to a generalized Bose-Hubbard model with two types of hopping and use DMRG to analyze the different phases and the phase-transitions.

preprint2020arXiv

Statistical localization: from strong fragmentation to strong edge modes

Certain disorder-free Hamiltonians can be non-ergodic due to a \emph{strong fragmentation} of the Hilbert space into disconnected sectors. Here, we characterize such systems by introducing the notion of `statistically localized integrals of motion' (SLIOM), whose eigenvalues label the connected components of the Hilbert space. SLIOMs are not spatially localized in the operator sense, but appear localized to sub-extensive regions when their expectation value is taken in typical states with a finite density of particles. We illustrate this general concept on several Hamiltonians, both with and without dipole conservation. Furthermore, we demonstrate that there exist perturbations which destroy these integrals of motion in the bulk of the system, while keeping them on the boundary. This results in statistically localized \emph{strong zero modes}, leading to infinitely long-lived edge magnetizations along with a thermalizing bulk, constituting the first example of such strong edge modes in a non-integrable model. We also show that in a particular example, these edge modes lead to the appearance of topological string order in a certain subset of highly excited eigenstates. Some of our suggested models can be realized in Rydberg quantum simulators.

preprint2019arXiv

Dynamics of strongly interacting systems: From Fock-space fragmentation to Many-Body Localization

We study the $t{-}V$ disordered spinless fermionic chain in the strong coupling regime, $t/V\rightarrow 0$. Strong interactions highly hinder the dynamics of the model, fragmenting its Hilbert space into exponentially many blocks in system size. Macroscopically, these blocks can be characterized by the number of new degrees of freedom, which we refer to as movers. We focus on two limiting cases: Blocks with only one mover and the ones with a finite density of movers. The former many-particle block can be exactly mapped to a single-particle Anderson model with correlated disorder in one dimension. As a result, these eigenstates are always localized for any finite amount of disorder. The blocks with a finite density of movers, on the other side, show an MBL transition that is tuned by the disorder strength. Moreover, we provide numerical evidence that its ergodic phase is diffusive at weak disorder. Approaching the MBL transition, we observe sub-diffusive dynamics at finite time scales and find indications that this might be only a transient behavior before crossing over to diffusion.

preprint2019arXiv

Fractional corner charges in a 2D super-lattice Bose-Hubbard model

We study a two dimensional super-lattice Bose-Hubbard model with alternating hoppings in the limit of strong on-site interactions. We evaluate the phase diagram of the model around half-filling using the density matrix renormalization group method and find two gapped phases separated by a gapless superfluid region. We demonstrate that the gapped states realize two distinct higher order symmetry protected topological phases, which are protected by a combination of charge conservation and $C_4$ lattice symmetry. The phases are distinguished in terms of a quantized fractional corner charge and a many-body topological invariant that is robust against arbitrary, symmetry preserving edge manipulations. We support our claims by numerically studying the full counting statistics of the corner charge, finding a sharp distribution peaked around the quantized values. These results are experimentally observable in ultracold atomic settings using state of the art quantum gas microscopy.

preprint2019arXiv

Isometric Tensor Network representation of string-net liquids

Recently, a class of tensor networks called isometric tensor network states (isoTNS) was proposed which generalizes the canonical form of matrix product states to tensor networks in higher dimensions. While this ansatz allows for efficient numerical computations, it remained unclear which phases admit an isoTNS representation. In this work, we show that two-dimensional string-net liquids, which represent a wide variety of topological phases including discrete gauge theories, admit an exact isoTNS representation. We further show that the isometric form can be preserved after applying a finite depth local quantum circuit. Taken together, these results show that long-range entanglement by itself is not an obstruction to isoTNS representation and suggest that all two-dimensional gapped phases with gappable edges admit an isoTNS representation.

preprint2019arXiv

Many-body topological invariants from randomized measurements

The classification of symmetry-protected topological (SPT) phases in one dimension has been recently achieved, and had a fundamental impact in our understanding of quantum phases in condensed matter physics. In this framework, SPT phases can be identified by many-body topological invariants, which are quantized non-local correlators for the many-body wavefunction. While SPT phases can now be realized in interacting synthethic quantum systems, the direct measurement of quantized many-body topological invariants has remained so far elusive. Here, we propose measurement protocols for many-body topological invariants for all types of protecting symmetries of one-dimensional interacting bosonic systems. Our approach relies on randomized measurements implemented with local random unitaries, and can be applied to any spin system with single-site addressability and readout. Our scheme thus provides a versatile toolbox to experimentally classify interacting SPT phases.

preprint2016arXiv

Bridging coupled wires and lattice Hamiltonian for two-component bosonic quantum Hall states

We investigate a model of hard-core bosons with correlated hopping on the honeycomb lattice in an external magnetic field by means of a coupled-wire approach. It has been numerically shown that this model exhibits at half filling the bosonic integer quantum Hall (BIQH) state, which is a symmetry-protected topological phase protected by the $U(1)$ particle conservation [Y.-C. He et al., Phys. Rev. Lett. 115, 116803 (2015)]. By combining the bosonization approach and a coupled-wire construction, we analytically confirm this finding and show that it even holds in the strongly anisotropic (quasi-one-dimensional) limit. We discuss the stability of the BIQH phase against tunneling that break the separate particle conservations on different sublattices down to a global particle conservation. We further argue that a phase transition between two different BIQH phases is in a deconfined quantum critical point described by two copies of the (2+1)-dimensional $O(4)$ nonlinear sigma model with the topological $θ$ term at $θ=π$. Finally we predict a possible fractional quantum Hall state, the Halperin (221) state, at 1/6 filling.

preprint2016arXiv

Density matrix renormalization group on a cylinder in mixed real and momentum space

We develop a variant of the density matrix renormalization group (DMRG) algorithm for two-dimensional cylinders that uses a real space representation along the cylinder and a momentum space representation in the perpendicular direction. The mixed representation allows us to use the momentum around the circumference as a conserved quantity in the DMRG algorithm. Compared with the traditional purely real-space approach, we find a significant speedup in computation time and a considerable reduction in memory usage. Applying the method to the interacting fermionic Hofstadter model, we demonstrate a reduction in computation time by over 20-fold, in addition to a sixfold memory reduction.

preprint2016arXiv

Domain-wall melting as a probe of many-body localization

Motivated by a recent optical-lattice experiment by Choi et al.[Science 352, 1547 (2016)], we discuss how domain-wall melting can be used to investigate many-body localization. First, by considering noninteracting fermion models, we demonstrate that experimentally accessible measures are sensitive to localization and can thus be used to detect the delocalization-localization transition, including divergences of characteristic length scales. Second, using extensive time-dependent density matrix renormalization group simulations, we study fermions with repulsive interactions on a chain and a two-leg ladder. The extracted critical disorder strengths agree well with the ones found in existing literature.

preprint2016arXiv

How periodic driving heats a disordered quantum spin chain

We study the energy absorption in real time of a disordered quantum spin chain subjected to coherent monochromatic periodic driving. We determine characteristic fingerprints of the well-known ergodic (Floquet-ETH for slow driving/weak disorder) and many-body localized (Floquet-MBL for fast driving/strong disorder) phases. In addition, we identify an intermediate regime, where the energy density of the system -- unlike the entanglement entropy a local and bounded observable -- grows logarithmically slowly over a very large time window.

preprint2016arXiv

Kagome chiral spin liquid as a gauged U(1) symmetry protected topological phase

While the existence of a chiral spin liquid (CSL) on a class of spin-1/2 kagome antiferromagnets is by now well-established numerically, a controlled theoretical path from the lattice model leading to a low energy topological field theory is still lacking. This we provide via an explicit construction, starting from reformulating a microscopic model for a CSL as a lattice gauge theory, and deriving the low-energy form of its continuum limit. A crucial ingredient is the realisation that the bosonic spinons of the gauge theory exhibit a $U(1)$ symmetry-protected topological (SPT) phase, which upon promoting its $U(1)$ global symmetry to a local gauge structure ("gauging") yields the CSL. We suggest that such an explicit lattice-based construction involving gauging of an SPT phase can be applied more generally to understand topological spin liquids.

preprint2016arXiv

Momentum-space entanglement and Loschmidt echo in Luttinger liquids after a quantum quench

Luttinger liquids (LLs) arise by coupling left- and right-moving particles through interactions in one dimension. This most natural partitioning of LLs is investigated by the momentum-space entanglement after a quantum quench using analytical and numerical methods. We show that the momentum-space entanglement spectrum of a LL possesses many universal features both in equilibrium and after a quantum quench. The largest entanglement eigenvalue is identical to the Loschmidt echo, i.e. the overlap of the disentangled and final wavefunctions of the system. The second largest eigenvalue is the overlap of the first excited state of the disentangled system with zero total momentum and the final wavefunction. The entanglement gap is universal both in equilibrium and after a quantum quench. The momentum-space entanglement entropy is always extensive and saturates fast to a time independent value after the quench, in sharp contrast to a spatial bipartitioning.

preprint2016arXiv

Obtaining highly-excited eigenstates of many-body localized Hamiltonians by the density matrix renormalization group

The eigenstates of many-body localized (MBL) Hamiltonians exhibit low entanglement. We adapt the highly successful density-matrix renormalization group method, which is usually used to find modestly entangled ground states of local Hamiltonians, to find individual highly excited eigenstates of many body localized Hamiltonians. The adaptation builds on the distinctive spatial structure of such eigenstates. We benchmark our method against the well studied random field Heisenberg model in one dimension. At moderate to large disorder, we find that the method successfully obtains excited eigenstates with high accuracy, thereby enabling a study of MBL systems at much larger system sizes than those accessible to exact-diagonalization methods.

preprint2016arXiv

Typicality approach to the optical conductivity in thermal and many-body localized phases

We study the frequency dependence of the optical conductivity $\text{Re} \, σ(ω)$ of the Heisenberg spin-$1/2$ chain in the thermal and near the transition to the many-body localized phase induced by the strength of a random $z$-directed magnetic field. Using the method of dynamical quantum typicality, we calculate the real-time dynamics of the spin-current autocorrelation function and obtain the Fourier transform $\text{Re} \, σ(ω)$ for system sizes much larger than accessible to standard exact-diagonalization approaches. We find that the low-frequency behavior of $\text{Re} \, σ(ω)$ is well described by $\text{Re} \, σ(ω) \approx σ_\text{dc} + a \, |ω|^α$, with $α\approx 1$ in a wide range within the thermal phase and close to the transition. We particularly detail the decrease of $σ_\text{dc}$ in the thermal phase as a function of increasing disorder for strong exchange anisotropies. We further find that the temperature dependence of $σ_\text{dc}$ is consistent with the existence of a mobility edge.

preprint2015arXiv

Absence of orthogonality catastrophe after a spatially inhomogeneous interaction quench in Luttinger liquids

We investigate the Loschmidt echo, the overlap of the initial and final wavefunctions of Luttinger liquids after a spatially inhomogeneous interaction quench. In studying the Luttinger model, we obtain an analytic solution of the bosonic Bogoliubov-de Gennes equations after quenching the interactions within a finite spatial region. As opposed to the power law temporal decay following a potential quench, the interaction quench in the Luttinger model leads to a finite, hardly time dependent overlap, therefore no orthogonality catastrophe occurs. The steady state value of the Loschmidt echo after a sudden inhomogeneous quench is the square of the respective adiabatic overlaps. Our results are checked and validated numerically on the XXZ Heisenberg chain.

preprint2015arXiv

Bosonic Integer Quantum Hall effect in an interacting lattice model

We study a bosonic model with correlated hopping on a honeycomb lattice, and show that its ground state is a bosonic integer quantum Hall (BIQH) phase, a prominent example of a symmetry protected topological (SPT) phase. By using the infinite density matrix renormalization group method, we establish the existence of the BIQH phase by providing clear numerical evidence: (i) a quantized Hall conductance with $|σ_{xy}|= 2$ (ii) two counter propagating gapless edge modes. Our simple model is an example of a novel class of systems that can stabilize SPT phases protected by a continuous symmetry on lattices and opens up new possibilities for the experimental realization of these exotic phases.

preprint2015arXiv

Distinct trivial phases protected by a point-group symmetry in quantum spin chains

The ground state of the $S=1$ antiferromagnetic Heisenberg chain belongs to the Haldane phase -- a well known example of symmetry-protected topological phase. A staggered field applied to the $S=1$ antiferromagnetic chain breaks all the symmetries that protect the Haldane phase as a topological phase, reducing it to a trivial phase. That is, the Haldane phase is then connected adiabatically to an antiferromagnetic product state. Nevertheless, as long as the symmetry under site-centered inversion combined with a spin rotation is preserved, the phase is still distinct from another trivial phase. We demonstrate the existence of such distinct symmetry-protected trivial phases using a field-theoretical approach and numerical calculations. Furthermore, a general proof and a non-local order parameter are given in terms of an matrix-product state formulation.

preprint2015arXiv

Efficient variational diagonalization of fully many-body localized Hamiltonians

We introduce a unitary matrix-product operator (UMPO) based variational method that approximately finds all the eigenstates of fully many-body localized (fMBL) one-dimensional Hamiltonians. The computational cost of the variational optimization scales linearly with system size for a fixed bond dimension of the UMPO ansatz. We demonstrate the usefulness of our approach by considering the Heisenberg chain in a strongly disordered magnetic field for which we compare the approximation to exact diagonalization results.

preprint2015arXiv

Infinite density matrix renormalization group for multicomponent quantum Hall systems

While the simplest quantum Hall plateaus, such as the $ν= 1/3$ state in GaAs, can be conveniently analyzed by assuming only a single active Landau level participates, for many phases the spin, valley, bilayer, subband, or higher Landau level indices play an important role. These `multi-component' problems are difficult to study using exact diagonalization because each component increases the difficulty exponentially. An important example is the plateau at $ν= 5/2$, where scattering into higher Landau levels chooses between the competing non-Abelian Pfaffian and anti-Pfaffian states. We address the methodological issues required to apply the infinite density matrix renormalization group to quantum Hall systems with multiple components and long-range Coulomb interactions, greatly extending accessible system sizes. As an initial application we study the problem of Landau level mixing in the $ν= 5/2$ state. Within the approach to Landau level mixing used here, we find that at the Coulomb point the anti-Pfaffian is preferred over the Pfaffian state over a range of Landau level mixing up to the experimentally relevant values.

preprint2015arXiv

Interaction driven phases in the half-filled honeycomb lattice: an infinite density matrix renormalization group study

The emergence of the Haldane Chern insulator state due to strong short range repulsive interactions in the half-filled fermionic spinless honeycomb lattice model has been proposed and challenged with different methods and yet it still remains controversial. In this work we revisit the problem using the infinite density matrix renormalization group method and report numerical evidence supporting i) the absence of the Chern insulator state, ii) two previously unnoticed charge ordered phases and iii) the existence and stability of all the non-topological competing orders that were found previously within mean field. In addition, we discuss the nature of the corresponding phase transitions based on our numerical data. Our work establishes the phase diagram of the half-filled honeycomb lattice model tilting the balance towards the absence of a Chern insulator phase for this model.

preprint2015arXiv

Of chain-based order and quantum spin liquids in dipolar spin ice

Recent experiments on the spin-ice material Dy2Ti2O7 suggest that the Pauling "ice entropy", characteristic of its classical Coulombic spin-liquid state, may be lost at low temperatures [D. Pomaranski et al., Nature Phys. 9, 353 (2013)]. However, despite nearly two decades of intensive study, the nature of the equilibrium ground state of spin ice remains uncertain. Here we explore how long-range dipolar interactions D, short-range exchange interactions, and quantum fluctuations combine to determine the ground state of dipolar spin ice. We identify a new organisational principle, namely that ordered ground states are selected from a set of "chain states" in which dipolar interactions are exponentially screened. Using both quantum and classical Monte Carlo simulation, we establish phase diagrams as a function of quantum tunneling g, and temperature T, and find that only a very small g_c << D is needed to stabilize a quantum spin-liquid ground state. We discuss the implications of these results for Dy2Ti2O7.

preprint2015arXiv

Quantum quench in 2D using the variational Baeriswyl wavefunction

By combining the Baeriswyl wavefunction with equilibrium and time-dependent variational principles, we develop a non-equilibrium formalism to study quantum quenches for two dimensional spinless fermions with nearest-neighbour hopping and repulsion. The variational ground state energy and the short time dynamics agree convincingly with the results of numerically exact simulations. We find that depending on the initial and final interaction strength, the quenched system either exhibits undamped oscillations or relaxes to a time independent steady state. The time averaged expectation value of the CDW order parameter rises sharply when crossing from the steady state regime to the oscillating regime, indicating that the system, being non-integrable, shows signs of thermalization with an effective temperature above or below the equilibrium critical temperature, respectively.

preprint2015arXiv

Signatures of the many-body localization transition in the dynamics of entanglement and bipartite fluctuations

The many-body localization transition is a dynamical quantum phase transition between a localized and an extended phase. We study this transition in the XXZ model with disordered magnetic field and focus on the time evolution following a global quantum quench. While the dynamics of the bipartite entanglement and spin fluctuations are already known to provide insights into the nature of the many-body localized phases, we discuss the relevance of these quantities in the context of the localization transition. In particular, we observe that near the transition the long time limits of both quantities show behavior similar to divergent thermodynamic fluctuations.

preprint2015arXiv

Sudden expansion and domain-wall melting of strongly interacting bosons in two-dimensional optical lattices and on multileg ladders

We numerically investigate the expansion of clouds of hard-core bosons in the two-dimensional square lattice using a matrix-product-state based method. This nonequilibrium set-up is induced by quenching the trapping potential to zero and our work is specifically motivated by a recent experiment with interacting bosons in an optical lattice [Ronzheimer et al., Phys. Rev. Lett. 110, 205301 (2013)]. As the anisotropy of the amplitudes $J_x$ and $J_y$ for hopping in different spatial directions is varied from the one- to the two-dimensional case, we observe a crossover from a fast ballistic expansion in the one-dimensional limit $J_x \gg J_y$ to much slower dynamics in the isotropic two-dimensional limit $J_x=J_y$. We further study the dynamics on multi-leg ladders and long cylinders. For these geometries we compare the expansion of a cloud to the melting of a domain wall, which helps us to identify several different regimes of the expansion as a function of time. By studying the dependence of expansion velocities on both the anisotropy $J_y/J_x$ and the number of legs, we observe that the expansion on two-leg ladders, while similar to the two-dimensional case, is slower than on wider ladders. We provide a qualitative explanation for this observation based on an analysis of the rung spectrum.

preprint2015arXiv

Variational Monte Carlo simulations using tensor-product projected states

We propose an efficient numerical method, which combines the advantages of recently developed tensor-network based methods and standard trial wave functions, to study the ground state properties of quantum many-body systems. In this approach, we apply a projector in the form of a tensor-product operator to an input wave function, such as a Jastrow-type or Hartree-Fock wave function, and optimize the tensor elements via variational Monte Carlo. The entanglement already contained in the input wave function can considerably reduce the bond dimensions compared to the regular tensor-product state representation. In particular, this allows us to also represent states that do not obey the area law of entanglement entropy. In addition, for fermionic systems, the fermion sign structure can be encoded in the input wave function. We show that the optimized states provide good approximations of the ground-state energy and correlation functions in the cases of two-dimensional bosonic and fermonic systems.

preprint2014arXiv

Berry phase induced dimerization in one-dimensional quadrupolar systems

We investigate the effect of the Berry phase on quadrupoles that occur for example in the low-energy description of spin models. Specifically we study here the one-dimensional bilinear-biquadratic spin-one model. An open question for many years about this model is whether it has a non-dimerized fluctuating nematic phase. The dimerization has recently been proposed to be related to Berry phases of the quantum fluctuations. We use an effective low-energy description to calculate the scaling of the dimerization according to this theory, and then verify the predictions using large scale density-matrix renormalization group (DMRG) simulations, giving good evidence that the state is dimerized all the way up to its transition into the ferromagnetic phase. We furthermore discuss the multiplet structure found in the entanglement spectrum of the ground state wave functions.

preprint2014arXiv

Distinct magnetic phase transition at the surface of an antiferromagnet

In the majority of magnetic systems the surface is required to order at the same temperature as the bulk. In the present study, we report a distinct and unexpected surface magnetic phase transition, uniquely at a lower temperature than the Néel temperature. Employing grazing incidence X-ray resonant magnetic scattering we have observed the near surface behavior of uranium dioxide. UO$_2$ is a non-collinear, triple-q, antiferromagnet with the U ions on an face-centered-cubic lattice. Theoretical investigations establish that at the surface the energy increase, due to the lost bonds, reduces when the spins near the surface rotate, gradually losing their normal to the surface component. At the surface the lowest-energy spin configuration has a double-q (planar) structure. With increasing temperature, thermal fluctuations saturate the in-plane crystal field anisotropy at the surface, leading to soft excitations that have ferromagnetic $XY$ character and are decoupled from the bulk. The structure factor of a finite two-dimensional $XY$ model, fits the experimental data well for several orders of magnitude of the scattered intensity. Our results support a distinct magnetic transition at the surface in the Kosterlitz-Thouless universality class.

preprint2014arXiv

Double Semion Phase in an Exactly Solvable Quantum Dimer Model on the Kagome Lattice

Quantum dimer models typically arise in various low energy theories like those of frustrated antiferromagnets. We introduce a quantum dimer model on the kagome lattice which stabilizes an alternative $\mathbb{Z}_2$ topological order, namely the so-called "double semion" order. For a particular set of parameters, the model is exactly solvable, allowing us to access the ground state as well as the excited states. We show that the double semion phase is stable over a wide range of parameters using numerical exact diagonalization. Furthermore, we propose a simple microscopic spin Hamiltonian for which the low-energy physics is described by the derived quantum dimer model.

preprint2014arXiv

Flux insertion, entanglement, and quantized responses

We discuss the construction of a ground state wavefunction on a finite ring given the ground state on a long chain. In the presence of local symmetries, we can obtain the ground state with arbitrary flux inserted through the ring. A key ingredient are the quantum numbers of the entanglement spectrum. This method allows us to characterize phases by measuring quantized responses, such as the Hall conductance, using data contained in the entanglement spectrum. As concrete examples, we show how the Berry phase allows us to map out the phase diagram of a spin-1 model and calculate the Hall conductivity of a quantum Hall system.

preprint2014arXiv

Fragility of Symmetry Protected Topological Order on a Hubbard Ladder

Anfuso and Rosch [Phys. Rev. B 75, 144420 (2007)] showed that the "topological" Haldane phase in a fermionic spin-1/2 ladder can be continuously deformed into a "trivial" phase without explicitly breaking symmetries when local charge fluctuations are taken into account. Within the framework of symmetry protected topological phases, we revisit the model and demonstrate how the Haldane phase can be adiabatically connected to a trivial phase due to charge fluctuations. Furthermore, we show that the Haldane phase remains stable as long as the system is symmetric under particular reflection symmetries.

preprint2014arXiv

From fractionally charged solitons to Majorana bound states in a one-dimensional interacting model

We consider one-dimensional topological insulators hosting fractionally charged midgap states in the presence and absence of induced superconductivity pairing. Under the protection of a discrete symmetry, relating positive and negative energy states, the solitonic midgap states remain pinned at zero energy when superconducting correlations are induced by proximity effect. When the superconducting pairing dominates the initial insulating gap, Majorana fermion phases develop for a class of insulators. As a concrete example, we study the Creutz model with induced s-wave superconductivity and repulsive Hubbard-type interactions. For a finite wire, without interactions, the solitonic modes originating from the nonsuperconducting model survive at zero energy, revealing a fourfold-degenerate ground state. However, interactions break the aforementioned discrete symmetry and completely remove this degeneracy, thereby producing a unique ground state which ischaracterized by a topological bulk invariant with respect to the product of fermion parity and bond inversion. In contrast, the Majorana edge modes are globally robust to interactions. Moreover, the parameter range for which a topological Majorana phase is stabilized expands when increasing the repulsive Hubbard interaction. The topological phase diagram of the interacting model is obtained using a combination of mean-field theory and density matrix renormalization group techniques.

preprint2014arXiv

Interplay of charge and spin fluctuations of strongly interacting electrons on the kagome lattice

We study electrons hopping on a kagome lattice at third filling described by an extended Hubbard Hamiltonian with on-site and nearest-neighbour repulsions in the strongly correlated limit. As a consequence of the commensurate filling and the large interactions, each triangle has precisely two electrons in the effective low energy description, and these electrons form chains of different lengths. The effective Hamiltonian includes the ring exchange around the hexagons as well as the nearest-neighbor Heisenberg interaction. Using large scale exact diagonalization, we find that the effective model exhibits two different phases: If the charge fluctuations are small, the magnetic fluctuations confine the charges to short loops around hexagons, yielding a gapped charge ordered phase. When the charge fluctuations dominate, the system undergoes a quantum phase transition to a resonating plaquette phase with ordered spins and gapless spin excitations. We find that a peculiar conservation law is fulfilled: the electron in the chains can be divided into two sublattices, and this division is conserved by the ring exchange term.

preprint2014arXiv

Many-body localization in a disordered quantum Ising chain

Many-body localization occurs in isolated quantum systems when Anderson localization persists in the presence of finite interactions. Despite strong evidence for the existence of a many-body localization transition a reliable extraction of the critical disorder strength is difficult due to a large drift with system size in the studied quantities. In this work we explore two entanglement properties that are promising for the study of the manybody localization transition: the variance of the half-chain entanglement entropy of exact eigenstates and the long time change in entanglement after a local quench from an exact eigenstate. We investigate these quantities in a disordered quantum Ising chain and use them to estimate the critical disorder strength and its energy dependence. In addition, we analyze a spin-glass transition at large disorder strength and provide evidence for it being a separate transition. We thereby give numerical support for a recently proposed phase diagram of many-body localization with localization protected quantum order [Huse et al. Phys. Rev. B 88, 014206 (2013)].

preprint2014arXiv

Numerical study of a transition between Z2 topologically ordered phases

Distinguishing different topologically ordered phases and characterizing phase transitions between them is a difficult task due to the absence of local order parameters. In this paper, we use a combination of analytical and numerical approaches to distinguish two such phases and characterize a phase transition between them. The "toric code" and "double semion" models are simple lattice models exhibiting Z2 topological order. Although both models express the same topological ground state degeneracies and entanglement entropies, they are distinct phases of matter because their emergent quasi-particles obey different statistics. For a 1D model, we tune a phase transition between these two phases and obtain an exact solution to the entire phase diagram, finding a second-order Ising x Ising transition. We then use exact diagonalization to study the 2D case and find indications of a first-order transition. We show that the quasi-particle statistics provides a robust indicator of the distinct topological orders throughout the whole phase diagram.

preprint2014arXiv

Real-time dynamics in the one-dimensional Hubbard model

We consider single-particle properties in the one-dimensional repulsive Hubbard model at commensurate fillings in the metallic phase. We determine the real-time evolution of the retarded Green's function by matrix-product state methods. We find that at sufficiently late times the numerical results are in good agreement with predictions of nonlinear Luttinger liquid theory. We argue that combining the two methods provides a way of determining the single-particle spectral function with very high frequency resolution.

preprint2014arXiv

Time-evolving a matrix product state with long-ranged interactions

We introduce a numerical algorithm to simulate the time evolution of a matrix product state under a long-ranged Hamiltonian. In the effectively one-dimensional representation of a system by matrix product states, long-ranged interactions are necessary to simulate not just many physical interactions but also higher-dimensional problems with short-ranged interactions. Since our method overcomes the restriction to short-ranged Hamiltonians of most existing methods, it proves particularly useful for studying the dynamics of both power-law interacting one-dimensional systems, such as Coulombic and dipolar systems, and quasi two-dimensional systems, such as strips or cylinders. First, we benchmark the method by verifying a long-standing theoretical prediction for the dynamical correlation functions of the Haldane-Shastry model. Second, we simulate the time evolution of an expanding cloud of particles in the two-dimensional Bose-Hubbard model, a subject of several recent experiments.

preprint2013arXiv

Coulombic Charge Ice

We consider a classical model of charges $\pm q$ on a pyrochlore lattice in the presence of long range Coulomb interactions. This model first appeared in the early literature on charge order in magnetite. In the limit where the interactions become short-ranged, the model has a ground state with an extensive entropy and dipolar charge-charge correlations. When long range interactions are introduced, the exact degeneracy is broken. We study the thermodynamics of the model and show the presence of a correlated charge liquid within a temperature window in which the physics is well described as a liquid of screened charged defects. The structure factor in this phase, which has smeared pinch points at the reciprocal lattice points, may be used to detect charge ice experimentally. In addition, the model exhibits fractionally charged excitations $\pm q/2$ which are shown to interact via a $1/r$ potential. At lower temperatures, the model exhibits a transition to a long-range ordered phase. We are able to treat the Coulombic charge ice model and the dipolar spin ice model on an equal footing by mapping both to a constrained charge model on the diamond lattice. We find that states of the two ice models are related by a staggering field which is reflected in the energetics of these two models. From this perspective we can understand the origin of the spin ice and charge ice ground states as coming from a dipolar model on a diamond lattice. We study the properties of charge ice in an external electric field finding that the correlated liquid is robust to the presence of a field in contrast to the case of spin ice in a magnetic field. Finally, we comment on the transport properties of Coulombic charge ice in the correlated liquid phase.

preprint2013arXiv

Detection of Symmetry Enriched Topological Phases

Topologically ordered systems in the presence of symmetries can exhibit new structures which are referred to as symmetry enriched topological (SET) phases. We introduce simple methods to detect the SET order directly from a complete set of topologically degenerate ground state wave functions. In particular, we first show how to directly determine the characteristic symmetry fractionalization of the quasiparticles from the reduced density matrix of the minimally entangled states. Second, we show how a simple generalization of a non-local order parameter can be measured to detect SETs. The usefulness of the proposed approached is demonstrated by examining two concrete model states which exhibit SET: (i) a spin-1 model on the honeycomb lattice and (ii) the resonating valence bond state on a kagome lattice. We conclude that the spin-1 model and the RVB state are in the same SET phases.

preprint2013arXiv

Exotic Ising dynamics in a Bose-Hubbard model

We explore the dynamical properties of a one-dimensional Bose-Hubbard model, where two bosonic species interact via Feshbach resonance. We focus on the region in the phase diagram which is described by an effective, low-energy ferromagnetic Ising model in both transverse and longitudinal fields. In this regime, we numerically calculate the dynamical structure factor of the Bose-Hubbard model using the time-evolving block decimation method. In the ferromagnetic phase, we observe both the continuum of excitations and the bound states in the presence of a longitudinal field. Near the Ising critical point, we observe the celebrated E8 mass spectrum in the excited states. We also point out possible measurements which could be used to detect these excitations in an optical lattice experiment.

preprint2013arXiv

Loschmidt echo and the many-body orthogonality catastrophe in a qubit-coupled Luttinger liquid

We investigate the many-body generalization of the orthogonality catastrophe by studying the generalized Loschmidt echo of Luttinger liquids (LLs) after a global change of interaction. It decays exponentially with system size and exhibits universal behaviour: the steady state exponent after quenching back and forth n-times between 2 LLs (bang-bang protocol) is 2n-times bigger than that of the adiabatic overlap, and depends only on the initial and final LL parameters. These are corroborated numerically by matrix-product state based methods of the XXZ Heisenberg model. An experimental setup consisting of a hybrid system containing cold atoms and a flux qubit coupled to a Feshbach resonance is proposed to measure the Loschmidt echo using rf spectroscopy or Ramsey interferometry.

preprint2013arXiv

Phase diagram of the isotropic spin-3/2 model on the z=3 Bethe lattice

We study an SU(2) symmetric spin-3/2 model on the z=3 Bethe lattice using the infinite Time Evolving Block Decimation (iTEBD) method. This model is shown to exhibit a rich phase diagram. We compute the expectation values of several order parameters which allow us to identify a ferromagnetic, a ferrimagnetic, a anti-ferromagnetic as well as a dimerized phase. We calculate the entanglement spectra from which we conclude the existence of a symmetry protected topological phase that is characterized by S=1/2 edge spins. Details of the iTEBD algorithm used for the simulations are included.

preprint2013arXiv

Topological characterization of fractional quantum Hall ground states from microscopic Hamiltonians

We show how to numerically calculate several quantities that characterize topological order starting from a microscopic fractional quantum Hall (FQH) Hamiltonian. To find the set of degenerate ground states, we employ the infinite density matrix renormalization group (iDMRG) method based on the matrix-product state (MPS) representation of FQH states on an infinite cylinder. To study localized quasiparticles of a chosen topological charge, we use pairs of degenerate ground states as boundary conditions for the iDMRG. We then show that the wave function obtained on the infinite cylinder geometry can be adapted to a torus of arbitrary modular parameter, which allows us to explicitly calculate the non-Abelian Berry connection associated with the modular T-transformation. As a result, the quantum dimensions, topological spins, quasiparticle charges, chiral central charge, and Hall viscosity of the phase can be obtained using data contained entirely in the entanglement spectrum of an infinite cylinder.

preprint2013arXiv

Topological phases in gapped edges of fractionalized systems

Recently, it has been proposed that exotic one-dimensional phases can be realized by gapping out the edge states of a fractional topological insulator. The low-energy edge degrees of freedom are described by a chain of coupled parafermions. We introduce a classification scheme for the phases that can occur in parafermionic chains. We find that the parafermions support both topological symmetry fractionalized phases as well as phases in which the parafermions condense. In the presence of additional symmetries, the phases form a non-Abelian group. As a concrete example of the classification, we consider the effective edge model for a $ν= 1/3$ fractional topological insulator for which we calculate the entanglement spectra numerically and show that all possible predicted phases can be realized.

preprint2012arXiv

Detection of Symmetry Protected Topological Phases in 1D

A topological phase is a phase of matter which cannot be characterized by a local order parameter. It has been shown that gapped phases in 1D systems can be completely characterized using tools related to projective representations of the symmetry groups. We show how to determine the matrices of these representations in a simple way in order to distinguish between different phases directly. From these matrices we also point out how to derive several different types of non-local order parameters for time reversal, inversion symmetry and $Z_2 \times Z_2$ symmetry, as well as some more general cases (some of which have been obtained before by other methods). Using these concepts, the ordinary string order for the Haldane phase can be related to a selection rule that changes at the critical point. We furthermore point out an example of a more complicated internal symmetry for which the ordinary string order cannot be applied.

preprint2012arXiv

Linear quantum quench in the Heisenberg XXZ chain: time dependent Luttinger model description of a lattice system

We study variable-rate linear quenches in the anisotropic Heisenberg (XXZ) chain, starting at the XX point. This is equivalent to switching on a nearest neighbour interaction for hard-core bosons or an interaction quench for free fermions. The physical observables we investigate are: the energy pumped into the system during the quench, the spin-flip correlation function, and the bipartite fluctuations of the z component of the spin in a box. We find excellent agreement between exact numerics (infinite system time-evolving block decimation, iTEBD) and analytical results from bosonization, as a function of the quench time, spatial coordinate and interaction strength. This provides a stringent and much-needed test of Luttinger liquid theory in a non-equilibrium situation.

preprint2012arXiv

Matrix-Product based Projected Wave Functions Ansatz for Quantum Many-Body Ground States

We develop a new projected wave function approach which is based on projection operators in the form of matrix-product operators (MPOs). Our approach allows to variationally improve the short range entanglement of a given trial wave function by optimizing the matrix elements of the MPOs while the long range entanglement is contained in the initial guess of the wave function. The optimization is performed using standard variational Monte Carlo techniques. We demonstrate the efficiency of our approach by considering a one-dimension model of interacting spinless fermions. In addition, we indicate how to generalize this approach to higher dimensions using projection operators which are based on tensor products.

preprint2012arXiv

Symmetry protection of topological order in one-dimensional quantum spin systems

We discuss the characterization and stability of the Haldane phase in integer spin chains on the basis of simple, physical arguments. We find that an odd-$S$ Haldane phase is a topologically non-trivial phase which is protected by any one of the following three global symmetries: (i) the dihedral group of $π$-rotations about $x,y$ and $z$ axes; (ii) time-reversal symmetry $S^{x,y,z} \rightarrow - S^{x,y,z}$; (iii) link inversion symmetry (reflection about a bond center), consistently with previous results [Phys. Rev. B \textbf{81}, 064439 (2010)]. On the other hand, an even-$S$ Haldane phase is not topologically protected (i.e., it is indistinct from a trivial, site-factorizable phase). We show some numerical evidence that supports these claims, using concrete examples.

preprint2012arXiv

The phase diagram of the anisotropic spin-2 XXZ model: an infinite system DMRG study

We study the ground state phase diagram of the quantum spin-2 XXZ chain in the presence of on-site anisotropy using a matrix-product state based infinite system density-matrix-renormalization-group (iDMRG) algorithm. One of the interests in this system is in connecting the highly quantum mechanical spin-1 phase diagram with the classical S=infinity phase diagram. Several of the recent advances within DMRG make it possible to perform a detailed analysis of the whole phase diagram. We consider different types of on-site anisotropies which allows us to establish the validity of the following statements: One, the spin-2 model can be tuned into a phase which is equivalent to the "topologically nontrivial" spin-1 Haldane phase. Two, the spin-2 Haldane phase at the isotropic Heisenberg point is adiabatically connected to the "trivial" large-D phase, with a continuous change of the Hamiltonian parameters. Furthermore, we study the spin-3 XXZ chain to help explain the development of the classical phase diagram. We present details on how to use the iDMRG method to map out the phase diagram and include an extensive discussion of the numerical methods.

preprint2012arXiv

Topological Phases of One-Dimensional Fermions: An Entanglement Point of View

The effect of interactions on topological insulators and superconductors remains, to a large extent, an open problem. Here, we describe a framework for classifying phases of one-dimensional interacting fermions, focusing on spinless fermions with time-reversal symmetry and particle number parity conservation, using concepts of entanglement. In agreement with an example presented by Fidkowski \emph{et. al.} (Phys. Rev. B 81, 134509 (2010)), we find that in the presence of interactions there are only eight distinct phases, which obey a $\mathbb{Z}_8$ group structure. This is in contrast to the $\mathbb{Z}$ classification in the non-interacting case. Each of these eight phases is characterized by a unique set of bulk invariants, related to the transformation laws of its entanglement (Schmidt) eigenstates under symmetry operations, and has a characteristic degeneracy of its entanglement levels. If translational symmetry is present, the number of distinct phases increases to 16.

preprint2012arXiv

Unbounded growth of entanglement in models of many-body localization

An important and incompletely answered question is whether a closed quantum system of many interacting particles can be localized by disorder. The time evolution of simple (unentangled) initial states is studied numerically for a system of interacting spinless fermions in one dimension described by the random-field XXZ Hamiltonian. Interactions induce a dramatic change in the propagation of entanglement and a smaller change in the propagation of particles. For even weak interactions, when the system is thought to be in a many-body localized phase, entanglement shows neither localized nor diffusive behavior but grows without limit in an infinite system: interactions act as a singular perturbation on the localized state with no interactions. The significance for proposed atomic experiments is that local measurements will show a large but nonthermal entropy in the many-body localized state. This entropy develops slowly (approximately logarithmically) over a diverging time scale as in glassy systems.

preprint2011arXiv

Bound states and E_8 symmetry effects in perturbed quantum Ising chains

In a recent experiment on CoNb_2O_6, Coldea et al. [Science 327, 177 (2010)] found for the first time experimental evidence of the exceptional Lie algebra E_8. The emergence of this symmetry was theoretically predicted long ago for the transverse quantum Ising chain in the presence of a weak longitudinal field. We consider an accurate microscopic model of CoNb_2O_6 incorporating additional couplings and calculate numerically the dynamical structure function using a recently developed matrix-product-state method. The excitation spectra show bound states characteristic of the weakly broken E_8 symmetry. We compare the observed bound state signatures in this model to those found in the transverse Ising chain in a longitudinal field and to experimental data.

preprint2011arXiv

Extended quantum U(1)-liquid phase in a three-dimensional quantum dimer model

Recently, quantum dimer models, in which the system can tunnel between different classical dimer configurations, have attracted a great deal of interest as a paradigm for the study of exotic quantum phases. Much of this excitement has centred on the claim that a certain class of quantum dimer model, defined on a bipartite lattice, can support a quantum U(1)-liquid phase with deconfined fractional excitations in three dimensions. These fractional monomer excitations are quantum analogues of the magnetic monopoles found in spin ice. In this article we use extensive quantum Monte Carlo simulations to establish the ground-state phase diagram of the quantum dimer model on the three-dimensional, bipartite, diamond lattice as a function of the ratio μ of the potential to kinetic energy terms in the Hamiltonian. We find that, for μ_c = 0.75 +/- 0.04, the model undergoes a first-order quantum phase transition from an ordered "R-state" into an extended quantum U(1)-liquid phase, which terminates in a quantum critical "RK point" for μ=1. This confirms the published field-theoretical scenario. We present detailed evidence for the existence of the U(1)-liquid phase, and indirect evidence for the existence of its photon and monopole excitations. We also explore some of the technical ramifications of this analysis, benchmarking quantum Monte Carlo against a variety of exact and perturbative results, comparing different variational wave functions. The ergodicity of the quantum dimer model on a diamond lattice is discussed in detail. These results complete and extend the analysis previously published in [O. Sikora et al., Phys. Rev. Lett. 103, 247001 (2009)].

preprint2011arXiv

Fermionic quantum dimer and fully-packed loop models on the square lattice

We consider fermionic fully-packed loop and quantum dimer models which serve as effective low-energy models for strongly correlated fermions on a checkerboard lattice at half and quarter filling, respectively. We identify a large number of fluctuationless states specific to each case, due to the fermionic statistics. We discuss the symmetries and conserved quantities of the system and show that for a class of fluctuating states in the half-filling case, the fermionic sign problem can be gauged away. This claim is supported by numerical evaluation of the low-lying states and can be understood by means of an algebraic construction. The elimination of the sign problem then allows us to analyze excitations at the Rokhsar-Kivelson point of the models using the relation to the height model and its excitations, within the single-mode approximation. We then discuss a mapping to a U(1) lattice gauge theory which relates the considered low-energy model to the compact quantum electrodynamics in 2+1 dimensions. Furthermore, we point out consequences and open questions in the light of these results.

preprint2011arXiv

Quantum Ice : a quantum Monte Carlo study

Ice states, in which frustrated interactions lead to a macroscopic ground-state degeneracy, occur in water ice, in problems of frustrated charge order on the pyrochlore lattice, and in the family of rare-earth magnets collectively known as spin ice. Of particular interest at the moment are "quantum spin ice" materials, where large quantum fluctuations may permit tunnelling between a macroscopic number of different classical ground states. Here we use zero-temperature quantum Monte Carlo simulations to show how such tunnelling can lift the degeneracy of a spin or charge ice, stabilising a unique "quantum ice" ground state --- a quantum liquid with excitations described by the Maxwell action of 3+1-dimensional quantum electrodynamics. We further identify a competing ordered "squiggle" state, and show how both squiggle and quantum ice states might be distinguished in neutron scattering experiments on a spin ice material.

preprint2011arXiv

Supersolid phase and magnetization plateaus observed in anisotropic spin-3/2 Heisenberg model on bipartite lattices

We study the spin-3/2 Heisenberg model including easy-plane and exchange anisotropies in one and two dimensions. In the Ising limit, when the off-diagonal exchange interaction J is zero, the phase diagram in magnetic field is characterized by magnetization plateaus that are either translationally invariant or have a two-sublattice order, with phase boundaries that are macroscopically degenerate. Using a site factorized variational wave function and perturbational expansion around the Ising limit, we find that superfluid and supersolid phases emerge between the plateaus for small finite values of J. The variational approach is complemented by a Density Matrix Renormalization Group study of a one-dimensional chain and exact diagonalization calculations on small clusters of a square lattice. The studied model may serve as a minimal model for the layered Ba2CoGe2O7 material compound, and we believe that the vicinity of the uniform 1/3 plateau in the model parameter space can be observed as an anomaly in the measured magnetization curve.

preprint2010arXiv

Entanglement spectrum of a topological phase in one dimension

We show that the Haldane phase of S=1 chains is characterized by a double degeneracy of the entanglement spectrum. The degeneracy is protected by a set of symmetries (either the dihedral group of $π$-rotations about two orthogonal axes, time-reversal symmetry, or bond centered inversion symmetry), and cannot be lifted unless either a phase boundary to another, "topologically trivial", phase is crossed, or the symmetry is broken. More generally, these results offer a scheme to classify gapped phases of one dimensional systems. Physically, the degeneracy of the entanglement spectrum can be observed by adiabatically weakening a bond to zero, which leaves the two disconnected halves of the system in a finitely entangled state.

preprint2010arXiv

Strongly correlated fermions on a kagome lattice

We study a model of strongly correlated spinless fermions on a kagome lattice at 1/3 filling, with interactions described by an extended Hubbard Hamiltonian. An effective Hamiltonian in the desired strong correlation regime is derived, from which the spectral functions are calculated by means of exact diagonalization techniques. We present our numerical results with a view to discussion of possible signatures of confinement/deconfinement of fractional charges.

preprint2009arXiv

Dynamics after a sweep through a quantum critical point

The coherent quantum evolution of a one-dimensional many-particle system after sweeping the Hamiltonian through a critical point is studied using a generalized quantum Ising model containing both integrable and non-integrable regimes. It is known from previous work that universal power laws appear in such quantities as the mean number of excitations created by the sweep. Several other phenomena are found that are not reflected by such averages: there are two scaling regimes of the entanglement entropy and a relaxation that is power-law rather than exponential. The final state of evolution after the quench is not well characterized by any effective temperature, and the Loschmidt echo converges algebraically to a constant for long times, with cusplike singularities in the integrable case that are dynamically broadened by nonintegrable perturbations.

preprint2009arXiv

Entanglement spectra of critical and near-critical systems in one dimension

The entanglement spectrum of a pure state of a bipartite system is the full set of eigenvalues of the reduced density matrix obtained from tracing out one part. Such spectra are known in several cases to contain important information beyond that in the entanglement entropy. This paper studies the entanglement spectrum for a variety of critical and near-critical quantum lattice models in one dimension, chiefly by the iTEBD numerical method, which enables both integrable and non-integrable models to be studied. We find that the distribution of eigenvalues in the entanglement spectra agrees with an approximate result derived by Calabrese and Lefevre to an accuracy of a few percent for all models studied. This result applies whether the correlation length is intrinsic or generated by the finite matrix size accessible in iTEBD. For the transverse Ising model, the known exact results for the entanglement spectrum are used to confirm the validity of the iTEBD approach. For more general models, no exact result is available but the iTEBD results directly test the hypothesis that all moments of the reduced density matrix are determined by a single parameter.

preprint2009arXiv

Theory of finite-entanglement scaling at one-dimensional quantum critical points

Studies of entanglement in many-particle systems suggest that most quantum critical ground states have infinitely more entanglement than non-critical states. Standard algorithms for one-dimensional many-particle systems construct model states with limited entanglement, which are a worse approximation to quantum critical states than to others. We give a quantitative theory of previously observed scaling behavior resulting from finite entanglement at quantum criticality: the scaling theory of finite entanglement is only superficially similar to finite-size scaling, and has a different physical origin. We find that finite-entanglement scaling is governed not by the scaling dimension of an operator but by the "central charge" of the critical point, which counts its universal degrees of freedom. An important ingredient is the recently obtained universal distribution of density-matrix eigenvalues at a critical point\cite{calabrese1}. The parameter-free theory is checked against numerical scaling at several quantum critical points.