Publications

Order by: publication date number of citations

• Amorphous topological phases protected by continuous rotation symmetry

  Helene Spring, Anton R. Akhmerov, Daniel Varjas

  SciPost Phys. 11, 022 (2021) · published 5 August 2021 | Toggle abstract · pdf

  Protection of topological surface states by reflection symmetry breaks down when the boundary of the sample is misaligned with one of the high symmetry planes of the crystal. We demonstrate that this limitation is removed in amorphous topological materials, where the Hamiltonian is invariant on average under reflection over any axis due to continuous rotation symmetry. While the local disorder caused by the amorphous structure weakens the topological protection, we demonstrate that the edge remains protected from localization. In order to classify such phases we perform a systematic search over all the possible symmetry classes in two dimensions and construct the example models realizing each of the proposed topological phases. Finally, we compute the topological invariant of these phases as an integral along a meridian of the spherical Brillouin zone of an amorphous Hamiltonian.

  21 citations

• Low rank compression in the numerical solution of the nonequilibrium Dyson equation

  Jason Kaye, Denis Golež

  SciPost Phys. 10, 091 (2021) · published 26 April 2021 | Toggle abstract · pdf

  We propose a method to improve the computational and memory efficiency of numerical solvers for the nonequilibrium Dyson equation in the Keldysh formalism. It is based on the empirical observation that the nonequilibrium Green's functions and self energies arising in many problems of physical interest, discretized as matrices, have low rank off-diagonal blocks, and can therefore be compressed using a hierarchical low rank data structure. We describe an efficient algorithm to build this compressed representation on the fly during the course of time stepping, and use the representation to reduce the cost of computing history integrals, which is the main computational bottleneck. For systems with the hierarchical low rank property, our method reduces the computational complexity of solving the nonequilibrium Dyson equation from cubic to near quadratic, and the memory complexity from quadratic to near linear. We demonstrate the full solver for the Falicov-Kimball model exposed to a rapid ramp and Floquet driving of system parameters, and are able to increase feasible propagation times substantially. We present examples with 262144 time steps, which would require approximately five months of computing time and 2.2 TB of memory using the direct time stepping method, but can be completed in just over a day on a laptop with less than 4 GB of memory using our method. We also confirm the hierarchical low rank property for the driven Hubbard model in the weak coupling regime within the GW approximation, and in the strong coupling regime within dynamical mean-field theory.

  21 citations

• Efficient and flexible approach to simulate low-dimensional quantum lattice models with large local Hilbert spaces

  Thomas Köhler, Jan Stolpp, Sebastian Paeckel

  SciPost Phys. 10, 058 (2021) · published 9 March 2021 | Toggle abstract · pdf

  Quantum lattice models with large local Hilbert spaces emerge across various fields in quantum many-body physics. Problems such as the interplay between fermions and phonons, the BCS-BEC crossover of interacting bosons, or decoherence in quantum simulators have been extensively studied both theoretically and experimentally. In recent years, tensor network methods have become one of the most successful tools to treat such lattice systems numerically. Nevertheless, systems with large local Hilbert spaces remain challenging. Here, we introduce a mapping that allows to construct artificial $U(1)$ symmetries for any type of lattice model. Exploiting the generated symmetries, numerical expenses that are related to the local degrees of freedom decrease significantly. This allows for an efficient treatment of systems with large local dimensions. Further exploring this mapping, we reveal an intimate connection between the Schmidt values of the corresponding matrix\hyp product\hyp state representation and the single\hyp site reduced density matrix. Our findings motivate an intuitive physical picture of the truncations occurring in typical algorithms and we give bounds on the numerical complexity in comparison to standard methods that do not exploit such artificial symmetries. We demonstrate this new mapping, provide an implementation recipe for an existing code, and perform example calculations for the Holstein model at half filling. We studied systems with a very large number of lattice sites up to $L=501$ while accounting for $N_{\rm ph}=63$ phonons per site with high precision in the CDW phase.

  20 citations

• Automatic differentiation applied to excitations with projected entangled pair states

  Boris Ponsioen, Fakher F. Assaad, Philippe Corboz

  SciPost Phys. 12, 006 (2022) · published 6 January 2022 | Toggle abstract · pdf

  The excitation ansatz for tensor networks is a powerful tool for simulating the low-lying quasiparticle excitations above ground states of strongly correlated quantum many-body systems. Recently, the two-dimensional tensor network class of infinite entangled pair states gained new ground state optimization methods based on automatic differentiation, which are at the same time highly accurate and simple to implement. Naturally, the question arises whether these new ideas can also be used to optimize the excitation ansatz, which has recently been implemented in two dimensions as well. In this paper, we describe a straightforward way to reimplement the framework for excitations using automatic differentiation, and demonstrate its performance for the Hubbard model at half filling.

  20 citations

• Topological entanglement properties of disconnected partitions in the Su-Schrieffer-Heeger model

  Tommaso Micallo, Vittorio Vitale, Marcello Dalmonte, Pierre Fromholz

  SciPost Phys. Core 3, 012 (2020) · published 1 December 2020 | Toggle abstract · pdf

  We study the disconnected entanglement entropy, $S^D$, of the Su-Schrieffer-Heeger model. $S^D$ is a combination of both connected and disconnected bipartite entanglement entropies that removes all area and volume law contributions, and is thus only sensitive to the non-local entanglement stored within the ground state manifold. Using analytical and numerical computations, we show that $S^D$ behaves as a topological invariant, i.e., it is quantized to either $0$ or $2 \log (2)$ in the topologically trivial and non-trivial phases, respectively. These results also hold in the presence of symmetry-preserving disorder. At the second-order phase transition separating the two phases, $S^D$ displays a system-size scaling behavior akin to those of conventional order parameters, that allows us to compute entanglement critical exponents. To corroborate the topological origin of the quantized values of $S^D$, we show how the latter remain quantized after applying unitary time evolution in the form of a quantum quench, a characteristic feature of topological invariants.

  20 citations

• A non-Abelian parton state for the $ν=2+3/8$ fractional quantum Hall effect

  Ajit Coimbatore Balram

  SciPost Phys. 10, 083 (2021) · published 15 April 2021 | Toggle abstract · pdf

  Fascinating structures have arisen from the study of the fractional quantum Hall effect (FQHE) at the even denominator fraction of $5/2$. We consider the FQHE at another even denominator fraction, namely $\nu=2+3/8$, where a well-developed and quantized Hall plateau has been observed in experiments. We examine the non-Abelian state described by the "$\bar{3}\bar{2}^{2}1^{4}$" parton wave function and numerically demonstrate it to be a feasible candidate for the ground state at $\nu=2+3/8$. We make predictions for experimentally measurable properties of the $\bar{3}\bar{2}^{2}1^{4}$ state that can reveal its underlying topological structure.

  20 citations

• Competing Spin Liquid Phases in the S=$\frac{1}{2}$ Heisenberg Model on the Kagome Lattice

  Shenghan Jiang, Panjin Kim, Jung Hoon Han, Ying Ran

  SciPost Phys. 7, 006 (2019) · published 9 July 2019 | Toggle abstract · pdf

  The properties of ground state of spin-$\frac{1}{2}$ kagome antiferromagnetic Heisenberg (KAFH) model have attracted considerable interest in the past few decades, and recent numerical simulations reported a spin liquid phase. The nature of the spin liquid phase remains unclear. For instance, the interplay between symmetries and $Z_2$ topological order leads to different types of $Z_2$ spin liquid phases. In this paper, we develop a numerical simulation method based on symmetric projected entangled-pair states (PEPS), which is generally applicable to strongly correlated model systems in two spatial dimensions. We then apply this method to study the nature of the ground state of the KAFH model. Our results are consistent with that the ground state is a $U(1)$ Dirac spin liquid rather than a $Z_2$ spin liquid.

  19 citations

• Metallic and insulating stripes and their relation with superconductivity in the doped Hubbard model

  Luca F. Tocchio, Arianna Montorsi, Federico Becca

  SciPost Phys. 7, 021 (2019) · published 12 August 2019 | Toggle abstract · pdf

  The dualism between superconductivity and charge/spin modulations (the so-called stripes) dominates the phase diagram of many strongly-correlated systems. A prominent example is given by the Hubbard model, where these phases compete and possibly coexist in a wide regime of electron dopings for both weak and strong couplings. Here, we investigate this antagonism within a variational approach that is based upon Jastrow-Slater wave functions, including backflow correlations, which can be treated within a quantum Monte Carlo procedure. We focus on clusters having a ladder geometry with $M$ legs (with $M$ ranging from $2$ to $10$) and a relatively large number of rungs, thus allowing us a detailed analysis in terms of the stripe length. We find that stripe order with periodicity $\lambda=8$ in the charge and $2\lambda=16$ in the spin can be stabilized at doping $\delta=1/8$. Here, there are no sizable superconducting correlations and the ground state has an insulating character. A similar situation, with $\lambda=6$, appears at $\delta=1/6$. Instead, for smaller values of dopings, stripes can be still stabilized, but they are weakly metallic at $\delta=1/12$ and metallic with strong superconducting correlations at $\delta=1/10$, as well as for intermediate (incommensurate) dopings. Remarkably, we observe that spin modulation plays a major role in stripe formation, since it is crucial to obtain a stable striped state upon optimization. The relevance of our calculations for previous density-matrix renormalization group results and for the two-dimensional case is also discussed.

  19 citations

• Fermionic Gaussian states: an introduction to numerical approaches

  Jacopo Surace, Luca Tagliacozzo

  SciPost Phys. Lect. Notes 54 (2022) · published 16 May 2022 | Toggle abstract · pdf

  This document is meant to be a practical introduction to the analytical and numerical manipulation of Fermionic Gaussian systems. Starting from the basics, we move to relevant modern results and techniques, presenting numerical examples and studying relevant Hamiltonians, such as the transverse field Ising Hamiltonian, in detail. We finish introducing novel algorithms connecting Fermionic Guassian states with matrix product states techniques. All the numerical examples make use of the free Julia package F_utilities.

  19 citations

• Light-induced topological magnons in two-dimensional van der Waals magnets

  Emil Viñas Boström, Martin Claassen, James W. McIver, Gregor Jotzu, Angel Rubio, Michael A. Sentef

  SciPost Phys. 9, 061 (2020) · published 30 October 2020 | Toggle abstract · pdf

  Driving a two-dimensional Mott insulator with circularly polarized light breaks time-reversal and inversion symmetry, which induces an optically-tunable synthetic scalar spin chirality interaction in the effective low-energy spin Hamiltonian. Here, we show that this mechanism can stabilize topological magnon excitations in honeycomb ferromagnets and in optical lattices. We find that the irradiated quantum magnet is described by a Haldane model for magnons that hosts topologically-protected edge modes. We study the evolution of the magnon spectrum in the Floquet regime and via time propagation of the magnon Hamiltonian for a slowly varying pulse envelope. Compared to similar but conceptually distinct driving schemes based on the Aharanov-Casher effect, the dimensionless light-matter coupling parameter $\lambda = eEa/\hbar\omega$ at fixed electric field strength is enhanced by a factor $\sim 10^5$. This increase of the coupling parameter allows to induce a topological gap of the order of $\Delta \approx 2$ meV with realistic laser pulses, bringing an experimental realization of light-induced topological magnon edge states within reach.

  19 citations