Source author record

M. B. Hastings

M. 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

39works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

39 published item(s)

preprint2020arXiv

Classical and Quantum Algorithms for Tensor Principal Component Analysis

We present classical and quantum algorithms based on spectral methods for a problem in tensor principal component analysis. The quantum algorithm achieves a quartic speedup while using exponentially smaller space than the fastest classical spectral algorithm, and a super-polynomial speedup over classical algorithms that use only polynomial space. The classical algorithms that we present are related to, but slightly different from those presented recently in Ref. 1. In particular, we have an improved threshold for recovery and the algorithms we present work for both even and odd order tensors. These results suggest that large-scale inference problems are a promising future application for quantum computers.

preprint2020arXiv

How Quantum Are Non-Negative Wavefunctions?

We consider wavefunctions which are non-negative in some tensor product basis. We study what possible teleportation can occur in such wavefunctions, giving a complete answer in some cases (when one system is a qubit) and partial answers elsewhere. We use this to show that a one-dimensional wavefunction which is non-negative and has zero correlation length can be written in a "coherent Gibbs state" form, as explained later. We conjecture that such holds in higher dimensions. Additionally, some results are provided on possible teleportation in general wavefunctions, explaining how Schmidt coefficients before measurement limit the possible Schmidt coefficients after measurement, and on the absence of a "generalized area law"\cite{genarealaw} even for Hamiltonians with no sign problem. One of the motivations for this work is an attempt to prove a conjecture about ground state wavefunctions which have an "intrinsic" sign problem that cannot be removed by any quantum circuit. We show a weaker version of this, showing that the sign problem is intrinsic for commuting Hamiltonians in the same phase as the double semion model under the technical assumption that TQO-2 holds\cite{tqo2}.

preprint2016arXiv

Local Maxima and Improved Exact Algorithm for MAX-2-SAT

Given a MAX-2-SAT instance, we define a local maximum to be an assignment such that changing any single variable reduces the number of satisfied clauses. We consider the question of the number of local maxima that an instance of MAX-2-SAT can have. We give upper bounds in both the sparse and nonsparse case, where the sparse case means that there is a bound $d$ on the average number of clauses involving any given variable. The bounds in the nonsparse case are tight up to polylogarithmic factors, while in the sparse case the bounds are tight up to a multiplicative factor in $d$ for large $d$. Additionally, we generalize to the question of assignments which are maxima up to changing $k> 1$ variables simultaneously; in this case, we give explicit constructions with large (in a sense explained below) numbers of such maxima in the sparse case. The basic idea of the upper bound proof is to consider a random assignment to some subset of the variables and determine the probability that some fraction of the remaining variables can be fixed without considering interactions between them. The bounded results hold in the case of weighted MAX-2-SAT as well. Using this technique and combining with ideas from Ref. 6, we find an algorithm for weighted MAX-2-SAT which is faster for large $d$ than previous algorithms which use polynomial space; this algorithm does require an additional bounds on maximum weights and degree.

preprint2016arXiv

Quantum Codes from High-Dimensional Manifolds

We construct toric codes on various high-dimensional manifolds. Assuming a conjecture in geometry we find families of quantum CSS stabilizer codes on $N$ qubits with logarithmic weight stabilizers and distance $N^{1-ε}$ for any $ε>0$. The conjecture is that there is a constant $C>0$ such that for any $n$-dimensional torus ${\mathbb T}^n={\mathbb R}^n/Λ$, where $Λ$ is a lattice, the least volume unoriented $n/2$-dimensional surface (using the Euclidean metric) representing nontrivial homology has volume at least $C^n$ times the volume of the least volume $n/2$-dimensional hyperplane representing nontrivial homology; in fact, it would suffice to have this result for $Λ$ an integral lattice with the surface restricted to faces of a cubulation by unit hypercubes. The main technical result is an estimate of Rankin invariants\cite{rankin} for certain random lattices, showing that in a certain sense they are optimal. Additionally, we construct codes with square-root distance, logarithmic weight stabilizers, and inverse polylogarithmic soundness factor (considered as quantum locally testable codes\cite{qltc}). We also provide an short, alternative proof that the shortest vector in the exterior power of a lattice may be non-split\cite{coulangeon}.

preprint2016arXiv

Random MERA States and the Tightness of the Brandao-Horodecki Entropy Bound

We construct a random MERA state with a bond dimension that varies with the level of the MERA. This causes the state to exhibit a very different entanglement structure from that usually seen in MERA, with neighboring intervals of length $l$ exhibiting a mutual information proportional to $εl$ for some constant $ε$, up to a length scale exponentially large in $ε$. We express the entropy of a random MERA in terms of sums over cuts through the MERA network, with the entropy in this case controlled by the cut minimizing bond dimensions cut through. One motivation for this construction is to investigate the tightness of the Brandao-Horodecki\cite{bh} entropy bound relating entanglement to correlation decay. Using the random MERA, we show that at least part of the proof is tight: there do exist states with the required property of having linear mutual information between neighboring intervals at all length scales. We conjecture that this state has exponential correlation decay and that it demonstrates that the Brandao-Horodecki bound is tight (at least up to constant factors), and we provide some numerical evidence for this as well as a sketch of how a proof of correlation decay might proceed.

preprint2016arXiv

The Asymptotics of Quantum Max-Flow Min-Cut

The quantum max-flow min-cut conjecture relates the rank of a tensor network to the minimum cut in the case that all tensors in the network are identical\cite{mfmc1}. This conjecture was shown to be false in Ref. \onlinecite{mfmc2} by an explicit counter-example. Here, we show that the conjecture is almost true, in that the ratio of the quantum max-flow to the quantum min-cut converges to $1$ as the dimension $N$ of the degrees of freedom on the edges of the network tends to infinity. The proof is based on estimating moments of the singular values of the network. We introduce a generalization of "rainbow diagrams"\cite{rainbow} to tensor networks to estimate the dominant diagrams. A direct comparison of second and fourth moments lower bounds the ratio of the quantum max-flow to the quantum min-cut by a constant. To show the tighter bound that the ratio tends to $1$, we consider higher moments. In addition, we show that the limiting moments as $N \rightarrow \infty$ agree with that in a different ensemble where tensors in the network are chosen independently, this is used to show that the distributions of singular values in the two different ensembles weakly converge to the same limiting distribution. We present also a numerical study of one particular tensor network, which shows a surprising dependence of the rank deficit on $N \mod 4$ and suggests further conjecture on the limiting behavior of the rank.

preprint2016arXiv

Training A Quantum Optimizer

We study a variant of the quantum approximate optimization algorithm [ E. Farhi, J. Goldstone, and S. Gutmann, arXiv:1411.4028] with slightly different parametrization and different objective: rather than looking for a state which approximately solves an optimization problem, our goal is to find a quantum algorithm that, given an instance of MAX-2-SAT, will produce a state with high overlap with the optimal state. Using a machine learning approach, we chose a "training set" of instances and optimized the parameters to produce large overlap for the training set. We then tested these optimized parameters on a larger instance set. As a training set, we used a subset of the hard instances studied by E. Crosson, E. Farhi, C. Yen-Yu Lin, H.-H. Lin, and P. Shor (CFLLS) [arXiv:1401.7320]. When tested on the full set, the parameters that we find produce significantly larger overlap than the optimized annealing times of CFLLS. Testing on other random instances from $20$ to $28$ bits continues to show improvement over annealing, with the improvement being most notable on the hardest instances. Further tests on instances of MAX-3-SAT also showed improvement on the hardest instances. This algorithm may be a possible application for near-term quantum computers with limited coherence times.

preprint2016arXiv

Turning Gate Synthesis Errors into Incoherent Errors

Using error correcting codes and fault tolerant techniques, it is possible, at least in theory, to produce logical qubits with significantly lower error rates than the underlying physical qubits. Suppose, however, that the gates that act on these logical qubits are only approximation of the desired gate. This can arise, for example, in synthesizing a single qubit unitary from a set of Clifford and $T$ gates; for a generic such unitary, any finite sequence of gates only approximates the desired target. In this case, errors in the gate can add coherently so that, roughly, the error $ε$ in the unitary of each gate must scale as $ε\lesssim 1/N$, where $N$ is the number of gates. If, however, one has the option of synthesizing one of several unitaries near the desired target, and if an average of these options is closer to the target, we give some elementary bounds showing cases in which the errors can be made to add incoherently by averaging over random choices, so that, roughly, one needs $ε\lesssim 1/\sqrt{N}$. We remark on one particular application to distilling magic states where this effect happens automatically in the usual circuits.

preprint2016arXiv

Weight Reduction for Quantum Codes

We present an algorithm that takes a CSS stabilizer code as input, and outputs another CSS stabilizer code such that the stabilizer generators all have weights $O(1)$ and such that $O(1)$ generators act on any given qubit. The number of logical qubits is unchanged by the procedure, while we give bounds on the increase in number of physical qubits and in the effect on distance and other code parameters, such as soundness (as a locally testable code) and "cosoundness" (defined later). Applications are discussed, including to codes from high-dimensional manifolds which have logarithmic weight stabilizers. Assuming a conjecture in geometry\cite{hdm}, this allows the construction of CSS stabilizer codes with generator weight $O(1)$ and almost linear distance. Another application of the construction is to increasing the distance to $X$ or $Z$ errors, whichever is smaller, so that the two distances are equal.

preprint2015arXiv

Reduced Space-Time and Time Costs Using Dislocation Codes and Arbitrary Ancillas

We propose two distinct methods of improving quantum computing protocols based on surface codes. First, we analyze the use of dislocations instead of holes to produce logical qubits, potentially reducing spacetime volume required. Dislocations induce defects which, in many respects, behave like Majorana quasi-particles. We construct circuits to implement these codes and present fault-tolerant measurement methods for these and other defects which may reduce spatial overhead. One advantage of these codes is that Hadamard gates take exactly $0$ time to implement. We numerically study the performance of these codes using a minimum weight and a greedy decoder using finite-size scaling. Second, we consider state injection of arbitrary ancillas to produce arbitrary rotations. This avoids the logarithmic (in precision) overhead in online cost required if $T$ gates are used to synthesize arbitrary rotations. While this has been considered before, we consider also the parallel performance of this protocol. Arbitrary ancilla injection leads to a probabilistic protocol in which there is a constant chance of success on each round; we use an amortized analysis to show that even in a parallel setting this leads to only a constant factor slowdown as opposed to the logarithmic slowdown that might be expected naively.

preprint2015arXiv

Towards Practical Quantum Variational Algorithms

The preparation of quantum states using short quantum circuits is one of the most promising near-term applications of small quantum computers, especially if the circuit is short enough and the fidelity of gates high enough that it can be executed without quantum error correction. Such quantum state preparation can be used in variational approaches, optimizing parameters in the circuit to minimize the energy of the constructed quantum state for a given problem Hamiltonian. For this purpose we propose a simple-to-implement class of quantum states motivated by adiabatic state preparation. We test its accuracy and determine the required circuit depth for a Hubbard model on ladders with up to 12 sites (24 spin-orbitals), and for small molecules. We find that this ansatz converges faster than previously proposed schemes based on unitary coupled clusters. While the required number of measurements is astronomically large for quantum chemistry applications to molecules, applying the variational approach to the Hubbard model (and related models) is found to be far less demanding and potentially practical on small quantum computers. We also discuss another application of quantum state preparation using short quantum circuits, to prepare trial ground states of models faster than using adiabatic state preparation.

preprint2014arXiv

Connecting Entanglement in Time and Space: Improving the Folding Algorithm

The "folding algorithm"\cite{fold1} is a matrix product state algorithm for simulating quantum systems that involves a spatial evolution of a matrix product state. Hence, the computational effort of this algorithm is controlled by the temporal entanglement. We show that this temporal entanglement is, in many cases, equal to the spatial entanglement of a modified Hamiltonian. This inspires a modification to the folding algorithm, that we call the "hybrid algorithm". We find that this leads to improved accuracy for the same numerical effort. We then use these algorithms to study relaxation in a transverse plus parallel field Ising model, finding persistent quasi-periodic oscillations for certain choices of initial conditions.

preprint2014arXiv

Improving Quantum Algorithms for Quantum Chemistry

We present several improvements to the standard Trotter-Suzuki based algorithms used in the simulation of quantum chemistry on a quantum computer. First, we modify how Jordan-Wigner transformations are implemented to reduce their cost from linear or logarithmic in the number of orbitals to a constant. Our modification does not require additional ancilla qubits. Then, we demonstrate how many operations can be parallelized, leading to a further linear decrease in the parallel depth of the circuit, at the cost of a small constant factor increase in number of qubits required. Thirdly, we modify the term order in the Trotter-Suzuki decomposition, significantly reducing the error at given Trotter-Suzuki timestep. A final improvement modifies the Hamiltonian to reduce errors introduced by the non-zero Trotter-Suzuki timestep. All of these techniques are validated using numerical simulation and detailed gate counts are given for realistic molecules.

preprint2014arXiv

Notes on Some Questions in Mathematical Physics and Quantum Information

This is a set of notes on some unrelated topics in mathematical physics, at varying levels of detail. First, I consider certain questions relating to the decay of correlation functions in matrix product states, in particular those generated by quantum expanders. This is discussed in relation to recent results of Brandao and Horodecki on area laws on systems with exponentially decaying correlation function\cite{areaexp}. Second, I consider some difficulties in trying to construct a tensor product state (or PEPS) describing a two-dimensional fermionic system with non-vanishing Hall conductance. Third, I present some relations between the theory of almost commuting matrices and that of vector bundles, making the connection between the classifications more explicit. Fourth, I present an open question about quantum channels, and some partial results.

preprint2014arXiv

The Trotter Step Size Required for Accurate Quantum Simulation of Quantum Chemistry

The simulation of molecules is a widely anticipated application of quantum computers. However, recent studies \cite{WBCH13a,HWBT14a} have cast a shadow on this hope by revealing that the complexity in gate count of such simulations increases with the number of spin orbitals $N$ as $N^8$, which becomes prohibitive even for molecules of modest size $N\sim 100$. This study was partly based on a scaling analysis of the Trotter step required for an ensemble of random artificial molecules. Here, we revisit this analysis and find instead that the scaling is closer to $N^6$ in worst case for real model molecules we have studied, indicating that the random ensemble fails to accurately capture the statistical properties of real-world molecules. Actual scaling may be significantly better than this due to averaging effects. We then present an alternative simulation scheme and show that it can sometimes outperform existing schemes, but that this possibility depends crucially on the details of the simulated molecule. We obtain further improvements using a version of the coalescing scheme of \cite{WBCH13a}; this scheme is based on using different Trotter steps for different terms. The method we use to bound the complexity of simulating a given molecule is efficient, in contrast to the approach of \cite{WBCH13a,HWBT14a} which relied on exponentially costly classical exact simulation.

preprint2013arXiv

Matrix Product Operators and Central Elements: Classical Description of a Quantum State

We study planar two-dimensional quantum systems on a lattice whose Hamiltonian is a sum of local commuting projectors of bounded range. We consider whether or not such a system has a zero energy ground state. To do this, we consider the problem as a one-dimensional problem, grouping all sites along a column into "supersites"; using $C^*$-algebraic methods (Bravyi and Vyalyi), we can solve this problem if we can characterize the central elements of the interaction algebra on these supersite. Unfortunately, these central elements may be very complex, making brute force impractical. Instead, we show a characterization of these elements in terms of matrix product operators with bounded bond dimension. This bound can be interpreted as a bound on the number of particle types in lattice theories with bounded Hilbert space dimension on each site. Topological order in this approach is related to the existence of certain central elements which cannot be "broken" into smaller pieces without creating an end excitation. Using this bound on bond dimension, we prove that several special cases of this problem are in NP, and we give part of a proof that the general case is in NP. Further, we characterize central elements that appear in certain specific models, including toric code and Levin-Wen models, as either product operators in the Abelian case or matrix product operators with low bond dimension in the non-Abelian case; this matrix product operator representation may have practical application in engineering the complicated multi-spin interactions in the Levin-Wen models.

preprint2013arXiv

Obstructions To Classically Simulating The Quantum Adiabatic Algorithm

We consider the adiabatic quantum algorithm for systems with "no sign problem", such as the transverse field Ising mode, and analyze the equilibration time for quantum Monte Carlo (QMC) on these systems. We ask: if the spectral gap is only inverse polynomially small, will equilibration methods based on slowly changing the Hamiltonian parameters in the QMC simulation succeed in a polynomial time? We show that this is not true, by constructing counter-examples. Some examples are Hamiltonians where the space of configurations where the wavefunction has non-negligible amplitude has a nontrivial fundamental group, causing the space of trajectories in imaginary time to break into disconnected components, with only negligible probability outside these components. For the simplest example we give with an abelian fundamental group, QMC does not equilibrate but still solves the optimization problem. More severe effects leading to failure to solve the optimization can occur when the fundamental group is a free group on two generators. Other examples where QMC fails have a trivial fundamental group, but still use ideas from topology relating group presentations to simplicial complexes. We define gadgets to realize these Hamiltonians as the effective low-energy dynamics of a transverse field Ising model. We present some analytic results on equilibration times which may be of some independent interest in the theory of equilibration of Markov chains. Conversely, we show that a small spectral gap implies slow equilibration at low temperature for some initial conditions and for a natural choice of local QMC updates.

preprint2013arXiv

Quantum Systems on Non-$k$-Hyperfinite Complexes: A Generalization of Classical Statistical Mechanics on Expander Graphs

We construct families of cell complexes that generalize expander graphs. These families are called non-$k$-hyperfinite, generalizing the idea of a non-hyperfinite (NH) family of graphs. Roughly speaking, such a complex has the property that one cannot remove a small fraction of points and be left with an object that looks $k-1$-dimensional at large scales. We then consider certain quantum systems on these complexes. A future goal is to construct a family of Hamiltonians such that every low energy state has topological order as part of an attempt to prove the quantum PCP conjecture. This goal is approached by constructing a toric code Hamiltonian with the property that every low energy state without vertex defects has topological order, a property that would not hold for any local system in any lattice $Z^d$ or indeed on any 1-hyperfinite complex. Further, such NH complexes find application in quantum coding theory. The hypergraph product codes[1] of Tillich and Zémor are generalized using NH complexes.

preprint2012arXiv

Solving Gapped Hamiltonians Locally

We show that any short-range Hamiltonian with a gap between the ground and excited states can be written as a sum of local operators, such that the ground state is an approximate eigenvector of each operator separately. We then show that the ground state of any such Hamiltonian is close to a generalized matrix product state. The range of the given operators needed to obtain a good approximation to the ground state is proportional to the square of the logarithm of the system size times a characteristic "factorization length". Applications to many-body quantum simulation are discussed. We also consider density matrices of systems at non-zero temperature.

preprint2012arXiv

Trivial Low Energy States for Commuting Hamiltonians, and the Quantum PCP Conjecture

We consider whether or not Hamiltonians which are sums of commuting projectors have "trivial" ground states which can be constructed by a local quantum circuit of bounded depth and range acting on a product state. While the toric code only has nontrivial ground states, commuting projector Hamiltonians which are sums of two-body interactions have trivial ground states. We define an "interaction complex" for a Hamiltonian, generalizing the interaction graph, and we show that if this complex can be continuously mapped to a 1-complex using a map with bounded diameter of pre-images then the Hamiltonian has a trivial ground state assuming one technical condition on the Hamiltonian (this condition holds for all stabilizer Hamiltonians, and we also prove the result for all Hamiltonians under an assumption on the 1-complex). While this includes cases considered by Ref., it also includes other Hamiltonians whose interaction complexes cannot be coarse-grained into the case of Ref. One motivation for this is the quantum PCP conjecture. Many commonly studied interaction complexes can be mapped to a 1-complex after removing a small fraction of sites. For commuting projector Hamiltonians on such complexes, a trivial ground state for the Hamiltonian with those sites removed is a low energy trivial state for the original Hamiltonian. Such states can act as a classical witness to the existence of a low energy state. While this result applies only to commuting Hamiltonians, it suggests that to prove a quantum PCP conjecture one should consider interaction complexes which cannot be mapped to 1-complexes after removing a small fraction of cells. We define this more precisely below; in a sense this generalizes the idea of an expander graph. Surprisingly, such complexes do exist as will be shown elsewhere, and have useful properties in quantum coding theory.

preprint2011arXiv

Making Almost Commuting Matrices Commute

Suppose two Hermitian matrices $A,B$ almost commute ($\Vert [A,B] \Vert \leq δ$). Are they close to a commuting pair of Hermitian matrices, $A',B'$, with $\Vert A-A' \Vert,\Vert B-B'\Vert \leq ε$? A theorem of H. Lin shows that this is uniformly true, in that for every $ε>0$ there exists a $δ>0$, independent of the size $N$ of the matrices, for which almost commuting implies being close to a commuting pair. However, this theorem does not specify how $δ$ depends on $ε$. We give uniform bounds relating $δ$ and $ε$. We provide tighter bounds in the case of block tridiagonal and tridiagonal matrices and a fully constructive method in that case. Within the context of quantum measurement, this implies an algorithm to construct a basis in which we can make a {\it projective} measurement that approximately measures two approximately commuting operators simultaneously. Finally, we comment briefly on the case of approximately measuring three or more approximately commuting operators using POVMs (positive operator-valued measures) instead of projective measurements.

preprint2011arXiv

Topological Order at Non-zero Temperature

We propose a definition for topological order at nonzero temperature in analogy to the usual zero temperature definition that a state is topologically ordered, or "nontrivial", if it cannot be transformed into a product state (or a state close to a product state) using a local (or approximately local) quantum circuit. We prove that any two dimensional Hamiltonian which is a sum of commuting local terms is not topologically ordered at $T>0$. We show that such trivial states cannot be used to store quantum information using certain stringlike operators. This definition is not too restrictive, however, as the four dimensional toric code does have a nontrivial phase at nonzero temperature.

preprint2010arXiv

A short proof of stability of topological order under local perturbations

Recently, the stability of certain topological phases of matter under weak perturbations was proven. Here, we present a short, alternate proof of the same result. We consider models of topological quantum order for which the unperturbed Hamiltonian $H_0$ can be written as a sum of local pairwise commuting projectors on a $D$-dimensional lattice. We consider a perturbed Hamiltonian $H=H_0+V$ involving a generic perturbation $V$ that can be written as a sum of short-range bounded-norm interactions. We prove that if the strength of $V$ is below a constant threshold value then $H$ has well-defined spectral bands originating from the low-lying eigenvalues of $H_0$. These bands are separated from the rest of the spectrum and from each other by a constant gap. The width of the band originating from the smallest eigenvalue of $H_0$ decays faster than any power of the lattice size.

preprint2010arXiv

Disordered Topological Insulators via $C^*$-Algebras

The theory of almost commuting matrices can be used to quantify topological obstructions to the existence of localized Wannier functions with time-reversal symmetry in systems with time-reversal symmetry and strong spin-orbit coupling. We present a numerical procedure that calculates a Z_2 invariant using these techniques, and apply it to a model of HgTe. This numerical procedure allows us to access sizes significantly larger than procedures based on studying twisted boundary conditions. Our numerical results indicate the existence of a metallic phase in the presence of scattering between up and down spin components, while there is a sharp transition when the system decouples into two copies of the quantum Hall effect. In addition to the Z_2 invariant calculation in the case when up and down components are coupled, we also present a simple method of evaluating the integer invariant in the quantum Hall case where they are decoupled.

preprint2010arXiv

Entanglement vs. gap for one-dimensional spin systems

We study the relationship between entanglement and spectral gap for local Hamiltonians in one dimension. The area law for a one-dimensional system states that for the ground state, the entanglement of any interval is upper-bounded by a constant independent of the size of the interval. However, the possible dependence of the upper bound on the spectral gap Delta is not known, as the best known general upper bound is asymptotically much larger than the largest possible entropy of any model system previously constructed for small Delta. To help resolve this asymptotic behavior, we construct a family of one-dimensional local systems for which some intervals have entanglement entropy which is polynomial in 1/Delta, whereas previously studied systems, such as free fermion systems or systems described by conformal field theory, had the entropy of all intervals bounded by a constant times log(1/Delta).

preprint2010arXiv

Locality in Quantum Systems

These lecture notes focus on the application of ideas of locality, in particular Lieb-Robinson bounds, to quantum many-body systems. We consider applications including correlation decay, topological order, a higher dimensional Lieb-Schultz-Mattis theorem, and a nonrelativistic Goldstone theorem. The emphasis is on trying to show the ideas behind the calculations. As a result, the proofs are only sketched with an emphasis on the intuitive ideas behind them, and in some cases we use techniques that give very slightly weaker bounds for simplicity. This is a preliminary version of the lecture notes, with the goal of getting the notes out close to the end of the school. Comments welcome.

preprint2010arXiv

Quasi-adiabatic Continuation for Disordered Systems: Applications to Correlations, Lieb-Schultz-Mattis, and Hall Conductance

We present a possible definition of a mobility gap for a many-body quantum system, in analogy to definitions of dynamical localization for single particle systems. Using this definition, we construct "corrected" quasi-adiabatic continuation operators. Under an appropriate definition of a unique ground state, we show how to introduce virtual fluxes. Armed with these results, we can directly carry over previous results in the case of a spectral gap. We present a proof of decay of correlation functions and we present a proof of Hall conductance quantization under very mild density-of-states assumptions defined later. We also generalize these definitions to the case of a "bulk mobility gap", in the case of a system with boundaries, and present a proof of Hall conductance quantization on an annulus under appropriate assumptions. Further, we present a new "optimized" quasi-adiabatic continuation operator which simplifies previous estimates and tightens bounds in certain cases. This is presented in an appendix which can be read independently of the rest of the paper as it also improves estimates in the case of systems with a spectral gap. This filter function used decays in time at least as fast as ${\cal O}(\exp(-t^α))$ for all $α<1$, a class of decay called subexponential (a tighter description of what is possible is below). Using this function it is possible to tighten recent estimates of the Hall conductance quantization for gapped systems\cite{hall} to a decay which is subexponential in system size.

preprint2010arXiv

Topological Insulators and C^*-Algebras: Theory and Numerical Practice

We apply ideas from $C^*$-algebra to the study of disordered topological insulators. We extract certain almost commuting matrices from the free Fermi Hamiltonian, describing band projected coordinate matrices. By considering topological obstructions to approximating these matrices by exactly commuting matrices, we are able to compute invariants quantifying different topological phases. We generalize previous two dimensional results to higher dimensions; we give a general expression for the topological invariants for arbitrary dimension and several symmetry classes, including chiral symmetry classes, and we present a detailed $K$-theory treatment of this expression for time reversal invariant three dimensional systems. We can use these results to show non-existence of localized Wannier functions for these systems. We use this approach to calculate the index for time-reversal invariant systems with spin-orbit scattering in three dimensions, on sizes up to $12^3$, averaging over a large number of samples. The results show an interesting separation between the localization transition and the point at which the average index (which can be viewed as an "order parameter" for the topological insulator) begins to fluctuate from sample too sample, implying the existence of an unsuspected quantum phase transition separating two different delocalized phases in this system. One of the particular advantages of the $C^*$-algebraic technique that we present is that it is significantly faster in practice than other methods of computing the index, allowing the study of larger systems. In this paper, we present a detailed discussion of numerical implementation of our method.

preprint2009arXiv

Almost Commuting Matrices, Localized Wannier Functions, and the Quantum Hall Effect

For models of non-interacting fermions moving within sites arranged on a surface in three dimensional space, there can be obstructions to finding localized Wannier functions. We show that such obstructions are $K$-theoretic obstructions to approximating almost commuting, complex-valued matrices by commuting matrices, and we demonstrate numerically the presence of this obstruction for a lattice model of the quantum Hall effect in a spherical geometry. The numerical calculation of the obstruction is straightforward, and does not require translational invariance or introducing a flux torus. We further show that there is a $Z_2$ index obstruction to approximating almost commuting self-dual matrices by exactly commuting self-dual matrices, and present additional conjectures regarding the approximation of almost commuting real and self-dual matrices by exactly commuting real and self-dual matrices. The motivation for considering this problem is the case of physical systems with additional antiunitary symmetries such as time reversal or particle-hole conjugation. Finally, in the case of the sphere--mathematically speaking three almost commuting Hermitians whose sum of square is near the identity--we give the first quantitative result showing this index is the only obstruction to finding commuting approximations. We review the known non-quantitative results for the torus.

preprint2009arXiv

Light Cone Matrix Product

We show how to combine the light-cone and matrix product algorithms to simulate quantum systems far from equilibrium for long times. For the case of the XXZ spin chain at $Δ=0.5$, we simulate to a time of $\approx 22.5$. While part of the long simulation time is due to the use of the light-cone method, we also describe a modification of the iTEBD algorithm with improved numerical stability, and we describe how to incorporate symmetry into this algorithm. While statistical sampling error means that we are not yet able to make a definite statement, the behavior of the simulation at long times indicates the appearance of either "revivals" in the order parameter as predicted previously or of a distinct shoulder in the decay of the order parameter.

preprint2009arXiv

Quantum Adiabatic Computation With a Constant Gap is Not Useful in One Dimension

We show that it is possible to use a classical computer to efficiently simulate the adiabatic evolution of a quantum system in one dimension with a constant spectral gap, starting the adiabatic evolution from a known initial product state. The proof relies on a recently proven area law for such systems, implying the existence of a good matrix product representation of the ground state, combined with an appropriate algorithm to update the matrix product state as the Hamiltonian is changed. This implies that adiabatic evolution with such Hamiltonians is not useful for universal quantum computation. Therefore, adiabatic algorithms which are useful for universal quantum computation either require a spectral gap tending to zero or need to be implemented in more than one dimension (we leave open the question of the computational power of adiabatic simulation with a constant gap in more than one dimension).

preprint2001arXiv

Entropic Tightening of Vibrated Chains

We investigate experimentally the distribution of configurations of a ring with an elementary topological constraint, a ``figure-8'' twist. Using vibrated granular chains, which permit controlled preparation and direct observation of such a constraint, we show that configurations where one of the loops is tight and the second is large are strongly preferred. This agrees with recent predictions for equilibrium properties of topologically-constrained polymers. However, the dynamics of the tightening process weakly violate detailed balance, a signature of the nonequilibrium nature of this system.

preprint2000arXiv

Dirac, Anderson, and Goldstone on the Kagome

We show that there exists a long-range RVB state for the kagome lattice spin-1/2 Heisenberg antiferromagnet for which the spinons have a massless Dirac spectrum. By considering various perturbations of the RVB state which give mass to the fermions by breaking a symmetry, we are able to describe a wide-ranging class of known states on the kagome lattice, including spin-Peierls solid and chiral spin liquid states. Using an RG treatment of fluctuations about the RVB state, we propose yet a different symmetry breaking pattern and show how collective excitations about this state account for the gapless singlet modes seen experimentally and numerically. We make further comparison with numerics for Chern numbers, dimer-dimer correlation functions, the triplet gap, and other quantities. To accomplish these calculations, we propose a variant of the SU(N) theory which enables us to include many of the effects of Gutzwiller projection at the mean-field level.

preprint1999arXiv

Ground State and Spin Glass Phase of the Large N Infinite Range Spin Glass Via Supersymmetry

The large N infinite range spin glass is considered, in particular the number of spin components k needed to form the ground state and the sample-to-sample fluctuations in the Lagrange multiplier field on each site. The physical significance of k for the correlation functions is discussed. The difference between the large N and spherical spin glass is emphasized; a slight difference between the average Lagrange multiplier of the large N and spherical spin glasses is derived, leading to a slight increase in the energy of the ground state compared to the naive expectation. Further, there is a change in the low energy density of excitations in the large N system. A form of level repulsion, similar to that found in random matrix theory, is found to exist in this system, surviving interactions. Even though the system is an interacting one, a supersymmetric formalism is developed to deal with the problem of averaging over disorder.

preprint1999arXiv

Renormalization Group for Large N Strongly Commensurate Dirty Boson Model

The large N sigma model, in D<4 space-time dimensions, with disorder a function of d space dimensions, is analyzed via a renormalization group treatment. Critical exponents for average quantities are calculated, first to lowest order and then to all orders, in $ε=D-2 - d/2$. In particular, it is found that $νd =2$. When D=d+1, this model is equivalent to a large N limit of the strongly commensurate dirty boson problem.

preprint1997arXiv

Non-Hermitian Fermion Mapping for One-Component Plasma

The two-dimensional one-component logarithmic Coulomb gas is mapped onto a non-hermitian fermionic field theory. At $β=2$, the field theory is free. Correlation functions are calculated and a perturbation theory is discussed for extending to other $β$. A phase transition is found at the mean-field level at large $β$. Some results are extended to spaces of constant negative curvature.

preprint1996arXiv

Bragg resonances for tunneling between edges of a 2D Quantum Hall system

A theory is presented for tunneling between compressible regions on the sides of a narrow incompressible Quantum Hall strip. Assuming that electron interactions lead to formation of a Wigner crystal on the edges of the compressible regions, we consider the situation when the non-conservation of electron momentum required for transport is provided, in the absence of disorder, by umklapp scattering on the crystal. The momentum given to the crystal is quantized due to the Bragg condition, which leads to resonances in tunneling conductivity as a function of the incompressible strip width, similar to those reported recently by N. Zhitenev, M. Brodsky, R. Ashoori, and M. Melloch.