Publications

Order by: publication date number of 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

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

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

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

  18 citations

• Mean-field model of interacting quasilocalized excitations in glasses

  Corrado Rainone, Pierfrancesco Urbani, Francesco Zamponi, Edan Lerner, and Eran Bouchbinder

  SciPost Phys. Core 4, 008 (2021) · published 16 April 2021 | Toggle abstract · pdf

  Structural glasses feature quasilocalized excitations whose frequencies $\omega$ follow a universal density of states ${\cal D}(\omega)\!\sim\!\omega^4$. Yet, the underlying physics behind this universality is not yet fully understood. Here we study a mean-field model of quasilocalized excitations in glasses, viewed as groups of particles embedded inside an elastic medium and described collectively as anharmonic oscillators. The oscillators, whose harmonic stiffness is taken from a rather featureless probability distribution (of upper cutoff $\kappa_0$) in the absence of interactions, interact among themselves through random couplings (characterized by strength $J$) and with the surrounding elastic medium (an interaction characterized by a constant force $h$). We first show that the model gives rise to a gapless density of states ${\cal D}(\omega)\!=\!A_{\rm g}\,\omega^4$ for a broad range of model parameters, expressed in terms of the strength of stabilizing anharmonicity, which plays a decisive role in the model. Then --- using scaling theory and numerical simulations --- we provide a complete understanding of the non-universal prefactor $A_{\rm g}(h,J,\kappa_0)$, of the oscillators' interaction-induced mean square displacement and of an emerging characteristic frequency, all in terms of properly identified dimensionless quantities. In particular, we show that $A_{\rm g}(h,J,\kappa_0)$ is a non-monotonic function of $J$ for a fixed $h$, varying predominantly exponentially with $-(\kappa_0 h^{2/3}\!/J^2)$ in the weak interactions (small $J$) regime --- reminiscent of recent observations in computer glasses --- and predominantly decays as a power-law for larger $J$, in a regime where $h$ plays no role. We discuss the physical interpretation of the model and its possible relations to available observations in structural glasses, along with delineating some future research directions.

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

  18 citations

• Quantitative functional renormalization for three-dimensional quantum Heisenberg models

  Nils Niggemann, Johannes Reuther, Björn Sbierski

  SciPost Phys. 12, 156 (2022) · published 11 May 2022 | Toggle abstract · pdf

  We employ a recently developed variant of the functional renormalization group method for spin systems, the so-called pseudo Majorana functional renormalization group, to investigate three-dimensional spin-1/2 Heisenberg models at finite temperatures. We study unfrustrated and frustrated Heisenberg systems on the simple cubic and pyrochlore lattices. Comparing our results with other quantum many-body techniques, we demonstrate a high quantitative accuracy of our method. Particularly, for the unfrustrated simple cubic lattice antiferromagnet ordering temperatures obtained from finite-size scaling of one-loop data deviate from error controlled quantum Monte Carlo results by $\sim5\%$ and we further confirm the established values for the critical exponent $\nu$ and the anomalous dimension $\eta$. As the PMFRG yields results in good agreement with QMC, but remains applicable when the system is frustrated, we next treat the pyrochlore Heisenberg antiferromagnet as a paradigmatic magnetically disordered system and find nearly perfect agreement of our two-loop static homogeneous susceptibility with other methods. We further investigate the broadening of pinch points in the spin structure factor as a result of quantum and thermal fluctuations and confirm a finite width in the extrapolated limit $T\rightarrow0$. While extensions towards higher loop orders $\ell$ seem to systematically improve our approach for magnetically disordered systems we also discuss subtleties when increasing $\ell$ in the presence of magnetic order. Overall, the pseudo Majorana functional renormalization group is established as a powerful many-body technique in quantum magnetism with a wealth of possible future applications.

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

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

  17 citations

• Impurities in a one-dimensional Bose gas: the flow equation approach

  Fabian Brauneis, Hans-Werner Hammer, Mikhail Lemeshko, Artem G. Volosniev

  SciPost Phys. 11, 008 (2021) · published 13 July 2021 | Toggle abstract · pdf

  A few years ago, flow equations were introduced as a technique for calculating the ground-state energies of cold Bose gases with and without impurities. In this paper, we extend this approach to compute observables other than the energy. As an example, we calculate the densities, and phase fluctuations of one-dimensional Bose gases with one and two impurities. For a single mobile impurity, we use flow equations to validate the mean-field results obtained upon the Lee-Low-Pines transformation. We show that the mean-field approximation is accurate for all values of the boson-impurity interaction strength as long as the phase coherence length is much larger than the healing length of the condensate. For two static impurities, we calculate impurity-impurity interactions induced by the Bose gas. We find that leading order perturbation theory fails when boson-impurity interactions are stronger than boson-boson interactions. The mean-field approximation reproduces the flow equation results for all values of the boson-impurity interaction strength as long as boson-boson interactions are weak.

  17 citations