Source author record

Matthew B. Hastings

Matthew B. Hastings 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

36works
13topics
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

36 published item(s)

preprint2022arXiv

Adiabatic paths of Hamiltonians, symmetries of topological order, and automorphism codes

The recent "honeycomb code" is a fault-tolerant quantum memory defined by a sequence of checks which implements a nontrivial automorphism of the toric code. We argue that a general framework to understand this code is to consider continuous adiabatic paths of gapped Hamiltonians and we give a conjectured description of the fundamental group and second and third homotopy groups of this space in two spatial dimensions. A single cycle of such a path can implement some automorphism of the topological order of that Hamiltonian. We construct such paths for arbitrary automorphisms of two-dimensional doubled topological order. Then, realizing this in the case of the toric code, we turn this path back into a sequence of checks, constructing an automorphism code closely related to the honeycomb code.

preprint2022arXiv

How Dynamical Quantum Memories Forget

Motivated by recent work showing that a quantum error correcting code can be generated by hybrid dynamics of unitaries and measurements, we study the long time behavior of such systems. We demonstrate that even in the "mixed" phase, a maximally mixed initial density matrix is purified on a time scale equal to the Hilbert space dimension (i.e., exponential in system size), albeit with noisy dynamics at intermediate times which we connect to Dyson Brownian motion. In contrast, we show that free fermion systems -- i.e., ones where the unitaries are generated by quadratic Hamiltonians and the measurements are of fermion bilinears -- purify in a time quadratic in the system size. In particular, a volume law phase for the entanglement entropy cannot be sustained in a free fermion system.

preprint2022arXiv

On Lieb-Robinson Bounds for the Double Bracket Flow

We consider the possibility of developing a Lieb-Robinson bound for the double bracket flow. This is a differential equation $$\partial_B H(B)=[[V,H(B)],H(B)]$$ which may be used to diagonalize Hamiltonians. Here, $V$ is fixed and $H(0)=H$. We argue (but do not prove) that $H(B)$ need not converge to a limit for nonzero real $B$ in the infinite volume limit, even assuming several conditions on $H(0)$. However, we prove Lieb-Robinson bounds for all $B$ for the double-bracket flow for free fermion systems, but the range increases \emph{exponentially} with the control parameter $B$.

preprint2021arXiv

Measurement sequences for magic state distillation

Magic state distillation uses special codes to suppress errors in input states, which are often tailored to a Clifford-twirled error model. We present detailed measurement sequences for magic state distillation protocols which can suppress arbitrary errors on any part of a protocol, assuming the independence of errors across qubits. Provided with input magic states, our protocol operates on a two-dimensional square grid by measurements of $ZZ$ on horizontal pairs of qubits, $XX$ on vertical pairs, and $Z,X$ on single qubits.

preprint2021arXiv

Union-Find Decoders For Homological Product Codes

Homological product codes are a class of codes that can have improved distance while retaining relatively low stabilizer weight. We show how to build union-find decoders for these codes, using a union-find decoder for one of the codes in the product and a brute force decoder for the other code. We apply this construction to the specific case of the product of a surface code with a small code such as a $[[4,2,2]]$ code, which we call an augmented surface code. The distance of the augmented surface code is the product of the distance of the surface code with that of the small code, and the union-find decoder, with slight modifications, can decode errors up to half the distance. We present numerical simulations, showing that while the threshold of these augmented codes is lower than that of the surface code, the low noise performance is improved.

preprint2020arXiv

An Exactly Solvable Model for a $4+1D$ Beyond-Cohomology Symmetry Protected Topological Phase

We construct an exactly solvable commuting projector model for a $4+1$ dimensional ${\mathbb Z}_2$ symmetry-protected topological phase (SPT) which is outside the cohomology classification of SPTs. The model is described by a decorated domain wall construction, with "three-fermion" Walker-Wang phases on the domain walls. We describe the anomalous nature of the phase in several ways. One interesting feature is that, in contrast to in-cohomology phases, the effective ${\mathbb Z}_2$ symmetry on a $3+1$ dimensional boundary cannot be described by a quantum circuit and instead is a nontrivial quantum cellular automaton (QCA). A related property is that a codimension-two defect (for example, the termination of a ${\mathbb Z}_2$ domain wall at a trivial boundary) will carry nontrivial chiral central charge $4$ mod $8$. We also construct a gapped symmetric topologically-ordered boundary state for our model, which constitutes an anomalous symmetry enriched topological phase outside of the classification of arXiv:1602.00187, and define a corresponding anomaly indicator.

preprint2019arXiv

The Group Structure of Quantum Cellular Automata

We consider the group structure of quantum cellular automata (QCA) modulo circuits and show that it is abelian even without assuming the presence of ancillas, at least for most reasonable choices of control space; this is a corollary of a general method of ancilla removal. Further, we show how to define a group of QCA that is well-defined without needing to use families, by showing how to construct a coherent family containing an arbitrary finite QCA; the coherent family consists of QCA on progressively finer systems of qudits where any two members are related by a shallow quantum circuit. This construction applied to translation invariant QCA shows that all translation invariant QCA in three dimensions and all translation invariant Clifford QCA in any dimension are coherent.

preprint2016arXiv

Double semions in arbitrary dimension

We present a generalization of the double semion topological quantum field theory to higher dimensions, as a theory of $d-1$ dimensional surfaces in a $d$ dimensional ambient space. We construct a local Hamiltonian which is a sum of commuting projectors and analyze the excitations and the ground state degeneracy. Defining a consistent set of local rules requires the sign structure of the ground state wavefunction to depend not just on the number of disconnected surfaces, but also upon their higher Betti numbers through the semicharacteristic. For odd $d$ the theory is related to the toric code by a local unitary transformation, but for even $d$ the dimension of the space of zero energy ground states is in general different from the toric code and for even $d>2$ it is also in general different from that of the twisted $Z_2$ Dijkgraaf-Witten model.

preprint2016arXiv

Hybrid quantum-classical approach to correlated materials

Recent improvements in control of quantum systems make it seem feasible to finally build a quantum computer within a decade. While it has been shown that such a quantum computer can in principle solve certain small electronic structure problems and idealized model Hamiltonians, the highly relevant problem of directly solving a complex correlated material appears to require a prohibitive amount of resources. Here, we show that by using a hybrid quantum-classical algorithm that incorporates the power of a small quantum computer into a framework of classical embedding algorithms, the electronic structure of complex correlated materials can be efficiently tackled using a quantum computer. In our approach, the quantum computer solves a small effective quantum impurity problem that is self-consistently determined via a feedback loop between the quantum and classical computation. Use of a quantum computer enables much larger and more accurate simulations than with any known classical algorithm, and will allow many open questions in quantum materials to be resolved once a small quantum computer with around one hundred logical qubits becomes available.

preprint2016arXiv

Optimal Circuit-Level Decoding for Surface Codes

Surface codes exploit topological protection to increase error resilience in quantum computing devices and can in principle be implemented in existing hardware. They are one of the most promising candidates for active error correction, not least due to a polynomial-time decoding algorithm which admits one of the highest predicted error thresholds. We consider the dependence of this threshold on underlying assumptions including different noise models, and analyze the performance of a minimum weight perfect matching (MWPM) decoding compared to a mathematically optimal maximum likelihood (ML) decoding. Our ML algorithm tracks the success probabilities for all possible corrections over time and accounts for individual gate failure probabilities and error propagation due to the syndrome measurement circuit. We present the very first evidence for the true error threshold of an optimal circuit level decoder, allowing us to draw conclusions about what kind of improvements are possible over standard MWPM.

preprint2015arXiv

From local to global ground states in Ising spin glasses

We consider whether it is possible to find ground states of frustrated spin systems by solving them locally. Using spin glass physics and Imry-Ma arguments in addition to numerical benchmarks we quantify the power of such local solution methods and show that for the average low-dimensional spin glass problem outside the spin- glass phase the exact ground state can be found in polynomial time. In the second part we present a heuristic, general-purpose hierarchical approach which for spin glasses on chimera graphs and lattices in two and three dimensions outperforms, to our knowledge, any other solver currently around, with significantly better scaling performance than simulated annealing.

preprint2015arXiv

Nonzero temperature Entanglement Negativity of quantum spin models: Area law, Linked Cluster Expansions and Sudden Death

We show that the bipartite logarithmic entanglement negativity (EN) of quantum spin models obeys an area law at all nonzero temperatures. We develop numerical linked cluster (NLC) expansions for the `area-law' logarithmic entanglement negativity as a function of temperature and other parameters. For one-dimensional models the results of NLC are compared with exact diagonalization on finite systems and are found to agree very well. The NLC results are also obtained for two dimensional XXZ and transverse-field Ising models. In all cases, we find a sudden onset (or sudden death) of negativity at a finite temperature above which the negativity is zero. We use perturbation theory to develop a physical picture for this sudden onset (or sudden death). The onset of EN or its magnitude are insensitive to classical finite-temperature phase transitions, supporting the argument for absence of any role of quantum mechanics at such transitions. On approach to a quantum critical point at $T=0$, negativity shows critical scaling in size and temperature.

preprint2015arXiv

Slowest local operators in quantum spin chains

We numerically construct slowly relaxing local operators in a nonintegrable spin-1/2 chain. Restricting the support of the operator to $M$ consecutive spins along the chain, we exhaustively search for the operator that minimizes the Frobenius norm of the commutator with the Hamiltonian. We first show that the Frobenius norm bounds the time scale of relaxation of the operator at high temperatures. We find operators with significantly slower relaxation than the slowest simple "hydrodynamic" mode due to energy diffusion. Then, we examine some properties of the nontrivial slow operators. Using both exhaustive search and tensor network techniques, we find similar slowly relaxing operators for a Floquet spin chain; this system is hydrodynamically "trivial", with no conservation laws restricting their dynamics. We argue that such slow relaxation may be a generic feature following from locality and unitarity.

preprint2015arXiv

Solving strongly correlated electron models on a quantum computer

One of the main applications of future quantum computers will be the simulation of quantum models. While the evolution of a quantum state under a Hamiltonian is straightforward (if sometimes expensive), using quantum computers to determine the ground state phase diagram of a quantum model and the properties of its phases is more involved. Using the Hubbard model as a prototypical example, we here show all the steps necessary to determine its phase diagram and ground state properties on a quantum computer. In particular, we discuss strategies for efficiently determining and preparing the ground state of the Hubbard model starting from various mean-field states with broken symmetry. We present an efficient procedure to prepare arbitrary Slater determinants as initial states and present the complete set of quantum circuits needed to evolve from these to the ground state of the Hubbard model. We show that, using efficient nesting of the various terms each time step in the evolution can be performed with just $\mathcal{O}(N)$ gates and $\mathcal{O}(\log N)$ circuit depth. We give explicit circuits to measure arbitrary local observables and static and dynamic correlation functions, both in the time and frequency domain. We further present efficient non-destructive approaches to measurement that avoid the need to re-prepare the ground state after each measurement and that quadratically reduce the measurement error.

preprint2014arXiv

Gate count estimates for performing quantum chemistry on small quantum computers

As quantum computing technology improves and quantum computers with a small but non-trivial number of N > 100 qubits appear feasible in the near future the question of possible applications of small quantum computers gains importance. One frequently mentioned application is Feynman's original proposal of simulating quantum systems, and in particular the electronic structure of molecules and materials. In this paper, we analyze the computational requirements for one of the standard algorithms to perform quantum chemistry on a quantum computer. We focus on the quantum resources required to find the ground state of a molecule twice as large as what current classical computers can solve exactly. We find that while such a problem requires about a ten-fold increase in the number of qubits over current technology, the required increase in the number of gates that can be coherently executed is many orders of magnitude larger. This suggests that for quantum computation to become useful for quantum chemistry problems, drastic algorithmic improvements will be needed.

preprint2014arXiv

On Metaplectic Modular Categories and their applications

For non-abelian simple objects in a unitary modular category, the density of their braid group representations, the #P-hard evaluation of their associated link invariants, and the BQP-completeness of their anyonic quantum computing models are closely related. We systematically study such properties of the non-abelian simple objects in metaplectic modular categories, which are unitary modular categories with fusion rules of SO(m)_2 for an odd integer m \geq 3. The simple objects with quantum dimensions \sqrt{m} have finite image braid group representations, and their link invariants are classically efficient to evaluate. We also provide classically efficient simulation of their braid group representations. These simulations of the braid group representations can be regarded as qudit generalizations of the Knill-Gottesmann theorem for the qubit case. The simple objects of dimension 2 give us a surprising result: while their braid group representations have finite images and are efficiently simulable classically after a generalized localization, their link invariants are #P-hard to evaluate exactly. We sharpen the #P-hardness by showing that any sufficiently accurate approximation of their associated link invariants is already #P-hard.

preprint2014arXiv

Quantization of Hall Conductance For Interacting Electrons on a Torus

We consider interacting, charged spins on a torus described by a gapped Hamiltonian with a unique groundstate and conserved local charge. Using quasi-adiabatic evolution of the groundstate around a flux-torus, we prove, without any averaging assumption, that the Hall conductance of the groundstate is quantized in integer multiples of e^2/h, up to exponentially small corrections in the linear size of the system. In addition, we discuss extensions to the fractional quantization case under an additional topological order assumption on the degenerate groundstate subspace.

preprint2013arXiv

Anomalies in the Entanglement Properties of the Square Lattice Heisenberg Model

We compute the bipartite entanglement properties of the spin-half square-lattice Heisenberg model by a variety of numerical techniques that include valence bond quantum Monte Carlo (QMC), stochastic series expansion QMC, high temperature series expansions and zero temperature coupling constant expansions around the Ising limit. We find that the area law is always satisfied, but in addition to the entanglement entropy per unit boundary length, there are other terms that depend logarithmically on the subregion size, arising from broken symmetry in the bulk and from the existence of corners at the boundary. We find that the numerical results are anomalous in several ways. First, the bulk term arising from broken symmetry deviates from an exact calculation that can be done for a mean-field Neel state. Second, the corner logs do not agree with the known results for non-interacting Boson modes. And, third, even the finite temperature mutual information shows an anomalous behavior as T goes to zero, suggesting that T->0 and L->infinity limits do not commute. These calculations show that entanglement entropy demonstrates a very rich behavior in d>1, which deserves further attention.

preprint2013arXiv

Faster Phase Estimation

We develop several algorithms for performing quantum phase estimation based on basic measurements and classical post-processing. We present a pedagogical review of quantum phase estimation and simulate the algorithm to numerically determine its scaling in circuit depth and width. We show that the use of purely random measurements requires a number of measurements that is optimal up to constant factors, albeit at the cost of exponential classical post-processing; the method can also be used to improve classical signal processing. We then develop a quantum algorithm for phase estimation that yields an asymptotic improvement in runtime, coming within a factor of log* of the minimum number of measurements required while still requiring only minimal classical post-processing. The corresponding quantum circuit requires asymptotically lower depth and width (number of qubits) than quantum phase estimation.

preprint2013arXiv

Homological Product Codes

Quantum codes with low-weight stabilizers known as LDPC codes have been actively studied recently due to their simple syndrome readout circuits and potential applications in fault-tolerant quantum computing. However, all families of quantum LDPC codes known to this date suffer from a poor distance scaling limited by the square-root of the code length. This is in a sharp contrast with the classical case where good families of LDPC codes are known that combine constant encoding rate and linear distance. Here we propose the first family of good quantum codes with low-weight stabilizers. The new codes have a constant encoding rate, linear distance, and stabilizers acting on at most $\sqrt{n}$ qubits, where $n$ is the code length. For comparison, all previously known families of good quantum codes have stabilizers of linear weight. Our proof combines two techniques: randomized constructions of good quantum codes and the homological product operation from algebraic topology. We conjecture that similar methods can produce good stabilizer codes with stabilizer weight $n^a$ for any $a>0$. Finally, we apply the homological product to construct new small codes with low-weight stabilizers.

preprint2013arXiv

Self-Correcting Quantum Memories Beyond the Percolation Threshold

We analyze several high dimensional generalizations of the toric code at nonzero temperature. We find that in large enough dimension, there can be a distinct separation between the critical temperature $T_c$, given by thermodynamic singularities, and the percolation temperature $T_p$, given by the percolation of defects. We argue that the regime $T_p<T<T_c$ is a range of temperatures where a self-correcting quantum memory can operate despite having percolating defects. We present analytic arguments and numerical evidence in support of this scenario, including a mean-field treatment and Monte Carlo simulations. Near $T_c$, simulations observe a large hysteretic behavior, which may have applications by allowing the self-correcting phase to survive in a "superheated" regime.

preprint2012arXiv

Effect of thermal fluctuations in topological p-wave superconductors

We study the effect of thermal fluctuations on the topological stability of chiral p-wave superconductors. We consider two models of superconductors: spinless and spinful with a focus on topological properties and Majorana zero-energy modes. We show that proliferation of vortex-antivortex pairs above the Kosterlitz-Thouless temperature T_KT drives the transition from a thermal Quantum Hall insulator to a thermal metal/insulator, and dramatically modifies the ground-state degeneracy splitting. Therefore, in order to utilize 2D chiral p-wave superconductors for topological quantum computing, the temperature should be much smaller than T_KT. Within the spinful chiral p-wave model, we also investigate the interplay between half-quantum vortices carrying Majorana zero-energy modes and full-quantum vortices having trivial topological charge, and discuss topological properties of half-quantum vortices in the background of proliferating full-quantum vortices.

preprint2012arXiv

Entanglement scaling in two-dimensional gapless systems

We numerically determine subleading scaling terms in the ground-state entanglement entropy of several two-dimensional (2D) gapless systems, including a Heisenberg model with Néel order, a free Dirac fermion in the π-flux phase, and the nearest-neighbor resonating-valence-bond wavefunction. For these models, we show that the entanglement entropy between cylindrical regions of length x and L - x, extending around a torus of length L, depends upon the dimensionless ratio x/L. This can be well-approximated on finite-size lattices by a function ln(sin(πx/L)), akin to the familiar chord-length dependence in one dimension. We provide evidence, however, that the precise form of this bulk-dependent contribution is a more general function in the 2D thermodynamic limit.

preprint2012arXiv

Metaplectic Anyons, Majorana Zero Modes, and their Computational Power

We introduce and study a class of anyon models that are a natural generalization of Ising anyons and Majorana fermion zero modes. These models combine an Ising anyon sector with a sector associated with $SO(m)_2$ Chern-Simons theory. We show how they can arise in a simple scenario for electron fractionalization and give a complete account of their quasiparticles types, fusion rules, and braiding. We show that the image of the braid group is finite for a collection of $2n$ fundamental quasiparticles and is a proper subgroup of the metaplectic representation of $Sp(2n-2,\mathbb{F}_m)\ltimes H(2n-2,\mathbb{F}_m)$, where $Sp(2n-2,\mathbb{F}_m)$ is the symplectic group over the finite field $\mathbb{F}_m$ and $H(2n-2,\mathbb{F}_m)$ is the extra special group (also called the $(2n-1)$-dimensional Heisenberg group) over $\mathbb{F}_m$. Moreover, the braiding of fundamental quasiparticles can be efficiently simulated classically. However, computing the result of braiding a certain type of composite quasiparticle is $# P$-hard, although it is not universal for quantum computation because it has a finite braid group image. This a rare example of a topological phase that is not universal for quantum computation through braiding but nevertheless has $# P$-hard link invariants. We argue that our models are closely related to recent analyses finding non-Abelian anyonic properties for defects in quantum Hall systems, generalizing Majorana zero modes in quasi-1D systems.

preprint2011arXiv

Finite Temperature Critical Behavior of Mutual Information

We study mutual information for Renyi entropy of arbitrary index n, in interacting quantum systems at finite-temperature critical points, using high-temperature expansion, quantum Monte Carlo simulations and scaling theory. We find that for n>1, the critical behavior is manifest at two temperatures T_c and n*T_c. For the XXZ model with Ising anisotropy, the coefficient of the area-law has a t*ln(t) singularity, whereas the subleading correction from corners has a logarithmic divergence, with a coefficient related to the exact results of Cardy and Peschel. For T<n*T_c there is a constant term associated with broken symmetries that jumps at both T_c and n*T_c, which can be understood in terms of a scaling function analogous to the boundary entropy of Affleck and Ludwig.

preprint2011arXiv

Galois Conjugates of Topological Phases

Galois conjugation relates unitary conformal field theories (CFTs) and topological quantum field theories (TQFTs) to their non-unitary counterparts. Here we investigate Galois conjugates of quantum double models, such as the Levin-Wen model. While these Galois conjugated Hamiltonians are typically non-Hermitian, we find that their ground state wave functions still obey a generalized version of the usual code property (local operators do not act on the ground state manifold) and hence enjoy a generalized topological protection. The key question addressed in this paper is whether such non-unitary topological phases can also appear as the ground states of Hermitian Hamiltonians. Specific attempts at constructing Hermitian Hamiltonians with these ground states lead to a loss of the code property and topological protection of the degenerate ground states. Beyond this we rigorously prove that no local change of basis can transform the ground states of the Galois conjugated doubled Fibonacci theory into the ground states of a topological model whose Hermitian Hamiltonian satisfies Lieb-Robinson bounds. These include all gapped local or quasi-local Hamiltonians. A similar statement holds for many other non-unitary TQFTs. One consequence is that the "Gaffnian" wave function cannot be the ground state of a gapped fractional quantum Hall state.

preprint2011arXiv

Topological Entanglement Entropy of a Bose-Hubbard Spin Liquid

The Landau paradigm of classifying phases by broken symmetries was demonstrated to be incomplete when it was realized that different quantum Hall states could only be distinguished by more subtle, topological properties. Today, the role of topology as an underlying description of order has branched out to include topological band insulators, and certain featureless gapped Mott insulators with a topological degeneracy in the groundstate wavefunction. Despite intense focus, very few candidates for these topologically ordered "spin liquids" exist. The main difficulty in finding systems that harbour spin liquid states is the very fact that they violate the Landau paradigm, making conventional order parameters non-existent. Here, we uncover a spin liquid phase in a Bose-Hubbard model on the kagome lattice, and measure its topological order directly via the topological entanglement entropy. This is the first smoking-gun demonstration of a non-trivial spin liquid, identified through its entanglement entropy as a gapped groundstate with emergent Z2 gauge symmetry.

preprint2011arXiv

Universal Signatures of Fractionalized Quantum Critical Points

Groundstates of certain materials can support exotic excitations with a charge that's a fraction of the fundamental electron charge. The condensation of these fractionalized particles has been predicted to drive novel quantum phase transitions, which haven't yet been observed in realistic systems. Through numerical and theoretical analysis of a physical model of interacting lattice bosons, we establish the existence of such an exotic critical point, called XY*. We measure a highly non-classical critical exponent eta = 1.49(2), and construct a universal scaling function of winding number distributions that directly demonstrates the distinct topological sectors of an emergent Z_2 gauge field. The universal quantities used to establish this exotic transition can be used to detect other fractionalized quantum critical points in future model and material systems.

preprint2011arXiv

Weakly-Coupled non-Abelian Anyons in Three Dimensions

We introduce a Hamiltonian coupling Majorana fermion degrees of freedom to a quantum dimer model. We argue that, in three dimensions, this model has deconfined quasiparticles supporting Majorana zero modes obeying nontrivial statistics. We introduce two effective field theory descriptions of this deconfined phase, in which the excitations have Coulomb interactions. A key feature of this system is the existence of topologically non-trivial fermionic excitations, called "Hopfions" because, in a suitable continuum limit of the dimer model, such excitations correspond to the Hopf map and are related to excitations identified in arXiv:1003.1964. We identify corresponding topological invariants of the quantum dimer model (with or without fermions) which are present even on lattices with trivial topology. The Hopfion energy gap depends upon the phase of the model. We briefly comment on the possibility of a phase with a gapped, deconfined $\mathbb{Z}_2$ gauge field, as may arise on the stacked triangular lattice.

preprint2010arXiv

Finite Size Scaling of Mutual Information: A Scalable Simulation

We develop a quantum Monte Carlo procedure to compute the Renyi mutual information of an interacting quantum many-body system at non-zero temperature. Performing simulations on a spin-1/2 XXZ model, we observe that for a subregion of fixed size embedded in a system of size L, the mutual information converges at large L to a limiting function which displays non-monotonic temperature behavior corresponding to the onset of correlations. For a region of size L/2 embedded in a system of size L, the mutual information divided by L converges to a limiting function of temperature, with apparently nontrivial corrections near critical points.

preprint2010arXiv

Markov entropy decomposition: a variational dual for quantum belief propagation

We present a lower bound for the free energy of a quantum many-body system at finite temperature. This lower bound is expressed as a convex optimization problem with linear constraints, and is derived using strong subadditivity of von Neumann entropy and a relaxation of the consistency condition of local density operators. The dual to this minimization problem leads to a set of quantum belief propagation equations, thus providing a firm theoretical foundation to that approach. The minimization problem is numerically tractable, and we find good agreement with quantum Monte Carlo for the spin-half Heisenberg anti-ferromagnet in two dimensions. This lower bound complements other variational upper bounds. We discuss applications to Hamiltonian complexity theory and give a generalization of the structure theorem of Hayden, Jozsa, Petz and Winter to trees in an appendix.

preprint2010arXiv

Measuring Renyi Entanglement Entropy with Quantum Monte Carlo

We develop a quantum Monte Carlo procedure, in the valence bond basis, to measure the Renyi entanglement entropy of a many-body ground state as the expectation value of a unitary {\it Swap} operator acting on two copies of the system. An improved estimator involving the ratio of {\it Swap} operators for different subregions enables simulations to converge the entropy in a time polynomial in the system size. We demonstrate convergence of the Renyi entropy to exact results for a Heisenberg chain. Finally, we calculate the scaling of the Renyi entropy in the two-dimensional Heisenberg model and confirm that the Néel groundstate obeys the expected area law for systems up to linear size L=28.

preprint2010arXiv

Projective Ribbon Permutation Statistics: a Remnant of non-Abelian Braiding in Higher Dimensions

In a recent paper, Teo and Kane proposed a 3D model in which the defects support Majorana fermion zero modes. They argued that exchanging and twisting these defects would implement a set R of unitary transformations on the zero mode Hilbert space which is a 'ghostly' recollection of the action of the braid group on Ising anyons in 2D. In this paper, we find the group T_{2n} which governs the statistics of these defects by analyzing the topology of the space K_{2n} of configurations of 2n defects in a slowly spatially-varying gapped free fermion Hamiltonian: T_{2n}\equiv {π_1}(K_{2n})$. We find that the group T_{2n}= Z \times T^r_{2n}, where the 'ribbon permutation group' T^r_{2n} is a mild enhancement of the permutation group S_{2n}: T^r_{2n} \equiv \Z_2 \times E((Z_2)^{2n}\rtimes S_{2n}). Here, E((Z_2)^{2n}\rtimes S_{2n}) is the 'even part' of (Z_2)^{2n} \rtimes S_{2n}, namely those elements for which the total parity of the element in (Z_2)^{2n} added to the parity of the permutation is even. Surprisingly, R is only a projective representation of T_{2n}, a possibility proposed by Wilczek. Thus, Teo and Kane's defects realize `Projective Ribbon Permutation Statistics', which we show to be consistent with locality. We extend this phenomenon to other dimensions, co-dimensions, and symmetry classes. Since it is an essential input for our calculation, we review the topological classification of gapped free fermion systems and its relation to Bott periodicity.

preprint2010arXiv

Random tensor theory: extending random matrix theory to random product states

We consider a problem in random matrix theory that is inspired by quantum information theory: determining the largest eigenvalue of a sum of p random product states in (C^d)^{otimes k}, where k and p/d^k are fixed while d grows. When k=1, the Marcenko-Pastur law determines (up to small corrections) not only the largest eigenvalue ((1+sqrt{p/d^k})^2) but the smallest eigenvalue (min(0,1-sqrt{p/d^k})^2) and the spectral density in between. We use the method of moments to show that for k>1 the largest eigenvalue is still approximately (1+sqrt{p/d^k})^2 and the spectral density approaches that of the Marcenko-Pastur law, generalizing the random matrix theory result to the random tensor case. Our bound on the largest eigenvalue has implications both for sampling from a particular heavy-tailed distribution and for a recently proposed quantum data-hiding and correlation-locking scheme due to Leung and Winter. Since the matrices we consider have neither independent entries nor unitary invariance, we need to develop new techniques for their analysis. The main contribution of this paper is to give three different methods for analyzing mixtures of random product states: a diagrammatic approach based on Gaussian integrals, a combinatorial method that looks at the cycle decompositions of permutations and a recursive method that uses a variant of the Schwinger-Dyson equations.

preprint2010arXiv

Strong and weak thermalization of infinite non-integrable quantum systems

When a non-integrable system evolves out of equilibrium for a long time, local observables are expected to attain stationary expectation values, independent of the details of the initial state. However, intriguing experimental results with ultracold gases have shown no thermalization in non-integrable settings, triggering an intense theoretical effort to decide the question. Here we show that the phenomenology of thermalization in a quantum system is much richer than its classical counterpart. Using a new numerical technique, we identify two distinct thermalization regimes, strong and weak, occurring for different initial states. Strong thermalization, intrinsically quantum, happens when instantaneous local expectation values converge to the thermal ones. Weak thermalization, well-known in classical systems, happens when local expectation values converge to the thermal ones only after time averaging. Remarkably, we find a third group of states showing no thermalization, neither strong nor weak, to the time scales one can reliably simulate.