Steven H. Simon

Global Phase Diagram of the Normal State of Twisted Bilayer Graphene

We investigate the full doping and strain-dependent phase diagram of the normal state of magic-angle twisted bilayer graphene (TBG). Using comprehensive Hartree-Fock calculations, we show that at temperatures where superconductivity is absent the global phase structure can be understood based on the competition and coexistence between three types of intertwined orders: a fully symmetric phase, spatially uniform flavor-symmetry-breaking states, and an incommensurate Kekulé spiral (IKS) order. For small strain, the IKS phase, recently proposed as a candidate order at all non-zero integer fillings of the moiré unit cell, is found to be ubiquitous for non-integer doping as well. We demonstrate that the corresponding electronic compressibility and Fermi surface structure are consistent with the 'cascade' physics and Landau fans observed experimentally. This suggests a unified picture of the phase diagram of TBG in terms of IKS order.

Quantifying information scrambling via Classical Shadow Tomography on Programmable Quantum Simulators

We develop techniques to probe the dynamics of quantum information, and implement them experimentally on an IBM superconducting quantum processor. Our protocols adapt shadow tomography for the study of time evolution channels rather than of quantum states, and rely only on single-qubit operations and measurements. We identify two unambiguous signatures of quantum information scrambling, neither of which can be mimicked by dissipative processes, and relate these to many-body teleportation. By realizing quantum chaotic dynamics in experiment, we measure both signatures, and support our results with numerical simulations of the quantum system. We additionally investigate operator growth under this dynamics, and observe behaviour characteristic of quantum chaos. As our methods require only a single quantum state at a time, they can be readily applied on a wide variety of quantum simulators.

Statistical mechanics of dimers on quasiperiodic Ammann-Beenker tilings

We study classical dimers on two-dimensional quasiperiodic Ammann-Beenker (AB) tilings. Despite the lack of periodicity we prove that each infinite tiling admits 'perfect matchings' in which every vertex is touched by one dimer. We introduce an auxiliary 'AB $^*$ ' tiling obtained from the AB tiling by deleting all 8-fold coordinated vertices. The AB $^*$ tiling is again two-dimensional, infinite, and quasiperiodic. The AB $^*$ tiling has a single connected component, which admits perfect matchings. We find that in all perfect matchings, dimers on the AB $^*$ tiling lie along disjoint one-dimensional loops and ladders, separated by 'membranes', sets of edges where dimers are absent. As a result, the dimer partition function of the AB $^*$ tiling factorizes into the product of dimer partition functions along these structures. We compute the partition function and free energy per edge on the AB $^*$ tiling using an analytic transfer matrix approach. Returning to the AB tiling, we find that membranes in the AB $^*$ tiling become 'pseudomembranes', sets of edges which collectively host at most one dimer. This leads to a remarkable discrete scale-invariance in the matching problem. The structure suggests that the AB tiling should exhibit highly inhomogenous and slowly decaying connected dimer correlations. Using Monte Carlo simulations, we find evidence supporting this supposition in the form of connected dimer correlations consistent with power law behaviour. Within the set of perfect matchings we find quasiperiodic analogues to the staggered and columnar phases observed in periodic systems.

Straightening Out the Frobenius-Schur Indicator

The Frobenius-Schur indicator is a parameter $κ_a=\pm 1$ assigned to each self-dual particle $a$ in a TQFT. If $κ_a$ is negative then straightening out a timelike zig-zag in the worldline of a particle of type $a$ can incur a minus sign and in this case the amplitude associated with the diagram is not invariant under deformation. This has caused some confusion about the topological invariance of even simple theories to space-time deformations. We clarify that, given a TQFT with negative Frobenius-Schur indicators, there are two distinct conventions commonly used to interpret a spacetime diagram as a physical amplitude, only one of which is isotopy invariant. We clarify in what sense TQFTs based on Chern-Simons theory with negative Frobenius-Schur indicators are isotopy invariant, and we explain how the Frobenius-Schur indicator is intimately linked with the need to frame world-lines in Chern-Simons theory. Further, in the non-isotopy-invariant interpretation of the diagram algebra we show how a trick of bookkeeping can usually be invoked to push minus signs onto the diagrammatic value of a loop (the "loop weight"), such that most of the evaluation of a diagram does not incur minus signs from straightening zig-zags, and only at the last step minus signs are added. We explain the conditions required for this to be possible. We then further examine what is required in order for a theory to have full isotopy invariance of planar spacetime diagrams, and discover that, if we have successfully pushed the signs from zig-zags onto the loop weight, the only possible obstruction to this is given by an object related to vertices, known as the "third Frobenius-Schur indicator". We finally discuss the extent to which this gives us full isotopy invariance for braided theories.

The Superconductivity of Sr $_2$ RuO $_4$ Under $c$ -Axis Uniaxial Stress

Applying in-plane uniaxial pressure to strongly correlated low-dimensional systems has been shown to tune the electronic structure dramatically. For example, the unconventional superconductor Sr $_2$ RuO $_4$ can be tuned through a single Van Hove singularity which results in a strong enhancement of both $T_\text{c}$ and $H_\text{c2}$ . Out-of-plane ( $c$ axis) uniaxial pressure is expected to tune the quasi-two-dimensional structure even more strongly, causing it to approach two Van Hove singularities simultaneously. Here we achieve a record value of $3.2$ GPa compression along the $c$ axis of Sr $_2$ RuO $_4$ . Although the rise in $H_\text{c2}$ shows that we are indeed approaching the van Hove points, $T_\text{c}$ is suppressed, a result that contradicts expectations based on simple two-dimensional models. As a first attempt to take the third dimension into account, we present three-dimensional calculations in the weak interaction limit, and discuss the extent to which they are consistent with observation. Our experimental results highlight the importance of out-of-plane effects in low-dimensional systems in general and provide new constraints on theories of the pairing interaction in Sr $_2$ RuO $_4$ .

preprint2021arXiv

Kekulé spiral order at all nonzero integer fillings in twisted bilayer graphene

We study magic angle graphene in the presence of both strain and particle-hole symmetry breaking due to non-local inter-layer tunneling. We perform a self-consistent Hartree-Fock study that incorporates these effects alongside realistic interaction and substrate potentials, and explore a comprehensive set of competing orders including those that break translational symmetry at arbitrary wavevectors. We find that at all non-zero integer fillings very small strains, comparable to those measured in scanning tunneling experiments, stabilize a fundamentally new type of time-reversal symmetric and spatially non-uniform order. This order, which we dub the 'incommensurate Kekulé spiral' (IKS) order, spontaneously breaks both the emergent valley-charge conservation and moiré translation symmetries, but preserves a modified translation symmetry $\hat{T}'$ -- which simultaneously shifts the spatial coordinates and rotates the $U(1)$ angle which characterizes the spontaneous inter-valley coherence. We discuss the phenomenological and microscopic properties of this order. We argue that our findings are consistent with all experimental observations reported so far, suggesting a unified explanation of the global phase diagram in terms of the IKS order.

Approximating observables on eigenstates of large many-body localized systems

Eigenstates of fully many-body localized (FMBL) systems can be organized into spin algebras based on quasilocal operators called l-bits. These spin algebras define quasilocal l-bit measurement ( $τ^z_i$ ) and l-bit flip ( $τ^x_i$ ) operators. For a disordered Heisenberg spin chain in the MBL regime we approximate l-bit flip operators by finding them exactly on small windows of systems and extending them onto the whole system by exploiting their quasilocal nature. We subsequently use these operators to represent approximate eigenstates. We then describe a method to calculate products of local observables on these eigenstates for systems of size $L$ in $O(L^2)$ time. This algorithm is used to compute the error of the approximate eigenstates.

Behavior of l-bits near the many-body localization transition

Eigenstates of fully many-body localized (FMBL) systems are described by quasilocal operators $τ_i^z$ (l-bits), which are conserved exactly under Hamiltonian time evolution. The algebra of the operators $τ_i^z$ and $τ_i^x$ associated with l-bits ( $\boldsymbolτ_i$ ) completely defines the eigenstates and the matrix elements of local operators between eigenstates at all energies. We develop a non-perturbative construction of the full set of l-bit algebras in the many-body localized phase for the canonical model of MBL. Our algorithm to construct the Pauli-algebra of l-bits combines exact diagonalization and a tensor network algorithm developed for efficient diagonalization of large FMBL Hamiltonians. The distribution of localization lengths of the l-bits is evaluated in the MBL phase and used to characterize the MBL-to-thermal transition.

Signatures of the Many-body Localized Regime in Two Dimensions

Lessons from Anderson localization highlight the importance of dimensionality of real space for localization due to disorder. More recently, studies of many-body localization have focussed on the phenomenon in one dimension using techniques of exact diagonalization and tensor networks. On the other hand, experiments in two dimensions have provided concrete results going beyond the previously numerically accessible limits while posing several challenging questions. We present the first large-scale numerical examination of a disordered Bose-Hubbard model in two dimensions realized in cold atoms, which shows entanglement based signatures of many-body localization. By generalizing a low-depth quantum circuit to two dimensions we approximate eigenstates in the experimental parameter regimes for large systems, which is beyond the scope of exact diagonalization. A careful analysis of the eigenstate entanglement structure provides an indication of the putative phase transition marked by a peak in the fluctuations of entanglement entropy in a parameter range consistent with experiments.

Transport properties of multilayer graphene

We apply the semi-classical quantum Boltzmann formalism for the computation of transport properties to multilayer graphene. We compute the electrical conductivity as well as the thermal conductivity and thermopower for Bernal-stacked multilayers with an even number of layers. We show that the window for hydrodynamic transport in multilayer graphene is similar to the case of bilayer graphene. We introduce a simple hydrodynamic model which we dub the multi-fluid model and which can be used to reproduce the results for the electrical conductivity and thermopower from the quantum Boltzmann equation.

Classical Dimers on Penrose Tilings

We study the classical dimer model on rhombic Penrose tilings, whose edges and vertices may be identified with those of a bipartite graph. We find that Penrose tilings do not admit perfect matchings (defect-free dimer coverings). Instead, their maximum matchings have a monomer density of $81-50φ\approx 0.098$ in the thermodynamic limit, with $φ=\left(1+\sqrt{5}\right)/2$ the golden ratio. Maximum matchings divide the tiling into a fractal of nested closed regions bounded by loops that cannot be crossed by monomers. These loops connect second-nearest neighbour even-valence vertices, each of which lies on such a loop. Assigning a charge to each monomer with a sign fixed by its bipartite sublattice, we find that each bounded region has an excess of one charge, and a corresponding set of monomers, with adjacent regions having opposite net charge. The infinite tiling is charge neutral. We devise a simple algorithm for generating maximum matchings, and demonstrate that maximum matchings form a connected manifold under local monomer-dimer rearrangements. We show that dart-kite Penrose tilings feature an imbalance of charge between bipartite sub-lattices, leading to a minimum monomer density of $\left(7-4φ\right)/5\approx 0.106$ all of one charge.

Energetics of Pfaffian-AntiPfaffian Domains

In several recent works it has been proposed that, due to disorder, the experimentally observed nu=5/2 quantum Hall state could be microscopically composed of domains of Pfaffian order along with domains of AntiPfaffian order. We numerically examine the energetics required for forming such domains and conclude that for the parameters appropriate for recent experiments, such domains would not occur.

Partial Equilibration of the Anti-Pfaffian edge due to Majorana Disorder

We consider electrical and thermal equilibration of the edge modes of the Anti-Pfaffian quantum Hall state at $ν=5/2$ due to tunneling of the Majorana edge mode to trapped Majorana zero modes in the bulk. Such tunneling breaks translational invariance and allows scattering between Majorana and other edge modes in such a way that there is a parametric difference between the length scales for equilibration of charge and heat transport between integer and Bose mode on the one hand, and for thermal equilibration of the Majorana edge mode on the other hand. We discuss a parameter regime in which this mechanism could explain the recent observation of quantized heat transport [Banerjee et all, Nature 559, 7713 (2018)].

Quantum Boltzmann equation for bilayer graphene

A-B stacked bilayer graphene has massive electron and hole-like excitations with zero gap in the nearest-neighbor hopping approximation. In equilibrium, the quasiparticle occupation approximately follows the usual Fermi-Dirac distribution. In this paper we consider perturbing this equilibrium distribution so as to determine DC transport coefficients near charge neutrality. We consider the regime $β|μ| \lesssim 1$ (with $β$ the inverse temperature and $μ$ the chemical potential) where there is not a well formed Fermi surface. Starting from the Kadanoff-Baym equations, we obtain the quantum Boltzmann equation of the electron and hole distribution functions when the system is weakly perturbed out of equilibrium. The effect of phonons, disorder, and boundary scattering for finite sized systems are incorporated through a generalized collision integral. The transport coefficients, including the electrical and thermal conductivity, thermopower, and shear viscosity, are calculated in the linear response regime. We also extend the formalism to include an external magnetic field. We present results from numerical solutions of the quantum Boltzmann equation. Finally, we derive a simplified two-fluid hydrodynamic model appropriate for this system, which reproduces the salient results of the full numerical calculations.

Transport in bilayer graphene near charge neutrality: Which scattering mechanisms are important?

Using the semiclassical quantum Boltzmann equation (QBE), we numerically calculate the DC transport properties of bilayer graphene near charge neutrality. We find, in contrast to prior discussions, that phonon scattering is crucial even at temperatures below 40K. Nonetheless, electron-electron scattering still dominates over phonon collisions allowing a hydrodynamic approach. We introduce a simple two-fluid hydrodynamic model of electrons and holes interacting via Coulomb drag and compare our results to the full QBE calculation. We show that the two-fluid model produces quantitatively accurate results for conductivity, thermopower, and thermal conductivity.

preprint2016arXiv

Driven impurity in an ultracold 1D Bose gas with intermediate interaction strength

We study a single impurity driven by a constant force through a 1D Bose gas using a Lieb-Liniger based approach. Our calculaton is exact in the interaction amongst the particles in the Bose gas, and is perturbative in the interaction between the gas and the impurity. In contrast to previous studies of this problem, we are able to handle arbitrary interaction strength for the Bose gas. We find very good agreement with recent experiments [Phys. Rev. Lett. 103, 150601 (2009)].

Composite fermion model for entanglement spectrum of fractional quantum Hall states

We show that the entanglement spectrum associated with a certain class of strongly correlated many-body states --- the wave functions proposed by Laughlin and Jain to describe the fractional quantum Hall effect --- can be very well described in terms of a simple model of non-interacting (or weakly interacting) composite fermions.

Fractional Chern Insulators in Bands with Zero Berry Curvature

Even if a noninteracting system has zero Berry curvature everywhere in the Brillouin zone, it is possible to introduce interactions that stabilise a fractional Chern insulator. These interactions necessarily break time-reversal symmetry (either spontaneously or explicitly) and have the effect of altering the underlying band structure. We outline a number of ways in which this may be achieved, and show how similar interactions may also be used to create a (time-reversal symmetric) fractional topological insulator. While our approach is rigorous in the limit of long range interactions, we show numerically that even for short range interactions a fractional Chern insulator can be stabilised in a band with zero Berry curvature.

Fractional quantum Hall effect in bilayer graphene beyond the single Landau level approximation

Bilayer graphene has been predicted to give unprecedented tunability of the electron-electron interaction with the help of external parameters, allowing one to stabilize different fractional quantum Hall states. Recent experimental works make theoretical analysis of such systems extremely relevant. In this paper we describe a methodology for investigating the possibility of realizing specific fractional quantum Hall states in bilayer graphene taking into account polarization effects and virtual interband transitions. We apply this methodology to explore the possibility of realizing the Moore-Read Pfaffian state in bilayer graphene.

Hidden order and flux attachment in symmetry protected topological phases: a Laughlin-like approach

Topological phases of matter are distinct from conventional ones by their lack of a local order parameter. Still in the quantum Hall effect, hidden order parameters exist and constitute the basis for the celebrated composite-particle approach. Whether similar hidden orders exist in 2D and 3D symmetry protected topological phases (SPTs) is a largely open question. Here we introduce a new approach for generating SPT groundstates, based on a generalization of the Laughlin wavefunction. This approach gives a simple and unifying picture of some classes of SPTs in 1D and 2D, and reveals their hidden order and flux attachment structures. For the 1D case, we derive exact relations between the wavefunctions obtained in this manner and group cohomology wavefunctions, as well as matrix product state classification. For the 2D Ising SPT, strong analytical and numerical evidence is given to show that the wavefunction obtained indeed describes the desired SPT. The Ising SPT then appears as a state with quasi-long-range order in composite degrees of freedom consisting of Ising-symmetry charges attached to Ising-symmetry fluxes.

Large Chern Number and Edge Currents in Sr2RuO4

We show from a weak coupling microscopic calculation that the most favored chiral superconducting order parameter in Sr2RuO4 has Chern number |C|=7. The two dominant components of this order parameter are given by sin(3 kx) + i sin(3 ky) and sin(kx) cos(ky) + i sin(ky) cos(kx) and lie in the same irreducible representation Eu of the tetragonal point group as the usually assumed gap function sin(kx) + i sin(ky). While the latter gap function leads to C=1, the two former lead to C=-7, which is also allowed for an Eu gap function since the tetragonal symmetry only fixes C modulo 4. Since it was shown that the edge currents of a |C|>1 superconductor vanish exactly in the continuum limit, and can be strongly reduced on the lattice, this form of order parameter could help resolve the conflict between experimental observation of time-reversal symmetry breaking and yet the absence of observed edge currents in Sr2RuO4.

Comment on "Elementary formula for the Hall conductivity of interacting systems"

In a recent paper by Neupert, Santos, Chamon, and Mudry [Phys. Rev. B 86, 165133 (2012)] it is claimed that there is an elementary formula for the Hall conductivity of fractional Chern insulators. We show that the proposed formula cannot generally be correct, and we suggest one possible source of the error. Our reasoning can be generalized to show no quantity (such as Hall conductivity) expected to be constant throughout an entire phase of matter can possibly be given as the expectation of any time independent short ranged operator.

Exact Solutions of Fractional Chern Insulators: Interacting Particles in the Hofstadter Model at Finite Size

We show that all the bands of the Hofstadter model on the torus have an exactly flat dispersion and Berry curvature when a special system size is chosen. This result holds for any hopping and Chern number. Our analysis therefore provides a simple rule for choosing a particularly advantageous system size when designing a Hofstadter system whose size is controllable, like a qubit lattice or an optical cavity array. The density operators projected onto the flat bands obey exactly the Girvin-MacDonald-Platzman algebra, like for Landau levels in the continuum in the case of $C=1$ , or obey its straightforward generalization in the case of $C>1$ . This allows a mapping between density-density interaction Hamiltonians for particles in the Hofstatder model and in a continuum Landau level. By using the well-known pseudopotential construction in the latter case, we obtain fractional Chern insulator phases, the lattice counterpart of fractional quantum Hall phases, that are exact zero-energy ground states of the Hofstadter model with certain interactions. Finally, the addition of a harmonic trapping potential is shown to lead to an appealingly symmetric description in which a new Hofstadter model appears in momentum space.

Focus on topological quantum computation

Topological quantum computation started as a niche area of research aimed at employing particles with exotic statistics, called anyons, for performing quantum computation. Soon it evolved to include a wide variety of disciplines. Advances in the understanding of anyon properties inspired new quantum algorithms and helped in the characterisation of topological phases of matter and their experimental realisation. The conceptual appeal of topological systems as well as their promise for building fault-tolerant quantum technologies fuelled the fascination in this field. This `focus on' brings together several of the latest developments in the field and facilitates the synergy between different approaches.

Induced self-stabilization in fractional quantum Hall states of light

Recent progress in nanoscale quantum optics and superconducting qubits has made the creation of strongly correlated, and even topologically ordered, states of photons a real possibility. Many of these states are gapped and exhibit anyon excitations, which could be used for a robust form of quantum information processing. However, while numerous qubit array proposals exist to engineer the Hamiltonian for these systems, the question of how to stabilize the many-body ground state of these photonic quantum simulators against photon losses remains largely unanswered. We here propose a simple mechanism which achieves this goal for abelian and non-abelian fractional quantum Hall states of light. Our construction uses a uniform two-photon drive field to couple the qubits of the primary lattice with an auxiliary ``shadow" lattice, composed of qubits with a much faster loss rate than the qubits of the primary quantum simulator itself. This coupling causes hole states created by photon losses to be rapidly refilled, and the system's many-body gap prevents further photons from being added once the strongly correlated ground state is reached. The fractional quantum Hall state (with a small, transient population of quasihole excitations) is thus the most stable state of the system, and all other configurations will relax toward it over time. The physics described here could be implemented in a circuit QED architecture, and the device parameters needed for our scheme to succeed are in reach of current technology. We also propose a simple 6 qubit device, which could easily be built in the near future, that can act as a proof of principle for our scheme.

Josephson Coupled Moore-Read States

We study a quantum Hall bilayer system of bosons at total filling factor $ν= 1$ , and study the phase that results from short ranged pair-tunneling combined with short ranged interlayer interactions. We introduce two exactly solvable model Hamiltonians which both yield the coupled Moore-Read state [Phys.~Rev.~Lett.~{\bf 108}, 256809 (2012)] as a ground state, when projected onto fixed particle numbers in each layer. One of these Hamiltonians describes a gapped topological phase while the other is gapless. However, on introduction of a pair tunneling term, the second system becomes gapped and develops the same topological order as the gapped Hamiltonian. Supported by the exact solution of the full zero-energy quasihole spectrum and a conformal field theory approach, we develop an intuitive picture of this system as two coupled composite fermion superconductors. In this language, pair tunneling provides a Josephson coupling of the superconducting phases of the two layers, and gaps out the Goldstone mode associated with particle transport between the layers. In particular, this implies that quasiparticles are confined between the layers. In the bulk, the resulting phase has the topological order of the Halperin 220 phase with $U(1)_2\times U(1)_2$ topological order, but it is realized in the symmetric/antisymmetric-basis of the layer index. Consequently, the edge spectrum at a fixed particle number reveals an unexpected $U(1)_4 \times U(1)$ structure.

Non-equilibrium noise in the (non-)Abelian fractional quantum Hall effect

We analyse the noise of the edge current of a generic fractional quantum Hall state in a tunnelling point contact system. We show that the non-symmetrized noise in the edge current for the system out-of-equilibrium is completely determined by the noise in the tunnelling current and the Nyquist-Johnson (equilibrium) noise of the edge current. Simply put, the noise in the tunnelling current does not simply add up the equilibrium noise of the edge current. A correction term arises associated with the correlation between the tunnelling current and the edge current. We show, using a non-equilibrium Ward identity, that this correction term is determined by the anti-symmetric part of the noise in the tunnelling current. This leads to a non-equilibrium fluctuation-dissipation theorem and related expressions for the excess and shot noise of the noise in the edge current. Our approach makes use of simple properties of the edge, such as charge conservation and chirality, and applies to generic constructions of the edge theory which includes edges of non-Abelian states and edges with multiple charged channels. Two important tools we make use of are the non-equilibrium Kubo formula and the non-equilibrium Ward identity. We discuss these identities in the appendix.

Passive correction of quantum logical errors in a driven, dissipative system: a blueprint for an analog quantum code fabric

A physical realization of self correcting quantum code would be profoundly useful for constructing a quantum computer. In this theoretical work, we provide a partial solution to major challenges preventing self correcting quantum code from being engineered in realistic devices. We consider a variant of Kitaev's toric code coupled to propagating bosons, which induce a long-ranged interaction between anyonic defects. By coupling the primary quantum system to an engineered dissipation source through resonant energy transfer, we demonstrate a "rate barrier" which leads to a potentially enormous increase in the system's quantum state lifetime through purely passive quantum error correction, even when coupled to an infinite temperature bath. While our mechanism is not scalable to infinitely large systems, the maximum effective size can be very large, and it is fully compatible with active error correction schemes. Our model uses only on-site and nearest-neighbor interactions, and could be implemented in superconducting qubits. We sketch one such implementation at the end of this work.

Perturbative Approach to Flat Chern Bands in the Hofstadter Model

We present a perturbative approach to the study of the Hofstadter model for when the amount of flux per plaquette is close to a rational fraction. Within this approximation certain eigenstates of the system are shown to be multi-component wavefunctions that connect smoothly to the Landau levels of the continuum. The perturbative corrections to these are higher Landau level contributions that break rotational invariance and allow the perturbed states to adopt the symmetry of the lattice. In the presence of interactions, this approach allows for the calculation of generalised Haldane pseudopotentials, and in turn, the many-body properties of the system. The method is sufficiently general that it can apply to a wide variety of lattices, interactions and magnetic field strengths.

Signatures of Fractional Exclusion Statistics in the Spectroscopy of Quantum Hall Droplets

We show how spectroscopic experiments on a small Laughlin droplet of rotating bosons can directly demonstrate Haldane fractional exclusion statistics of quasihole excitations. The characteristic signatures appear in the single-particle excitation spectrum. We show that the transitions are governed by a "many-body selection rule" which allows one to relate the number of allowed transitions to the number of quasihole states on a finite geometry. We illustrate the theory with numerically exact simulations of small numbers of particles.

3- and 4-body Interactions from 2-body interactions in Spin Models: A Route to Abelian and Non-Abelian Fractional Chern Insulators

We describe a method for engineering local $k+1$ -body interactions ( $k=1,2,3$ ) from two-body couplings in spin- ${1}{2}$ systems. When implemented in certain systems with a flat single-particle band with a unit Chern number, the resulting many-body ground states are fractional Chern insulators which exhibit abelian and non-abelian anyon excitations. The most complex of these, with $k=3$ , has Fibonacci anyon excitations; our system is thus capable of universal topological quantum computation. We then demonstrate that an appropriately tuned circuit of qubits could faithfully replicate this model up to small corrections, and further, we describe the process by which one might create and manipulate non-abelian vortices in these circuits, allowing for direct control of the system's quantum information content.

Entanglement Spectrum of Composite Fermion States in Real Space

We study the entanglement spectra of many particle systems in states which are closely related to products of Slater determinants or products of permanents, or combinations of the two. Such states notably include the Laughlin and Jain composite fermion states which describe most of the observed conductance plateaus of the fractional quantum Hall effect. We identify a set of 'Entanglement Wave Functions' (EWF), for subsets of the particles, which completely describe the entanglement spectra of such product wave functions, both in real space and in particle space. A subset of the EWF for the Laughlin and Jain states can be recognized as Composite Fermion states. These states provide an exact description of the low angular momentum sectors of the real space entanglement spectrum (RSES) of these trial wave functions and a physical explanation of the branches of excitations observed in the RSES of the Jain states.

Landau Level Mixing in the Perturbative Limit

We study the effects of Landau level mixing in the limit of weak electron interaction. We use a numerical method to obtain the two- and three-body corrections to quantum Hall pseudopotentials, which are exact to lowest order in the Landau level mixing parameter. Our results are in general agreement with certain analytic results (some derived here, some derived by other authors) in the thermodynamic limit. We find that the convergence to this thermodynamic limit can be slow. This suggests that errors could occur if one tries to use pseudopotentials derived in a thermodynamic limit for numerical work on finite systems.

Multi-Particle Pseudopotentials for Multi-Component Quantum Hall Systems

The Haldane pseudopotential construction has been an extremely powerful concept in quantum Hall physics --- it not only gives a minimal description of the space of Hamiltonians but also suggests special model Hamiltonians (those where certain pseudopotential are set to zero) that may have exactly solvable ground states with interesting properties. The purpose of this paper is to generalize the pseudopotential construction to situations where interactions are N-body and where the particles may have internal degrees of freedom such as spin or valley index. Assuming a rotationally invariant Hamiltonian, the essence of the problem is to obtain a full basis of wavefunctions for N particles with fixed relative angular momentum L. This basis decomposes into representations of SU(n) with n the number of internal degrees of freedom. We give special attention to the case where the internal degree of freedom has n=2 states, which encompasses the important cases of spin-1/2 particles and quantum Hall bilayers. We also discuss in some detail the cases of spin-1 particles (n=3) and graphene (n=4, including two spin and two valley degrees of freedom).

Spin-singlet Gaffnian wave function for fractional quantum Hall systems

We characterize in detail a wave function conceivable in fractional quantum Hall systems where a spin or equivalent degree of freedom is present. This wave function combines the properties of two previously proposed quantum Hall wave functions, namely the non-Abelian spin-singlet state and the nonunitary Gaffnian wave function. This is a spin-singlet generalization of the spin-polarized Gaffnian, which we call the "spin-singlet Gaffnian" (SSG). In this paper we present evidence demonstrating that the SSG corresponds to the ground state of a certain local Hamiltonian, which we explicitly construct, and, further, we provide a relatively simple analytic expression for the unique ground-state wave functions, which we define as the zero energy eigenstates of that local Hamiltonian. In addition, we have determined a certain nonunitary, rational conformal field theory which provides an underlying description of the SSG and we thus conclude that the SSG is ungapped in the thermodynamic limit. In order to verify our construction, we implement two recently proposed techniques for the analysis of fractional quantum Hall trial states: The "spin dressed squeezing algorithm", and the "generalized Pauli principle".

Tunnelling current through fractional quantum Hall interferometers

We calculate the tunnelling current through a Fabry-Pérot interferometer in the fractional quantum Hall regime. Within linear response theory (weak tunnelling but arbitrary source-drain voltage) we find a general expression for the current due to tunnelling of quasiparticles in terms of Carlson's $R$ function. Our result is valid for fractional quantum Hall states with an edge theory consisting of a charged channel and any number of neutral channels, with possibly different edge velocities and different chiralities. We analyse the case with a single neutral channel in detail, which applies for instance to the edge of the Moore-Read state. In addition we consider an asymmetric interferometer with different edge lengths between the point contacts on opposite edges, and we study the behaviour of the current as a function of varying edge length. Recent experiments attempted to measure the Aharanov-Bohm effect by changing the area inside the interferometer using a plunger gate. Theoretical analyses of these experiments have so far not taken into account the accompanying change in the edge lengths. We show that the tunnelling current exhibits multiple osculations as a function of this edge length, with frequencies proportional to the injected edge current and inversely proportional to the edge velocities. In particular the edge velocities can be measured by looking at the Fourier spectrum of the edge current. We provide a numerical scheme to calculate and plot the $R$ function, and include sample plots for a variety of edge states with parameter values which are experimentally relevant.

Evaluation of ranks of real space and particle entanglement spectra for large systems

We devise a way to calculate the dimensions of symmetry sectors appearing in the Particle Entanglement Spectrum (PES) and Real Space Entanglement Spectrum (RSES) of multi-particle systems from their real space wave functions. We first note that these ranks in the entanglement spectra equal the dimensions of spaces of wave functions with a number of particles fixed. This also yields equality of the multiplicities in the PES and the RSES. Our technique allows numerical calculations for much larger systems than were previously feasible. For somewhat smaller systems, we can find approximate entanglement energies as well as multiplicities. We illustrate the method with results on the RSES and PES multiplicities for integer quantum Hall states, Laughlin and Jain composite fermion states and for the Moore-Read state at filling $ν=5/2$ , for system sizes up to 70 particles.

Exactly Solvable Lattice Models with Crossing Symmetry

We show how to compute the exact partition function for lattice statistical-mechanical models whose Boltzmann weights obey a special "crossing" symmetry. The crossing symmetry equates partition functions on different trivalent graphs, allowing a transformation to a graph where the partition function is easily computed. The simplest example is counting the number of nets without ends on the honeycomb lattice, including a weight per branching. Other examples include an Ising model on the Kagome' lattice with three-spin interactions, dimers on any graph of corner-sharing triangles, and non-crossing loops on the honeycomb lattice, where multiple loops on each edge are allowed. We give several methods for obtaining models with this crossing symmetry, one utilizing discrete groups and another anyon fusion rules. We also present results indicating that for models which deviate slightly from having crossing symmetry, a real-space decimation (renormalization-group-like) procedure restores the crossing symmetry.

Fractional Quantum Hall Effect of Lattice Bosons Near Commensurate Flux

We study interacting bosons on a lattice in a magnetic field. When the number of flux quanta per plaquette is close to a rational fraction, the low energy physics is mapped to a multi-species continuum model: bosons in the lowest Landau level where each boson is given an internal degree of freedom, or pseudospin. We find that the interaction potential between the bosons involves terms that do not conserve pseudospin, corresponding to umklapp processes, which in some cases can also be seen as BCS-type pairing terms. We argue that in experimentally realistic regimes for bosonic atoms in optical lattices with synthetic magnetic fields, these terms are crucial for determining the nature of allowed ground states. In particular, we show numerically that certain paired wavefunctions related to the Moore-Read Pfaffian state are stabilized by these terms, whereas certain other wavefunctions can be destabilized when umklapp processes become strong.

Importance of interband transitions for the fractional quantum Hall effect in bilayer graphene

Several recent works have proposed that electron-electron interactions in bilayer graphene can be tuned with the help of external parameters, making it possible to stabilize different fractional quantum Hall states. In these prior works, phase diagrams were calculated based on a single Landau level approximation. We go beyond this approximation and investigate the influence of polarization effects and virtual interband transitions on the stability of fractional quantum Hall states in bilayer graphene. We find that for realistic values of the dielectric constant, the phase diagram is strongly modified by these effects. We illustrate this by evaluating the region of stability of the Pfaffian state.

Spinful Composite Fermions in a Negative Effective Field

In this paper we study fractional quantum Hall composite fermion wavefunctions at filling fractions ν= 2/3, 3/5, and 4/7. At each of these filling fractions, there are several possible wavefunctions with different spin polarizations, depending on how many spin-up or spin-down composite fermion Landau levels are occupied. We calculate the energy of the possible composite fermion wavefunctions and we predict transitions between ground states of different spin polarizations as the ratio of Zeeman energy to Coulomb energy is varied. Previously, several experiments have observed such transitions between states of differing spin polarization and we make direct comparison of our predictions to these experiments. For more detailed comparison between theory and experiment, we also include finite-thickness effects in our calculations. We find reasonable qualitative agreement between the experiments and composite fermion theory. Finally, we consider composite fermion states at filling factors ν= 2+2/3, 2+3/5, and 2+4/7. The latter two cases we predict to be spin polarized even at zero Zeeman energy.

Telegraph noise and the Fabry-Perot quantum Hall interferometer

We consider signatures of abelian and nonabelian quasiparticle statistics in quantum Hall Fabry-Perot interferometers. When quasiparticles enter and exit the interference cell, for instance due to glassy motion in the dopant layer, the anyonic phase can be observed in phase jumps. In the case of the nonabelian nu=5/2 state, if the interferometer is small, we argue that free Majorana states in the interference cell are either strongly coupled to one another or are strongly coupled to the edge. We analyze the expected phase jumps and in particular suggest that changes in the parity of the Majorana ground state should gives rise to characteristic jumps of pi in the interference phase.

Three-dimensional topological lattice models with surface anyons

We study a class of three dimensional exactly solvable models of topological matter first put forward by Walker and Wang [arXiv:1104.2632v2]. While these are not models of interacting fermions, they may well capture the topological behavior of some strongly correlated systems. In this work we give a full pedagogical treatment of a special simple case of these models, which we call the 3D semion model: We calculate its ground state degeneracies for a variety of boundary conditions, and classify its low-lying excitations. While point defects in the bulk are confined in pairs connected by energetic strings, the surface excitations are more interesting: the model has deconfined point defects pinned to the boundary of the lattice, and these exhibit semionic braiding statistics. The surface physics is reminiscent of a $ν=1/2$ bosonic fractional quantum Hall effect in its topological limit, and these considerations help motivate an effective field theoretic description for the lattice models as variants of $bF$ theories. Our special example of the 3D semion model captures much of the behavior of more general `confined Walker-Wang models'. We contrast the 3D semion model with the closely related 3D version of the toric code (a lattice gauge theory) which has deconfined point excitations in the bulk and we discuss how more general models may have some confined and some deconfined excitations. Having seen that there exist lattice models whose surfaces have the same topological order as a bosonic fractional quantum Hall effect on a confining bulk, we construct a lattice model whose surface has similar topological order to a fermionic quantum hall effect. We find that in these models a fermion is always deconfined in the three dimensional bulk.

preprint2011arXiv

Condensation of achiral simple currents in topological lattice models: a Hamiltonian study of topological symmetry breaking

We describe a family of phase transitions connecting phases of differing non-trivial topological order by explicitly constructing Hamiltonians of the Levin-Wen[PRB 71, 045110] type which can be tuned between two solvable points, each of which realizes a different topologically ordered phase. We show that the low-energy degrees of freedom near the phase transition can be mapped onto those of a Potts model, and we discuss the stability of the resulting phase diagram to small perturbations about the model. We further explain how the excitations in the condensed phase are formed from those in the original topological theory, some of which are split into multiple components by condensation, and we discuss the implications of our results for understanding the nature of general achiral topological phases in 2+1 dimensions in terms of doubled Chern-Simons theories.

preprint2011arXiv

Wilson Line Picture of Levin-Wen Partition Functions

Levin and Wen [Phys. Rev. B 71, 045110 (2005)] have recently given a lattice Hamiltonian description of doubled Chern-Simons theories. We relate the partition function of these theories to an expectation of Wilson loops that form a link in 2+1 dimensional spacetime known in the mathematical literature as Chain-Mail. This geometric construction gives physical interpretation of the Levin-Wen Hilbert space and Hamiltonian, its topological invariance, exactness under coarse-graining, and how two opposite chirality sectors of the doubled theory arise.

A Random Matrix--Theoretic Approach to Handling Singular Covariance Estimates

In many practical situations we would like to estimate the covariance matrix of a set of variables from an insufficient amount of data. More specifically, if we have a set of $N$ independent, identically distributed measurements of an $M$ dimensional random vector the maximum likelihood estimate is the sample covariance matrix. Here we consider the case where $N<M$ such that this estimate is singular and therefore fundamentally bad. We present a radically new approach to deal with this situation. Let $X$ be the $M\times N$ data matrix, where the columns are the $N$ independent realizations of the random vector with covariance matrix $Σ$ . Without loss of generality, we can assume that the random variables have zero mean. We would like to estimate $Σ$ from $X$ . Let $K$ be the classical sample covariance matrix. Fix a parameter $1\leq L\leq N$ and consider an ensemble of $L\times M$ random unitary matrices, $\{Φ\}$ , having Haar probability measure. Pre and post multiply $K$ by $Φ$ , and by the conjugate transpose of $Φ$ respectively, to produce a non--singular $L\times L$ reduced dimension covariance estimate. A new estimate for $Σ$ , denoted by $\mathrm{cov}_L(K)$ , is obtained by a) projecting the reduced covariance estimate out (to $M\times M$ ) through pre and post multiplication by the conjugate transpose of $Φ$ , and by $Φ$ respectively, and b) taking the expectation over the unitary ensemble. Another new estimate (this time for $Σ^{-1}$ ), $\mathrm{invcov}_L(K)$ , is obtained by a) inverting the reduced covariance estimate, b) projecting the inverse out (to $M\times M$ ) through pre and post multiplication by the conjugate transpose of $Φ$ , and by $Φ$ respectively, and c) taking the expectation over the unitary ensemble. We have a closed analytical expression for $\mathrm{invcov}_L(K)$ and $\mathrm{cov}_L(K)$ in terms of its eigenvalue decomposition.

Phase transitions in topological lattice models via topological symmetry breaking

We study transitions between phases of matter with topological order. By studying these transitions in exactly solvable lattice models we show how universality classes may be identified and critical properties described. As a familiar example to elucidate our results concretely, we describe in detail a transition between a fully gapped achiral 2D $p$ -wave superconductor ( $p+ip$ for pseudospin up/ $p-ip$ for pseudospin down) to an $s$ -wave superconductor which we show to be in the 2D transverse field Ising universality class.

Space-Time Geometry of Topological phases

The 2+1 dimensional lattice models of Levin and Wen [PRB 71, 045110 (2005)] provide the most general known microscopic construction of topological phases of matter. Based heavily on the mathematical structure of category theory, many of the special properties of these models are not obvious. In the current paper, we present a geometrical space-time picture of the partition function of the Levin-Wen models which can be described as doubles (two copies with opposite chiralities) of underlying Anyon theories. Our space-time picture describes the partition function as a knot invariant of a complicated link, where both the lattice variables of the microscopic Levin-Wen model and the terms of the hamiltonian are represented as labeled strings of this link. This complicated link, previously studied in the mathematical literature, and known as Chain-Mail, can be related directly to known topological invariants of 3-manifolds such as the so called Turaev-Viro invariant and the Witten-Reshitikhin-Turaev invariant. We further consider quasi-particle excitations of the Levin-Wen models and we see how they can be understood by adding additional strings to the Chain-Mail link representing quasi-particle world-lines. Our construction gives particularly important new insight into how a doubled theory arises from these microscopic models.