The ergodicity postulate, a foundational pillar of Gibbsian statistical mechanics predicts that a periodically driven (Floquet) system in the absence of any conservation law heats to a featureless `infinite temperature' state. Here, we find--for a clean and interacting generic spin chain subject to a \it strong driving field--that this can be prevented by the emergence of \it approximate but stable conservation-laws not present in the undriven system. We identify their origin: they do not necessarily owe their stability to familiar protections by symmetry, topology, disorder, or even high energy costs. We show numerically, \it in the thermodynamic limit, that when required by these emergent conservation-laws, the entanglement-entropy density of an infinite subsystem remains zero over our entire simulation time of several decades in natural units. We further provide a recipe for designing such conservation laws with high accuracy. Finally, we present an ensemble description, which we call the strongly driven ensemble incorporating these constraints. This provides a way to control many-body chaos through stable Floquet-engineering. Strong signatures of these conservation-laws should be experimentally accessible since they manifest in all length and time scales. Variants of the spin model we have used, have already been realized using Rydberg-dressed atoms.
The realization of synthetic gauge fields for charge neutral ultracold atoms and the simulation of quantum Hall physics has witnessed remarkable experimental progress. Here, we establish key signatures of fractional quantum Hall systems in their non-equilibrium quantum dynamics. We show that in the lowest Landau level the system generically relaxes subdiffusively. The slow relaxation is understood from emergent conservation laws of the total charge and the associated dipole moment that arise from the effective Hamiltonian projected onto the lowest Landau level, leading to subdiffusive fracton hydrodynamics. We discuss the prospect of rotating quantum gases as well as ultracold atoms in optical lattices for observing this unconventional relaxation dynamics.
Tyler A. Cochran, Bernhard Jobst, Eliott Rosenberg, Yuri D. Lensky, Gaurav Gyawali, Norhan Eassa, Melissa Will, Dmitry Abanin, Rajeev Acharya, Laleh Aghababaie Beni, Trond I. Andersen, Markus Ansmann, Frank Arute, Kunal Arya, Abraham Asfaw, Juan Atalaya, Ryan Babbush, Brian Ballard, Joseph C. Bardin, Andreas Bengtsson, et al (172) Lattice gauge theories (LGTs) can be employed to understand a wide range of phenomena, from elementary particle scattering in high-energy physics to effective descriptions of many-body interactions in materials. Studying dynamical properties of emergent phases can be challenging as it requires solving many-body problems that are generally beyond perturbative limits. We investigate the dynamics of local excitations in a $\mathbb{Z}_2$ LGT using a two-dimensional lattice of superconducting qubits. We first construct a simple variational circuit which prepares low-energy states that have a large overlap with the ground state; then we create particles with local gates and simulate their quantum dynamics via a discretized time evolution. As the effective magnetic field is increased, our measurements show signatures of transitioning from deconfined to confined dynamics. For confined excitations, the magnetic field induces a tension in the string connecting them. Our method allows us to experimentally image string dynamics in a (2+1)D LGT from which we uncover two distinct regimes inside the confining phase: for weak confinement the string fluctuates strongly in the transverse direction, while for strong confinement transverse fluctuations are effectively frozen. In addition, we demonstrate a resonance condition at which dynamical string breaking is facilitated. Our LGT implementation on a quantum processor presents a novel set of techniques for investigating emergent particle and string dynamics.
An important result in the theory of quantum control is the "universality" of $2$-local unitary gates, i.e. the fact that any global unitary evolution of a system of $L$ qudits can be implemented by composition of $2$-local unitary gates. Surprisingly, recent results have shown that universality can break down in the presence of symmetries: in general, not all globally symmetric unitaries can be constructed using $k$-local symmetric unitary gates. This also restricts the dynamics that can be implemented by symmetric local Hamiltonians. In this paper, we show that obstructions to universality in such settings can in general be understood in terms of superoperator symmetries associated with unitary evolution by restricted sets of gates. These superoperator symmetries lead to block decompositions of the operator Hilbert space, which dictate the connectivity of operator space, and hence the structure of the dynamical Lie algebra. We demonstrate this explicitly in several examples by systematically deriving the superoperator symmetries from the gate structure using the framework of commutant algebras, which has been used to systematically derive symmetries in other quantum many-body systems. We clearly delineate two different types of non-universality, which stem from different structures of the superoperator symmetries, and discuss its signatures in physical observables. In all, our work establishes a comprehensive framework to explore the universality of unitary circuits and derive physical consequences of its absence.
The simulation of quantum many-body systems, relevant for quantum chemistry and condensed matter physics, is one of the most promising applications of near-term quantum computers before fault-tolerance. However, since the vast majority of quantum computing technologies are built around qubits and discrete gate-based operations, the translation of the physical problem into this framework is a crucial step. This translation will often be device specific, and a suboptimal implementation will be punished by the exponential compounding of errors on real devices. The importance of an efficient mapping is already revealed for models of spinful fermions in two or three dimensions, which naturally arise when the relevant physics relates to electrons. Using the most direct and well-known mapping, the Jordan-Wigner transformation, leads to a non-local representation of local degrees of freedom, and necessities efficient decompositions of non-local unitary gates into a sequence of hardware accessible local gates. In this paper, we provide a step-by-step recipe for simulating the paradigmatic two-dimensional Fermi-Hubbard model on a quantum computer using only local operations. To provide the ingredients for such a recipe, we briefly review the plethora of different approaches that have emerged recently but focus on the Derby-Klassen compact fermion mapping in order to make our discussion concrete. We provide a detailed recipe for an end-to-end simulation including embedding on a physical device, preparing initial states such as ground states, simulation of unitary time evolution, and measurement of observables and spectral functions. We explicitly compute the resource requirements for simulating a global quantum quench and conclude by discussing the challenges and future directions for simulating strongly-correlated fermionic matter on quantum computers.
We find that the presence of strong non-Abelian conserved quantities can lead to highly entangled stationary states even for unital quantum channels. We derive exact expressions for the bipartite logarithmic negativity, Rényi negativities, and operator space entanglement for stationary states restricted to one symmetric subspace, with focus on the trivial subspace. We prove that these apply to open quantum evolutions whose commutants, characterizing all strongly conserved quantities, correspond to either the universal enveloping algebra of a Lie algebra or to the Read-Saleur commutants. The latter provides an example of quantum fragmentation, whose dimension is exponentially large in system size. We find a general upper bound for all these quantities given by the logarithm of the dimension of the commutant on the smaller bipartition of the chain. As Abelian examples, we show that strong U($1$) symmetries and classical fragmentation lead to separable stationary states in any symmetric subspace. In contrast, for non-Abelian SU$(N)$ symmetries, both logarithmic and Rényi negativities scale logarithmically with system size. Finally, we prove that while Rényi negativities with $n>2$ scale logarithmically with system size, the logarithmic negativity (as well as generalized Rényi negativities with $n<2$) exhibits a volume law scaling for the Read-Saleur commutants. Our derivations rely on the commutant possessing a Hopf algebra structure in the limit of infinitely large systems, and hence also apply to finite groups and quantum groups.
The relaxation behaviour of isolated quantum systems taken out of equilibrium is among the most intriguing questions in many-body physics. Quantum systems out of equilibrium typically relax to thermal equilibrium states by scrambling local information and building up entanglement entropy. However, kinetic constraints in the Hamiltonian can lead to a breakdown of this fundamental paradigm due to a fragmentation of the underlying Hilbert space into dynamically decoupled subsectors in which thermalisation can be strongly suppressed. Here, we experimentally observe Hilbert space fragmentation (HSF) in a two-dimensional tilted Bose-Hubbard model. Using quantum gas microscopy, we engineer a wide variety of initial states and find a rich set of manifestations of HSF involving bulk states, interfaces and defects, i.e., d = 2, 1 and 0 dimensional objects. Specifically, uniform initial states with equal particle number and energy differ strikingly in their relaxation dynamics. Inserting controlled defects on top of a global, non-thermalising chequerboard state, we observe highly anisotropic, sub-dimensional dynamics, an immediate signature of their fractonic nature. An interface between localized and thermalising states in turn displays dynamics depending on its orientation. Our results mark the first observation of HSF beyond one dimension, as well as the concomitant direct observation of fractons, and pave the way for in-depth studies of microscopic transport phenomena in constrained systems
Entanglement of mixed quantum states can be quantified using the partial transpose and its corresponding entanglement measure, the logarithmic negativity. Recently, the notion of partial transpose has been extended to systems of anyons, which are exotic quasiparticles whose exchange statistics go beyond the bosonic and fermionic case. Studying the fundamental properties of this anyonic partial transpose, we first reveal that when applied to the special case of fermionic systems, it can be reduced to the fermionic partial transpose or its twisted variant depending on whether or not a boundary Majorana fermion is present. Focusing on ground state properties, we find that the anyonic partial transpose captures both the correct entanglement scaling for gapless systems, as predicted by conformal field theory, and the phase transition between a topologically trivial and a nontrivial phase. For non-abelian anyons and the bipartition geometry, we find a rich multiplet structure in the eigenvalues of the partial transpose, the so-called negativity spectrum, and reveal the possibility of defining both a charge- and an imbalance-resolved negativity.
The potential impact of quantum machine learning algorithms on industrial applications remains an exciting open question. Conventional methods for encoding classical data into quantum computers are not only too costly for a potential quantum advantage in the algorithms but also severely limit the scale of feasible experiments on current hardware. Therefore, recent works, despite claiming the near-term suitability of their algorithms, do not provide experimental benchmarking on standard machine learning datasets. We attempt to solve the data encoding problem by improving a recently proposed variational algorithm [1] that approximately prepares the encoded data, using asymptotically shallow circuits that fit the native gate set and topology of currently available quantum computers. We apply the improved algorithm to encode the Fashion-MNIST dataset [2], which can be directly used in future empirical studies of quantum machine learning algorithms. We deploy simple quantum variational classifiers trained on the encoded dataset on a current quantum computer ibmq-kolkata [3] and achieve moderate accuracies, providing a proof of concept for the near-term usability of our data encoding method.
The advent of near-term digital quantum computers could offer us an exciting opportunity to investigate quantum many-body phenomena beyond that of classical computing. To make the best use of the hardware available, it is paramount that we have methods that accurately simulate Hamiltonian dynamics for limited circuit depths. In this paper, we propose a method to classically optimise unitary brickwall circuits to approximate quantum time evolution operators. Our method is scalable in system size through the use of tensor networks. We demonstrate that, for various three-body Hamiltonians, our approach produces quantum circuits that can outperform Trotterization in both their accuracy and the quantum circuit depth needed to implement the dynamics, with the exact details being dependent on the Hamiltonian. We also explain how to choose an optimal time step that minimises the combined errors of the quantum device and the brickwall circuit approximation.
Efficient preparation of many-body ground states is key to harnessing the power of quantum computers in studying quantum many-body systems. In this work, we propose a simple method to design exact linear-depth parameterized quantum circuits which prepare a family of ground states across topological quantum phase transitions in 2D. We achieve this by constructing ground states represented by isometric tensor networks (isoTNS), which form a subclass of tensor network states that are efficiently preparable. By continuously tuning a parameter in the wavefunction, the many-body ground state undergoes quantum phase transitions, exhibiting distinct 2D quantum phases. We illustrate this by constructing an isoTNS path with bond dimension $D = 2$ interpolating between distinct symmetry-enriched topological (SET) phases. At the transition point, the wavefunction is related to a gapless point in the classical six-vertex model. Furthermore, the critical wavefunction supports a power-law correlation along one spatial direction while remaining long-range ordered in the other spatial direction. We provide an explicit parametrized local quantum circuit for the path and show that the 2D isoTNS can also be efficiently simulated by a holographic quantum algorithm requiring only an 1D array of qubits.
We explore the relationship between higher-form symmetries and entanglement properties in discrete lattice gauge theories, which can exhibit both topologically ordered phases and higher-form symmetry-protected topological (SPT) phases. Our study centers on generalizing the Fradkin-Shenker model, where the Gauss law constraint can be either emergent or exact. The phase diagram includes a topologically ordered phase and a non-trivial SPT phase protected by a 1-form and a 0-form symmetry. We obtain the following key findings: First, the entanglement properties depend on whether the 1-form symmetries and the Gauss law are exact or emergent. For the emergent Gauss law, the entanglement spectrum (ES) of the non-trivial SPT phase exhibits degeneracies, which are robust at low energies against weak perturbations that explicitly break the exact 1-form symmetry. When the Gauss law and the 1-form symmetry are both exact, the ES degeneracy is extensive. This extensive degeneracy is fragile and can be removed completely by infinitesimal perturbations that explicitly break the exact 1-form symmetry while keeping the Gauss law exact. Second, we consider the ES in the topologically ordered phase where 1-form symmetries are spontaneously broken. In contrast to the ES of the non-trivial SPT phase, we find that spontaneous higher-form symmetry breaking removes "half" of the ES levels, leading to a non-degenerate ES in the topologically ordered phase in general. Third, we derive a connection between spontaneous higher-form symmetry breaking and the topological entanglement entropy (TEE). Using this relation, we investigate the entanglement entropy that can be distilled in the deconfined phase of the original Fradkin-Shenker model using gauge-invariant measurements. We show that the TEE is robust against the measurement when the 1-form symmetry is emergent but fragile when the 1-form symmetry is exact.
Machine learning tasks are an exciting application for quantum computers, as it has been proven that they can learn certain problems more efficiently than classical ones. Applying quantum machine learning algorithms to classical data can have many important applications, as qubits allow for dealing with exponentially more data than classical bits. However, preparing the corresponding quantum states usually requires an exponential number of gates and therefore may ruin any potential quantum speedups. Here, we show that classical data with a sufficiently quickly decaying Fourier spectrum after being mapped to a quantum state can be well-approximated by states with a small Schmidt rank (i.e., matrix product states) and we derive explicit error bounds. These approximated states can, in turn, be prepared on a quantum computer with a linear number of nearest-neighbor two-qubit gates. We confirm our results numerically on a set of $1024\times1024$-pixel images taken from the 'Imagenette' dataset. Additionally, we consider different variational circuit ansätze and demonstrate numerically that one-dimensional sequential circuits achieve the same compression quality as more powerful ansätze.
Charge and energy are expected to diffuse in interacting systems of fermions at finite temperatures, even in the absence of disorder, with the interactions inducing a crossover from the coherent and ballistic streaming of quasi-particles at early times, to incoherent diffusive behavior at late times. The relevant crossover timescales and the transport coefficients are both controlled by the strength of interactions. In this work we develop a numerical method to simulate such systems at high temperatures, applicable in a wide range of interaction strengths, by adapting Dissipation-assisted Operator Evolution (DAOE) to fermions. Our fermion DAOE, which approximates the exact dynamics by systematically discarding information from high $n$-point functions, is tailored to capture non-interacting dynamics exactly, thus providing a good starting point for the weakly interacting problem. Applying our method to a microscopic model of weakly interacting fermions, we numerically demonstrate that the crossover from ballistic to diffusive transport happens at a time $t_D\sim1/\Delta^{2}$ and that the diffusion constant similarly scales as $D \sim 1/\Delta^2$, where $\Delta$ is the interaction strength. We substantiate this scaling with a Fermi's golden rule calculation in the operator spreading picture, interpreting $t_D$ as the fermion-fermion scattering time and lifetime of the single-particle Green's function.
In the study of the thermalization of closed quantum systems, the role of kinetic constraints on the temporal dynamics and the eventual thermalization is attracting significant interest. Kinetic constraints typically lead to long-lived metastable states depending on initial conditions. We consider a model of interacting hardcore bosons with an additional kinetic constraint that was originally devised to capture glassy dynamics at high densities. As a main result, we demonstrate that the system is highly prone to localization in the presence of uncorrelated disorder. Adding disorder quickly triggers long-lived dynamics as evidenced in the time evolution of density autocorrelations. Moreover, the kinetic constraint favors localization also in the eigenstates, where a finite-size transition to a many-body localized phase occurs for much lower disorder strengths than for the same model without a kinetic constraint. Our work sheds light on the intricate interplay of kinetic constraints and localization and may provide additional control over many-body localized phases in the time domain.
Tensor networks have recently found applications in machine learning for both supervised learning and unsupervised learning. The most common approaches for training these models are gradient descent methods. In this work, we consider an alternative training scheme utilizing basic tensor network operations, e.g., summation and compression. The training algorithm is based on compressing the superposition state constructed from all the training data in product state representation. The algorithm could be parallelized easily and only iterates through the dataset once. Hence, it serves as a pre-training algorithm. We benchmark the algorithm on the MNIST dataset and show reasonable results for generating new images and classification tasks. Furthermore, we provide an interpretation of the algorithm as a compressed quantum kernel density estimation for the probability amplitude of input data.
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.
Quantum phase transitions between different topologically ordered phases exhibit rich structures and are generically challenging to study in microscopic lattice models. In this work, we propose a tensor-network solvable model that allows us to tune between different symmetry enriched topological (SET) phases. Concretely, we consider a decorated two-dimensional toric code model for which the ground state can be expressed as a two-dimensional tensor-network state with bond dimension $D=3$ and two tunable parameters. We find that the time-reversal (TR) symmetric system exhibits three distinct phases (i) an SET toric code phase in which anyons transform non-trivially under TR, (ii) a toric code phase in which TR does not fractionalize, and (iii) a topologically trivial phase that is adiabatically connected to a product state. We characterize the different phases using the topological entanglement entropy and a membrane order parameter that distinguishes the two SET phases. Along the phase boundary between the SET toric code phase and the toric code phase, the model has an enhanced $U(1)$ symmetry and the ground state is a quantum critical loop gas wavefunction whose squared norm is equivalent to the partition function of the classical $O(2)$ model. By duality transformations, this tensor-network solvable model can also be used to describe transitions between SET double-semion phases and between $\mathbb{Z}_2\times\mathbb{Z}_2^T$ symmetry protected topological phases in two dimensions.
Repeated projective measurements in unitary circuits can lead to an entanglement phase transition as the measurement rate is tuned. In this work, we consider a different setting in which the projective measurements are replaced by dynamically chosen unitary gates that minimize the entanglement. This can be seen as a one-dimensional unitary circuit game in which two players get to place unitary gates on randomly assigned bonds at different rates: The "entangler" applies a random local unitary gate with the aim of generating extensive (volume law) entanglement. The "disentangler," based on limited knowledge about the state, chooses a unitary gate to reduce the entanglement entropy on the assigned bond with the goal of limiting to only finite (area law) entanglement. In order to elucidate the resulting entanglement dynamics, we consider three different scenarios: (i) a classical discrete height model, (ii) a Clifford circuit, and (iii) a general $U(4)$ unitary circuit. We find that both the classical and Clifford circuit models exhibit phase transitions as a function of the rate that the disentangler places a gate, which have similar properties that can be understood through a connection to the stochastic Fredkin chain. In contrast, the "entangler" always wins when using Haar random unitary gates and we observe extensive, volume law entanglement for all non-zero rates of entangling.
We study measurement-induced symmetry-protected topological (SPT) order in a wide class of quantum random circuit models by combining calculations within the stabilizer formalism with tensor network simulations. We construct a family of quantum random circuits, generating the out-of-equilibrium version of all generalized cluster models, and derive a set of non-local string order parameters to distinguish different SPT phases. We apply this framework to investigate a random circuit realization of the XZX cluster model, and use the string order parameter to demonstrate that the phase diagram is stable against extending the class of unitary gates in the circuit, from Clifford gates to Haar unitaries. We then turn to the XZZX generalized cluster model, and demonstrate the coexistence of SPT order and spontaneous symmetry breaking, by relying on string order parameters and a connected correlation function.
Tensor networks provide a useful tool to describe low-dimensional complex many-body systems. Finding efficient algorithms to use these methods for finite-temperature simulations in two dimensions is a continuing challenge. Here, we use the class of recently introduced isometric tensor network states, which can also be directly realized with unitary gates on a quantum computer. We utilize a purification ansatz to efficiently represent thermal states of the transverse field Ising model. By performing an imaginary-time evolution starting from infinite temperature, we find that this approach offers a different way with low computational complexity to represent states at finite temperatures.
We propose and benchmark a modified time evolution block decimation (TEBD) algorithm that uses a truncation scheme based on the QR decomposition instead of the singular value decomposition (SVD). The modification reduces the scaling with the dimension of the physical Hilbert space $d$ from $d^3$ down to $d^2$. Moreover, the QR decomposition has a lower computational complexity than the SVD and allows for highly efficient implementations on GPU hardware. In a benchmark simulation of a global quench in a quantum clock model, we observe a speedup of up to three orders of magnitude comparing QR and SVD based updates on an A100 GPU.
The exact contraction of a generic two-dimensional (2D) tensor network state (TNS) is known to be exponentially hard, making simulation of 2D systems difficult. The recently introduced class of isometric TNS (isoTNS) represents a subset of TNS that allows for efficient simulation of such systems on finite square lattices. The isoTNS ansatz requires the identification of an "orthogonality column" of tensors, within which one-dimensional matrix product state (MPS) methods can be used for calculation of observables and optimization of tensors. Here we extend isoTNS to infinitely long strip geometries and introduce an infinite version of the Moses Move algorithm for moving the orthogonality column around the network. Using this algorithm, we iteratively transform an infinite MPS representation of a 2D quantum state into a strip isoTNS and investigate the entanglement properties of the resulting state. In addition, we demonstrate that the local observables can be evaluated efficiently. Finally, we introduce an infinite time-evolving block decimation algorithm (iTEBD\textsuperscript2) and use it to approximate the ground state of the 2D transverse field Ising model on lattices of infinite strip geometry.
Quantum convolutional neural networks (QCNNs) have been introduced as classifiers for gapped quantum phases of matter. Here, we propose a model-independent protocol for training QCNNs to discover order parameters that are unchanged under phase-preserving perturbations. We initiate the training sequence with the fixed-point wavefunctions of the quantum phase and then add translation-invariant noise that respects the symmetries of the system to mask the fixed-point structure on short length scales. We illustrate this approach by training the QCNN on phases protected by time-reversal symmetry in one dimension, and test it on several time-reversal symmetric models exhibiting trivial, symmetry-breaking, and symmetry-protected topological order. The QCNN discovers a set of order parameters that identifies all three phases and accurately predicts the location of the phase boundary. The proposed protocol paves the way towards hardware-efficient training of quantum phase classifiers on a programmable quantum processor.
Simulating time evolution of generic quantum many-body systems using classical numerical approaches has an exponentially growing cost either with evolution time or with the system size. In this work, we present a polynomially scaling hybrid quantum-classical algorithm for time evolving a one-dimensional uniform system in the thermodynamic limit. This algorithm uses a layered uniform sequential quantum circuit as a variational ansatz to represent infinite translation-invariant quantum states. We show numerically that this ansatz requires a number of parameters polynomial in the simulation time for a given accuracy. Furthermore, this favourable scaling of the ansatz is maintained during our variational evolution algorithm. All steps of the hybrid optimization are designed with near-term digital quantum computers in mind. After benchmarking the evolution algorithm on a classical computer, we demonstrate the measurement of observables of this uniform state using a finite number of qubits on a cloud-based quantum processing unit. With more efficient tensor contraction schemes, this algorithm may also offer improvements as a classical numerical algorithm.
The advent of noisy-intermediate scale quantum computers has introduced the exciting possibility of achieving quantum speedups in machine learning tasks. These devices, however, are composed of a small number of qubits, and can faithfully run only short circuits. This puts many proposed approaches for quantum machine learning beyond currently available devices. We address the problem of efficiently compressing and loading classical data for use on a quantum computer. Our proposed methods allow both the required number of qubits and depth of the quantum circuit to be tuned. We achieve this by using a correspondence between matrix-product states and quantum circuits, and further propose a hardware-efficient quantum circuit approach, which we benchmark on the Fashion-MNIST dataset. Finally, we demonstrate that a quantum circuit based classifier can achieve competitive accuracy with current tensor learning methods using only 11 qubits.
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.
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.
We investigate the computational power of the recently introduced class of isometric tensor network states (isoTNSs), which generalizes the isometric conditions of the canonical form of one-dimensional matrix-product states to tensor networks in higher dimensions. We discuss several technical details regarding the implementation of isoTNSs-based algorithms and compare different disentanglers -- which are essential for an efficient handling of isoTNSs. We then revisit the time evolving block decimation for isoTNSs ($\text{TEBD}^2$) and explore its power for real time evolution of two-dimensional (2D) lattice systems. Moreover, we introduce a density matrix renormalization group algorithm for isoTNSs ($\text{DMRG}^2$) that allows to variationally find ground states of 2D lattice systems. As a demonstration and benchmark, we compute the dynamical spin structure factor of 2D quantum spin systems for two paradigmatic models: First, we compare our results for the transverse field Ising model on a square lattice with the prediction of the spin-wave theory. Second, we consider the Kitaev model on the honeycomb lattice and compare it to the result from the exact solution.
The Bose-Hubbard model subjected to an effective magnetic field hosts a plethora of phases with different topological orders when tuning the chemical potential. Using the density matrix renormalization group method, we identify several gapped phases near the first Mott lobe at strong interactions. They are connected by a particle-hole symmetry to a variety of quantum Hall states stabilized at low fillings. We characterize phases of both particle and hole type and identify signatures compatible with Laughlin, Moore-Read, and Bosonic Integer Quantum Hall states by calculating the quantized Hall conductance and by extracting the topological entanglement entropy. Furthermore, we analyze the entanglement spectrum of a Laughlin state of bosonic particles and holes for a range of interaction strengths, as well as the entanglement spectrum of a Moore-Read state. These results further corroborate the existence of topological states at high fillings, close to the first Mott lobe, as hole analogues of the respective low-filling states.
Finding physical realizations of topologically ordered states in experimental settings, from condensed matter to artificial quantum systems, has been the main challenge en route to utilizing their unconventional properties. We show how to realize a large class of topologically ordered states and simulate their quasiparticle excitations on a digital quantum computer. To achieve this we design a set of linear-depth quantum circuits to generate ground states of general string-net models together with unitary open string operators to simulate the creation and braiding of abelian and non-abelian anyons. We show that the abelian (non-abelian) unitary string operators can be implemented with a constant (linear) depth quantum circuit. Our scheme allows us to directly probe characteristic topological properties, including topological entanglement entropy, braiding statistics, and fusion channels of anyons. Moreover, this set of efficiently prepared topologically ordered states has potential applications in the development of fault-tolerant quantum computers.
Quantum many-body systems may defy thermalization even without disorder. Intriguingly, non-ergodicity may be caused by a fragmentation of the many-body Hilbert-space into dynamically disconnected subspaces. The tilted one-dimensional Fermi-Hubbard model was proposed as a platform to realize fragmented models perturbatively in the limit of large tilt. Here, we demonstrate the validity of this effective description for the transient dynamics using ultracold fermions. The effective analytic model allows for a detailed understanding of the emergent microscopic processes, which in our case exhibit a pronounced doublon-number dependence. We study this experimentally by tuning the doublon fraction in the initial state.
Models whose ground states can be written as an exact matrix product state (MPS) provide valuable insights into phases of matter. While MPS-solvable models are typically studied as isolated points in a phase diagram, they can belong to a connected network of MPS-solvable models, which we call the MPS skeleton. As a case study where we can completely unearth this skeleton, we focus on the one-dimensional BDI class -- non-interacting spinless fermions with time-reversal symmetry. This class, labelled by a topological winding number, contains the Kitaev chain and is Jordan-Wigner-dual to various symmetry-breaking and symmetry-protected topological (SPT) spin chains. We show that one can read off from the Hamiltonian whether its ground state is an MPS: defining a polynomial whose coefficients are the Hamiltonian parameters, MPS-solvability corresponds to this polynomial being a perfect square. We provide an explicit construction of the ground state MPS, its bond dimension growing exponentially with the range of the Hamiltonian. This complete characterization of the MPS skeleton in parameter space has three significant consequences: (i) any two topologically distinct phases in this class admit a path of MPS-solvable models between them, including the phase transition which obeys an area law for its entanglement entropy; (ii) we illustrate that the subset of MPS-solvable models is dense in this class by constructing a sequence of MPS-solvable models which converge to the Kitaev chain (equivalently, the quantum Ising chain in a transverse field); (iii) a subset of these MPS states can be particularly efficiently processed on a noisy intermediate-scale quantum computer.
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.
K. J. Satzinger, Y. Liu, A. Smith, C. Knapp, M. Newman, C. Jones, Z. Chen, C. Quintana, X. Mi, A. Dunsworth, C. Gidney, I. Aleiner, F. Arute, K. Arya, J. Atalaya, R. Babbush, J. C. Bardin, R. Barends, J. Basso, A. Bengtsson, et al (78) The discovery of topological order has revolutionized the understanding of quantum matter in modern physics and provided the theoretical foundation for many quantum error correcting codes. Realizing topologically ordered states has proven to be extremely challenging in both condensed matter and synthetic quantum systems. Here, we prepare the ground state of the toric code Hamiltonian using an efficient quantum circuit on a superconducting quantum processor. We measure a topological entanglement entropy near the expected value of $\ln2$, and simulate anyon interferometry to extract the braiding statistics of the emergent excitations. Furthermore, we investigate key aspects of the surface code, including logical state injection and the decay of the non-local order parameter. Our results demonstrate the potential for quantum processors to provide key insights into topological quantum matter and quantum error correction.
Pimonpan Sompet, Sarah Hirthe, Dominik Bourgund, Thomas Chalopin, Julian Bibo, Joannis Koepsell, Petar Bojović, Ruben Verresen, Frank Pollmann, Guillaume Salomon, Christian Gross, Timon A. Hilker, Immanuel Bloch 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.
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.
In this work we study the many-body localization (MBL) transition and relate it to the eigenstate structure in the Fock space. Besides the standard entanglement and multifractal probes, we introduce the radial probability distribution of eigenstate coefficients with respect to the Hamming distance in the Fock space from the wave function maximum and relate the cumulants of this distribution to the properties of the quasi-local integrals of motion in the MBL phase. We demonstrate non-self-averaging property of the many-body fractal dimension $D_q$ and directly relate it to the jump of $D_q$ as well as of the localization length of the integrals of motion at the MBL transition. We provide an example of the continuous many-body transition confirming the above relation via the self-averaging of $D_q$ in the whole range of parameters. Introducing a simple toy-model, which hosts ergodic thermal bubbles, we give analytical evidences both in standard probes and in terms of newly introduced radial probability distribution that the MBL transition in the Fock space is consistent with the avalanche mechanism for delocalization, i.e., the Kosterlitz-Thouless scenario. Thus, we show that the MBL transition can been seen as a transition between ergodic states to non-ergodic extended states and put the upper bound for the disorder scaling for the genuine Anderson localization transition with respect to the non-interacting case.
In this work, we present a quantum information framework for the entanglement behavior of the low energy quasiparticle (QP) excitations in various quantum phases in one-dimensional (1D) systems. We first establish an exact correspondence between the correlation matrix and the QP entanglement Hamiltonian for free fermions and find an extended in-gap state in the QP entanglement Hamiltonian as a consequence of the position uncertainty of the QP. A more general understanding of such an in-gap state can be extended to a Kramers theorem for the QP entanglement Hamiltonian, which also applies to strongly interacting systems. Further, we present a set of ubiquitous entanglement spectrum features, dubbed entanglement fragmentation, conditional mutual information, and measurement induced non-local entanglement for QPs in 1D symmetry protected topological phases. Our result thus provides a new framework to identify different phases of matter in terms of their QP entanglement.
The thermalization of isolated quantum many-body systems is deeply related to fundamental questions of quantum information theory. While integrable or many-body localized systems display non-ergodic behavior due to extensively many conserved quantities, recent theoretical studies have identified a rich variety of more exotic phenomena in between these two extreme limits. The tilted one-dimensional Fermi-Hubbard model, which is readily accessible in experiments with ultracold atoms, emerged as an intriguing playground to study non-ergodic behavior in a clean disorder-free system. While non-ergodic behavior was established theoretically in certain limiting cases, there is no complete understanding of the complex thermalization properties of this model. In this work, we experimentally study the relaxation of an initial charge-density wave and find a remarkably long-lived initial-state memory over a wide range of parameters. Our observations are well reproduced by numerical simulations of a clean system. Using analytical calculations we further provide a detailed microscopic understanding of this behavior, which can be attributed to emergent kinetic constraints.
In a recent paper (Commun. Phys. 3, 100) Znidaric studies the growth of higher Renyi entropies in diffusive systems and claims that they generically grow ballistically in time, except for spin-1/2 models in d=1 dimension. Here, we point out that the necessary conditions for sub-ballistic growth of Renyi entropies are in fact much more general, and apply to a large class of systems, including experimentally relevant ones in arbitrary dimension and with larger local Hilbert spaces.
The current generation of noisy intermediate scale quantum computers introduces new opportunities to study quantum many-body systems. In this paper, we show that quantum circuits can provide a dramatically more efficient representation than current classical numerics of the quantum states generated under non-equilibrium quantum dynamics. For quantum circuits, we perform both real- and imaginary-time evolution using an optimization algorithm that is feasible on near-term quantum computers. We benchmark the algorithms by finding the ground state and simulating a global quench of the transverse field Ising model with a longitudinal field on a classical computer. Furthermore, we implement (classically optimized) gates on a quantum processing unit and demonstrate that our algorithm effectively captures real time evolution.
Analog quantum simulation has the potential to be an indispensable technique in the investigation of complex quantum systems. In this work, we numerically investigate a one-dimensional, faithful, analog, quantum electronic circuit simulator built out of Josephson junctions for one of the paradigmatic models of an integrable quantum field theory: the quantum sine-Gordon (qSG) model in 1+1 space-time dimensions. We analyze the lattice model using the density matrix renormalization group technique and benchmark our numerical results with existing Bethe ansatz computations. Furthermore, we perform analytical form-factor calculations for the two-point correlation function of vertex operators, which closely agree with our numerical computations. Finally, we compute the entanglement spectrum of the qSG model. We compare our results with those obtained using the integrable lattice-regularization based on the quantum XYZ chain and show that the quantum circuit model is less susceptible to corrections to scaling compared to the XYZ chain. We provide numerical evidence that the parameters required to realize the qSG model are accessible with modern-day superconducting circuit technology, thus providing additional credence towards the viability of the latter platform for simulating strongly interacting quantum field theories.
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.
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.
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.
The presence of global conserved quantities in interacting systems generically leads to diffusive transport at late times. Here, we show that systems conserving the dipole moment of an associated global charge, or even higher moment generalizations thereof, escape this scenario, displaying subdiffusive decay instead. Modelling the time evolution as cellular automata for specific cases of dipole- and quadrupole-conservation, we numerically find distinct anomalous exponents of the late time relaxation. We explain these findings by analytically constructing a general hydrodynamic model that results in a series of exponents depending on the number of conserved moments, yielding an accurate description of the scaling form of charge correlation functions. We analyze the spatial profile of the correlations and discuss potential experimentally relevant signatures of higher moment conservation.
Tensor networks permit computational and entanglement resources to be concentrated in interesting regions of Hilbert space. Implemented on NISQ machines they allow simulation of quantum systems that are much larger than the computational machine itself. This is achieved by parallelising the quantum simulation. Here, we demonstrate this in the simplest case; an infinite, translationally invariant quantum spin chain. We provide Cirq and Qiskit code that translate infinite, translationally invariant matrix product state (iMPS) algorithms to finite-depth quantum circuit machines, allowing the representation, optimisation and evolution arbitrary one-dimensional systems. Illustrative simulated output of these codes for achievable circuit sizes is given.
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 \emphundoped 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.
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.
Certain disorder-free Hamiltonians can be non-ergodic due to a \emphstrong 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 \emphstrong 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.
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.
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.
We study the growth of entanglement in quantum systems with a conserved quantity exhibiting diffusive transport, focusing on how initial inhomogeneities are imprinted on the entropy. We propose a simple effective model, which generalizes the minimal cut picture of \textitJonay et al. in such a way that the `line tension' of the cut depends on the local entropy density. In the case of noisy dynamics, this is described by a Kardar-Parisi-Zhang (KPZ) equation coupled to a diffusing field. We investigate the resulting dynamics and find that initial inhomogeneities of the conserved charge give rise to features in the entanglement profile, whose width and height both grow in time as $\propto\sqrt{t}$. In particular, for a domain wall quench, diffusion restricts entanglement growth to be $S_\text{vN} \lesssim \sqrt{t}$. We find that for charge density wave initial states, these features in the entanglement profile are present even after the charge density has equilibrated. Our conclusions are supported by numerical results on random circuits and deterministic spin chains.
Universal quantum computers are potentially an ideal setting for simulating many-body quantum dynamics that is out of reach for classical digital computers. We use state-of-the-art IBM quantum computers to study paradigmatic examples of condensed matter physics -- we simulate the effects of disorder and interactions on quantum particle transport, as well as correlation and entanglement spreading. Our benchmark results show that the quality of the current machines is below what is necessary for quantitatively accurate continuous time dynamics of observables and reachable system sizes are small comparable to exact diagonalization. Despite this, we are successfully able to demonstrate clear qualitative behaviour associated with localization physics and many-body interaction effects.
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.
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.
We explore the effect of local constraints on one-dimensional bosonic and fermionic ground state phases. Motivated by recent experiments on Rydberg chains, we constrain the occupation of neighboring sites in known phases of matter. Starting from Kitaev's topological superconductor wire, we find that a soft constraint induces a stable gapless Luttinger liquid phase. While Luttinger and Fermi liquids are usually unstable to superconducting proximity effects, the constraint suppresses pair creations, allowing for an emergent $U(1)$ symmetry and gaplessness. We substantiate this intuitive picture using field theoretical and Bethe ansatz methods. In particular, in the hard constraint limit, the model is explicitly $U(1)$-symmetric and integrable. For the corresponding spin-$1/2$ chains related by a Jordan-Wigner transformation, the Luttinger liquid is stabilized by the $\mathbb Z_2$ spin flip symmetry. Longer-range constraints stabilize gapless phases even without $\mathbb Z_2$ symmetry, connecting to the seminal work by Fendley, Sengupta and Sachdev [Physical Review B 69, 075106 (2004)], clarifying how the gapless floating phase observed therein can be vastly extended.
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.
Quantum states of matter---such as solids, magnets and topological phases---typically exhibit collective excitations---phonons, magnons, anyons. These involve the motion of many particles in the system, yet, remarkably, act like a single emergent entity---a quasiparticle. Known to be long-lived at the lowest energies, common wisdom says that quasiparticles become unstable when they encounter the inevitable continuum of many-particle excited states at high energies. Whilst correct for weak interactions, we show that this is far from the whole story: strong interactions generically stabilise quasiparticles by pushing them out of the continuum. This general mechanism is straightforwardly illustrated in an exactly solvable model. Using state-of-the-art numerics, we find it at work also in the spin-$\frac{1}{2}$ triangular lattice Heisenberg antiferromagnet (TLHAF) near the isotropic point---this is surprising given the common expectation of magnon decay in this paradigmatic frustrated magnet. Turning to existing experimental data, we identify the detailed phenomenology of avoided decay in the TLHAF material Ba$_3$CoSb$_2$O$_9$, and even in liquid helium---one of the earliest instances of quasiparticle decay. Our work unifies various phenomena above the universal low-energy regime in a comprehensive description. This broadens our window of understanding of many-body excitations, and provides a new perspective for controlling and stabilising quantum matter in the strongly-interacting regime.
We analyze a disordered central spin model, where a central spin interacts equally with each spin in a periodic one dimensional random-field Heisenberg chain. If the Heisenberg chain is initially in the many-body localized (MBL) phase, we find that the coupling to the central spin suffices to delocalize the chain for a substantial range of coupling strengths. We calculate the phase diagram of the model and identify the phase boundary between the MBL and ergodic phase. Within the localized phase, the central spin significantly enhances the rate of the logarithmic entanglement growth and its saturation value. We attribute the increase in entanglement entropy to a non-extensive enhancement of magnetization fluctuations induced by the central spin. Finally, we demonstrate that correlation functions of the central spin can be utilized to distinguish between MBL and ergodic phases of the 1D chain. Hence, we propose the use of a central spin as a possible experimental probe to identify the MBL phase.
The scrambling of quantum information in closed many-body systems, as measured by out-of-time-ordered correlation functions (OTOCs), has lately received considerable attention. Recently, a hydrodynamical description of OTOCs has emerged from considering random local circuits, aspects of which are conjectured to be universal to ergodic many-body systems, even without randomness. Here we extend this approach to systems with locally conserved quantities (e.g., energy). We do this by considering local random unitary circuits with a conserved U$(1)$ charge and argue, with numerical and analytical evidence, that the presence of a conservation law slows relaxation in both time ordered \textitand out-of-time-ordered correlation functions, both can have a diffusively relaxing component or "hydrodynamic tail" at late times. We verify the presence of such tails also in a deterministic, peridocially driven system. We show that for OTOCs, the combination of diffusive and ballistic components leads to a wave front with a specific, asymmetric shape, decaying as a power law behind the front. These results also explain existing numerical investigations in non-noisy ergodic systems with energy conservation. Moreover, we consider OTOCs in Gibbs states, parametrized by a chemical potential $\mu$, and apply perturbative arguments to show that for $\mu\gg 1$ the ballistic front of information-spreading can only develop at times exponentially large in $\mu$ -- with the information traveling diffusively at earlier times. We also develop a new formalism for describing OTOCs and operator spreading, which allows us to interpret the saturation of OTOCs as a form of thermalization on the Hilbert space of operators.
Thermalization and scrambling are the subject of much recent study from the perspective of many-body quantum systems with locally bounded Hilbert spaces (`spin chains'), quantum field theory and holography. We tackle this problem in 1D spin-chains evolving under random local unitary circuits and prove a number of exact results on the behavior of out-of-time-ordered commutators (OTOCs), and entanglement growth in this setting. These results follow from the observation that the spreading of operators in random circuits is described by a `hydrodynamical' equation of motion, despite the fact that random unitary circuits do not have locally conserved quantities (e.g., no conserved energy). In this hydrodynamic picture quantum information travels in a front with a `butterfly velocity' $v_{\text{B}}$ that is smaller than the light cone velocity of the system, while the front itself broadens diffusively in time. The OTOC increases sharply after the arrival of the light cone, but we do \emphnot observe a prolonged exponential regime of the form $\sim e^{\lambda_\text{L}(t-x/v)}$ for a fixed Lyapunov exponent $\lambda_\text{L}$. We find that the diffusive broadening of the front has important consequences for entanglement growth, leading to an entanglement velocity that can be significantly smaller than the butterfly velocity. We conjecture that the hydrodynamical description applies to more generic ergodic systems and support this by verifying numerically that the diffusive broadening of the operator wavefront also holds in a more traditional non-random Floquet spin-chain. We also compare our results to Clifford circuits, which have less rich hydrodynamics and consequently trivial OTOC behavior, but which can nevertheless exhibit linear entanglement growth and thermalization.
We study the dynamics of thermalization following a quantum quench using tensor-network methods. Contrary to the common belief that the rapid growth of entanglement and the resulting exponential growth of the bond dimension restricts simulations to short times, we demonstrate that the long time limit of local observables can be well captured using the time-dependent variational principle. This allows to extract transport coefficients such as the energy diffusion constant from simulations with rather small bond dimensions. We further study the characteristic of the chaotic wave that precedes the emergence of hydrodynamics, to find a ballistic diffusively-broadening wave-front.
We benchmark a variant of the recently introduced DMRG-X algorithm against exact results for the localized random field XX chain. We find that the eigenstates obtained via DMRG-X exhibit a highly accurate l-bit description for system sizes much bigger than the direct, many body, exact diagonalization in the spin variables is able to access. We take advantage of the underlying free fermion description of the XX model to accurately test the strengths and limitations of this algorithm for large system sizes. We discuss the theoretical constraints on the performance of the algorithm from the entanglement properties of the eigenstates, and its actual performance at different values of disorder. A small but significant improvement to the algorithm is also presented, which helps significantly with convergence. We find that at high entanglement, DMRG-X shows a bias towards eigenstates with low entanglement, but can be improved with increased bond dimension. This result suggests that one must be careful when applying the algorithm for interacting many body localized spin models near a transition.
We investigate the stationary and dynamical behavior of an Anderson localized chain coupled to a single central bound state. The coupling to the central site partially dilutes the Anderson localized peak towards the nearly resonant sites. In particular, the number of resonantly coupled sites remains finite in the thermodynamic limit. This is further supported by a multifractal analysis of eigenstates that shows the frozen spectrum of fractal dimension, which is characteristic for localized phases in models with power-law hopping. Although the well-known Fano-resonance problem is seemingly similar to our system, it fails to describe it because of the absence of level repulsion within the energy spectrum. For weak coupling strengths to the central site, we identify a regime with a logarithmic in time transport of particles and information.
We show that the anyonic statistics of fractionalized excitations display characteristic signatures in threshold spectroscopic measurements. Drawing motivation from topologically ordered phases such as gapped quantum spin liquids and fractional chern insulators which possess fractionalized excitations, we consider gapped systems with abelian anyonic excitations. The low energy onset of associated correlation functions near the threshold show universal behaviour depending on the statistics of the anyons. This explains some recent theoretical results in spin systems and also provides a route towards detecting statistics in experiments such as neutron scattering and tunneling spectroscopy.
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.
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.
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.