Publications

Order by: publication date number of citations

• QuSpin: a Python package for dynamics and exact diagonalisation of quantum many body systems. Part II: bosons, fermions and higher spins

  Phillip Weinberg, Marin Bukov

  SciPost Phys. 7, 020 (2019) · published 9 August 2019 | Toggle abstract · pdf

  We present a major update to QuSpin, SciPostPhys.2.1.003 -- an open-source Python package for exact diagonalization and quantum dynamics of arbitrary boson, fermion and spin many-body systems, supporting the use of various (user-defined) symmetries in one and higher dimension and (imaginary) time evolution following a user-specified driving protocol. We explain how to use the new features of QuSpin using seven detailed examples of various complexity: (i) the transverse-field Ising chain and the Jordan-Wigner transformation, (ii) free particle systems: the Su-Schrieffer-Heeger (SSH) model, (iii) the many-body localized 1D Fermi-Hubbard model, (iv) the Bose-Hubbard model in a ladder geometry, (v) nonlinear (imaginary) time evolution and the Gross-Pitaevskii equation on a 1D lattice, (vi) integrability breaking and thermalizing dynamics in the translationally-invariant 2D transverse-field Ising model, and (vii) out-of-equilibrium Bose-Fermi mixtures. This easily accessible and user-friendly package can serve various purposes, including educational and cutting-edge experimental and theoretical research. The complete package documentation is available under http://weinbe58.github.io/QuSpin/index.html.

  176 citations

• Hyperbolic nodal band structures and knot invariants

  Marcus Stålhammar, Lukas Rødland, Gregory Arone, Jan Carl Budich, Emil J. Bergholtz

  SciPost Phys. 7, 019 (2019) · published 8 August 2019 | Toggle abstract · pdf

  We extend the list of known band structure topologies to include a large family of hyperbolic nodal links and knots, occurring both in conventional Hermitian systems where their stability relies on discrete symmetries, and in the dissipative non-Hermitian realm where the knotted nodal lines are generic and thus stable towards any small perturbation. We show that these nodal structures, taking the forms of Turk's head knots, appear in both continuum- and lattice models with relatively short-ranged hopping that is within experimental reach. To determine the topology of the nodal structures, we devise an efficient algorithm for computing the Alexander polynomial, linking numbers and higher order Milnor invariants based on an approximate and well controlled parameterisation of the knot.

  15 citations

• Exact real-time dynamics of single-impurity Anderson model from a single-spin hybridization-expansion

  Patryk Kubiczek, Alexey N. Rubtsov, Alexander I. Lichtenstein

  SciPost Phys. 7, 016 (2019) · published 2 August 2019 | Toggle abstract · pdf

  In this work we introduce a modified real-time continuous-time hybridization-expansion quantum Monte Carlo solver for a time-dependent single-orbital Anderson impurity model: CT-1/2-HYB-QMC. In the proposed method the diagrammatic expansion is performed only for one out of the two spin channels, while the resulting effective single-particle problem for the other spin is solved semi-analytically for each expansion diagram. CT-1/2-HYB-QMC alleviates the dynamical sign problem by reducing the order of sampled diagrams and makes it possible to reach twice as long time scales in comparison to the standard CT-HYB method. We illustrate the new solver by calculating an electric current through impurity in paramagnetic and spin-polarized cases.

  5 citations

• Out of time ordered effective dynamics of a quartic oscillator

  Bidisha Chakrabarty, Soumyadeep Chaudhuri

  SciPost Phys. 7, 013 (2019) · published 29 July 2019 | Toggle abstract · pdf

  We study the dynamics of a quantum Brownian particle weakly coupled to a thermal bath. Working in the Schwinger-Keldysh formalism, we develop an effective action of the particle up to quartic terms. We demonstrate that this quartic effective theory is dual to a stochastic dynamics governed by a non-linear Langevin equation. The Schwinger-Keldysh effective theory, or the equivalent non-linear Langevin dynamics, is insufficient to determine the out of time order correlators (OTOCs) of the particle. To overcome this limitation, we construct an extended effective action in a generalised Schwinger-Keldysh framework. We determine the additional quartic couplings in this OTO effective action and show their dependence on the bath's 4-point OTOCs. We analyse the constraints imposed on the OTO effective theory by microscopic reversibility and thermality of the bath. We show that these constraints lead to a generalised fluctuation-dissipation relation between the non-Gaussianity in the distribution of the thermal noise experienced by the particle and the thermal jitter in its damping coefficient. The quartic effective theory developed in this work provides extension of several results previously obtained for the cubic OTO dynamics of a Brownian particle.

  7 citations

• Density-operator evolution: Complete positivity and the Keldysh real-time expansion

  Viktor Reimer, Maarten Rolf Wegewijs

  SciPost Phys. 7, 012 (2019) · published 25 July 2019 | Toggle abstract · pdf

  We study the reduced time-evolution of open quantum systems by combining quantum-information and statistical field theory. Inspired by prior work [EPL 102, 60001 (2013) and Phys. Rev. Lett. 111, 050402 (2013)] we establish the explicit structure guaranteeing the complete positivity (CP) and trace-preservation (TP) of the real-time evolution expansion in terms of the microscopic system-environment coupling. This reveals a fundamental two-stage structure of the coupling expansion: Whereas the first stage defines the dissipative timescales of the system --before having integrated out the environment completely-- the second stage sums up elementary physical processes described by CP superoperators. This allows us to establish the nontrivial relation between the (Nakajima-Zwanzig) memory-kernel superoperator for the density operator and novel memory-kernel operators that generate the Kraus operators of an operator-sum. Importantly, this operational approach can be implemented in the existing Keldysh real-time technique and allows approximations for general time-nonlocal quantum master equations to be systematically compared and developed while keeping the CP and TP structure explicit. Our considerations build on the result that a Kraus operator for a physical measurement process on the environment can be obtained by 'cutting' a group of Keldysh real-time diagrams 'in half'. This naturally leads to Kraus operators lifted to the system plus environment which have a diagrammatic expansion in terms of time-nonlocal memory-kernel operators. These lifted Kraus operators obey coupled time-evolution equations which constitute an unraveling of the original Schr\"odinger equation for system plus environment. Whereas both equations lead to the same reduced dynamics, only the former explicitly encodes the operator-sum structure of the coupling expansion.

  8 citations

• Bootstrapping MN and tetragonal CFTs in three dimensions

  Andreas Stergiou

  SciPost Phys. 7, 010 (2019) · published 17 July 2019 | Toggle abstract · pdf

  Conformal field theories (CFTs) with MN and tetragonal global symmetry in $d=2+1$ dimensions are relevant for structural, antiferromagnetic and helimagnetic phase transitions. As a result, they have been studied in great detail with the $\varepsilon=4-d$ expansion and other field theory methods. The study of these theories with the nonperturbative numerical conformal bootstrap is initiated in this work. Bounds for operator dimensions are obtained and they are found to possess sharp kinks in the MN case, suggesting the existence of full-fledged CFTs. Based on the existence of a certain large-$N$ expansion in theories with MN symmetry, these are argued to be the CFTs predicted by the $\varepsilon$ expansion. In the tetragonal case no new kinks are found, consistently with the absence of such CFTs in the $\varepsilon$ expansion. Estimates for critical exponents are provided for a few cases describing phase transitions in actual physical systems. In two particular MN cases, corresponding to theories with global symmetry groups $O(2)^2\rtimes S_2$ and $O(2)^3\rtimes S_3$, a second kink is found. In the $O(2)^2\rtimes S_2$ case it is argued to be saturated by a CFT that belongs to a new universality class relevant for the structural phase transition of NbO$_2$ and paramagnetic-helimagnetic transitions of the rare-earth metals Ho and Dy. In the $O(2)^3\rtimes S_3$ case it is suggested that the CFT that saturates the second kink belongs to a new universality class relevant for the paramagnetic-antiferromagnetic phase transition of the rare-earth metal Nd.

  20 citations

• A web of 2d dualities: ${\bf Z}_2$ gauge fields and Arf invariants

  Andreas Karch, David Tong, Carl Turner

  SciPost Phys. 7, 007 (2019) · published 10 July 2019 | Toggle abstract · pdf

  We describe a web of well-known dualities connecting quantum field theories in $d=1+1$ dimensions. The web is constructed by gauging ${\bf Z}_2$ global symmetries and includes a number of perennial favourites such as the Jordan-Wigner transformation, Kramers-Wannier duality, bosonization of a Dirac fermion, and T-duality. There are also less-loved examples, such as non-modular invariant $c=1$ CFTs that depend on a background spin structure.

  69 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

• Entanglement evolution and generalised hydrodynamics: interacting integrable systems

  Vincenzo Alba, Bruno Bertini, Maurizio Fagotti

  SciPost Phys. 7, 005 (2019) · published 8 July 2019 | Toggle abstract · pdf

  We investigate the dynamics of bipartite entanglement after the sudden junction of two leads in interacting integrable models. By combining the quasiparticle picture for the entanglement spreading with Generalised Hydrodynamics we derive an analytical prediction for the dynamics of the entanglement entropy between a finite subsystem and the rest. We find that the entanglement rate between the two leads depends only on the physics at the interface and differs from the rate of exchange of thermodynamic entropy. This contrasts with the behaviour in free or homogeneous interacting integrable systems, where the two rates coincide.

  60 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