Researcher profile

Steven H. Simon

Steven H. Simon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2022arXiv

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.

preprint2022arXiv

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.

preprint2022arXiv

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.

preprint2022arXiv

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.

preprint2022arXiv

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.

preprint2020arXiv

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.

preprint2020arXiv

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.

preprint2020arXiv

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.

preprint2020arXiv

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.

preprint2019arXiv

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.

preprint2019arXiv

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.

preprint2019arXiv

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)].

preprint2019arXiv

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.

preprint2019arXiv

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.