Publications

Order by: publication date number of citations

• Role of Sb in the superconducting kagome metal CsV$_3$Sb$_5$ revealed by its anisotropic compression

  Alexander A. Tsirlin, Pierre Fertey, Brenden R. Ortiz, Berina Klis, Valentino Merkl, Martin Dressel, Stephen D. Wilson, Ece Uykur

  SciPost Phys. 12, 049 (2022) · published 2 February 2022 | Toggle abstract · pdf

  Pressure evolution of the superconducting kagome metal CsV$_3$Sb$_5$ is studied with single-crystal x-ray diffraction and density-functional band-structure calculations. A highly anisotropic compression observed up to 5 GPa is ascribed to the fast shrinkage of the Cs-Sb distances and suppression of Cs rattling motion. This prevents Sb displacements required to stabilize the three-dimensional charge-density-wave (CDW) order and elucidates the disappearance of the CDW already at 2 GPa despite only minor changes in the electronic structure of the normal state. At higher pressures, vanadium bands still change only marginally, whereas antimony bands undergo a major reconstruction caused by the gradual formation of the interlayer Sb-Sb bonds. Our results exclude pressure tuning of vanadium kagome bands as the main mechanism for the non-trivial evolution of superconductivity in real-world kagome metals. Concurrently, we establish the central role of Sb atoms in the stabilization of a three-dimensional CDW and Fermi surface reconstruction.

  30 citations

• Evaluation of time-dependent correlators after a local quench in iPEPS: hole motion in the t-J model

  Claudius Hubig, Annabelle Bohrdt, Michael Knap, Fabian Grusdt, J. Ignacio Cirac

  SciPost Phys. 8, 021 (2020) · published 6 February 2020 | Toggle abstract · pdf

  Infinite projected entangled pair states (iPEPS) provide a convenient variational description of infinite, translationally-invariant two-dimensional quantum states. However, the simulation of local excitations is not directly possible due to the translationally-invariant ansatz. Furthermore, as iPEPS are either identical or orthogonal, expectation values between different states as required during the evaluation of non-equal-time correlators are ill-defined. Here, we show that by introducing auxiliary states on each site, it becomes possible to simulate both local excitations and evaluate non-equal-time correlators in an iPEPS setting under real-time evolution. We showcase the method by simulating the t-J model after a single hole has been placed in the half-filled antiferromagnetic background and evaluating both return probabilities and spin correlation functions, as accessible in quantum gas microscopes.

  28 citations

• Investigating ultrafast quantum magnetism with machine learning

  Giammarco Fabiani, Johan H. Mentink

  SciPost Phys. 7, 004 (2019) · published 5 July 2019 | Toggle abstract · pdf

  We investigate the efficiency of the recently proposed Restricted Boltzmann Machine (RBM) representation of quantum many-body states to study both the static properties and quantum spin dynamics in the two-dimensional Heisenberg model on a square lattice. For static properties we find close agreement with numerically exact Quantum Monte Carlo results in the thermodynamical limit. For dynamics and small systems, we find excellent agreement with exact diagonalization, while for systems up to N=256 spins close consistency with interacting spin-wave theory is obtained. In all cases the accuracy converges fast with the number of network parameters, giving access to much bigger systems than feasible before. This suggests great potential to investigate the quantum many-body dynamics of large scale spin systems relevant for the description of magnetic materials strongly out of equilibrium.

  25 citations

• The vacua of dipolar cavity quantum electrodynamics

  Michael Schuler, Daniele De Bernardis, Andreas M. Läuchli, Peter Rabl

  SciPost Phys. 9, 066 (2020) · published 9 November 2020 | Toggle abstract · pdf

  The structure of solids and their phases is mainly determined by static Coulomb forces while the coupling of charges to the dynamical, i.e., quantized degrees of freedom of the electromagnetic field plays only a secondary role. Recently, it has been speculated that this general rule can be overcome in the context of cavity quantum electrodynamics (QED), where the coupling of dipoles to a single field mode can be dramatically enhanced. Here we present a first exact analysis of the ground states of a dipolar cavity QED system in the non-perturbative coupling regime, where electrostatic and dynamical interactions play an equally important role. Specifically, we show how strong and long-range vacuum fluctuations modify the states of dipolar matter and induce novel phases with unusual properties. Beyond a purely fundamental interest, these general mechanisms can be important for potential applications, ranging from cavity-assisted chemistry to quantum technologies based on ultrastrongly coupled circuit QED systems.

  24 citations

• The ALF (Algorithms for Lattice Fermions) project release 2.0. Documentation for the auxiliary-field quantum Monte Carlo code

  F. F. Assaad, M. Bercx, F. Goth, A. Götz, J. S. Hofmann, E. Huffman, Z. Liu, F. Parisen Toldin, J. S. E. Portela, J. Schwab

  SciPost Phys. Codebases 1 (2022) · published 22 August 2022 · licensed under CC BY-SA (4.0) | Toggle abstract · pdf

  The Algorithms for Lattice Fermions package provides a general code for the finite-temperature and projective auxiliary-field quantum Monte Carlo algorithm. The code is engineered to be able to simulate any model that can be written in terms of sums of single-body operators, of squares of single-body operators and single-body operators coupled to a bosonic field with given dynamics. The package includes five pre-defined model classes: SU(N) Kondo, SU(N) Hubbard, SU(N) t-V and SU(N) models with long range Coulomb repulsion on honeycomb, square and N-leg lattices, as well as $Z_2$ unconstrained lattice gauge theories coupled to fermionic and $Z_2$ matter. An implementation of the stochastic Maximum Entropy method is also provided. One can download the code from our Git instance at https://git.physik.uni-wuerzburg.de/ALF/ALF/-/tree/ALF-2.0 and sign in to file issues.

  23 citations

• Fast counting with tensor networks

  Stefanos Kourtis, Claudio Chamon, Eduardo R. Mucciolo, Andrei E. Ruckenstein

  SciPost Phys. 7, 060 (2019) · published 12 November 2019 | Toggle abstract · pdf

  We introduce tensor network contraction algorithms for counting satisfying assignments of constraint satisfaction problems (#CSPs). We represent each arbitrary #CSP formula as a tensor network, whose full contraction yields the number of satisfying assignments of that formula, and use graph theoretical methods to determine favorable orders of contraction. We employ our heuristics for the solution of #P-hard counting boolean satisfiability (#SAT) problems, namely monotone #1-in-3SAT and #Cubic-Vertex-Cover, and find that they outperform state-of-the-art solvers by a significant margin.

  23 citations

• Scars from protected zero modes and beyond in $U(1)$ quantum link and quantum dimer models

  Saptarshi Biswas, Debasish Banerjee, Arnab Sen

  SciPost Phys. 12, 148 (2022) · published 6 May 2022 | Toggle abstract · pdf

  We demonstrate the presence of anomalous high-energy eigenstates, or many-body scars, in $U(1)$ quantum link and quantum dimer models on square and rectangular lattices. In particular, we consider the paradigmatic Rokhsar-Kivelson Hamiltonian $H=\mathcal{O}_{\mathrm{kin}} + \lambda \mathcal{O}_{\mathrm{pot}}$ where $\mathcal{O}_{\mathrm{pot}}$ ($\mathcal{O}_{\mathrm{kin}}$) is defined as a sum of terms on elementary plaquettes that are diagonal (off-diagonal) in the computational basis. Both these interacting models possess an exponentially large number of mid-spectrum zero modes in system size at $\lambda=0$ that are protected by an index theorem preventing any mixing with the nonzero modes at this coupling. We classify different types of scars for $|\lambda| \lesssim \mathcal{O}(1)$ both at zero and finite winding number sectors complementing and significantly generalizing our previous work [Banerjee and Sen, Phys. Rev. Lett. 126, 220601 (2021)]. The scars at finite $\lambda$ show a rich variety with those that are composed solely from the zero modes of $\mathcal{O}_{\mathrm{kin}}$, those that contain an admixture of both the zero and the nonzero modes of $\mathcal{O}_{\mathrm{kin}}$, and finally those composed solely from the nonzero modes of $\mathcal{O}_{\mathrm{kin}}$. We give analytic expressions for certain "lego scars" for the quantum dimer model on rectangular lattices where one of the linear dimensions can be made arbitrarily large, with the building blocks (legos) being composed of emergent singlets and other more complicated entangled structures.

  22 citations

• Spectral response of Josephson junctions with low-energy quasiparticles

  Anna Keselman, Chaitanya Murthy, Bernard van Heck, Bela Bauer

  SciPost Phys. 7, 050 (2019) · published 16 October 2019 | Toggle abstract · pdf

  We study nanowire-based Josephson junctions shunted by a capacitor and take into account the presence of low-energy quasiparticle excitations. These are treated by extending conventional models used to describe superconducting qubits to include the coherent coupling between fermionic quasiparticles, in particular the Majorana zero modes that emerge in topological superconductors, and the plasma mode of the junction. Using accurate, unbiased matrix-product state techniques, we compute the energy spectrum and response function of the system across the topological phase transition. Furthermore, we develop a perturbative approach, valid in the harmonic limit with small charging energy, illustrating how the presence of low-energy quasiparticles affects the spectrum and response of the junction. Our results are of direct interest to on-going experimental investigations of nanowire-based superconducting qubits.

  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.

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

  20 citations