Publications

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

• Matrix product operator symmetries and intertwiners in string-nets with domain walls

  Laurens Lootens, Jürgen Fuchs, Jutho Haegeman, Christoph Schweigert, Frank Verstraete

  SciPost Phys. 10, 053 (2021) · published 1 March 2021 | Toggle abstract · pdf

  We provide a description of virtual non-local matrix product operator (MPO) symmetries in projected entangled pair state (PEPS) representations of string-net models. Given such a PEPS representation, we show that the consistency conditions of its MPO symmetries amount to a set of six coupled equations that can be identified with the pentagon equations of a bimodule category. This allows us to classify all equivalent PEPS representations and build MPO intertwiners between them, synthesising and generalising the wide variety of tensor network representations of topological phases. Furthermore, we use this generalisation to build explicit PEPS realisations of domain walls between different topological phases as constructed by Kitaev and Kong [Commun. Math. Phys. 313 (2012) 351-373]. While the prevailing abstract categorical approach is sufficient to describe the structure of topological phases, explicit tensor network representations are required to simulate these systems on a computer, such as needed for calculating thresholds of quantum error-correcting codes based on string-nets with boundaries. Finally, we show that all these string-net PEPS representations can be understood as specific instances of Turaev-Viro state-sum models of topological field theory on three-manifolds with a physical boundary, thereby putting these tensor network constructions on a mathematically rigorous footing.

  27 citations

• Floquet conformal field theories with generally deformed Hamiltonians

  Ruihua Fan, Yingfei Gu, Ashvin Vishwanath, Xueda Wen

  SciPost Phys. 10, 049 (2021) · published 25 February 2021 | Toggle abstract · pdf

  In this work, we study non-equilibrium dynamics in Floquet conformal field theories (CFTs) in 1+1D, in which the driving Hamiltonian involves the energy-momentum density spatially modulated by an arbitrary smooth function. This generalizes earlier work which was restricted to the sine-square deformed type of Floquet Hamiltonians, operating within a $\mathfrak{sl}_2$ sub-algebra. Here we show remarkably that the problem remains soluble in this generalized case which involves the full Virasoro algebra, based on a geometrical approach. It is found that the phase diagram is determined by the stroboscopic trajectories of operator evolution. The presence/absence of spatial fixed points in the operator evolution indicates that the driven CFT is in a heating/non-heating phase, in which the entanglement entropy grows/oscillates in time. Additionally, the heating regime is further subdivided into a multitude of phases, with different entanglement patterns and spatial distribution of energy-momentum density, which are characterized by the number of spatial fixed points. Phase transitions between these different heating phases can be achieved simply by changing the duration of application of the driving Hamiltonian. We demonstrate the general features with concrete CFT examples and compare the results to lattice calculations and find remarkable agreement.

  16 citations

• Phase-space methods for simulating the dissipative many-body dynamics of collective spin systems

  Julian Huber, Peter Kirton, Peter Rabl

  SciPost Phys. 10, 045 (2021) · published 22 February 2021 | Toggle abstract · pdf

  We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses an extension of the truncated Wigner approximation to map the exact open system dynamics onto stochastic differential equations for the corresponding phase space distribution. This approach is most effective in the limit of very large spin quantum numbers, where exact numerical simulations and other approximation methods are no longer applicable. We benchmark this numerical technique for known superradiant decay and spin-squeezing processes and illustrate its application for the simulation of non-equilibrium phase transitions in dissipative spin lattice models.

  16 citations

• Probing many-body localization in a disordered quantum dimer model on the honeycomb lattice

  Francesca Pietracaprina, Fabien Alet

  SciPost Phys. 10, 044 (2021) · published 22 February 2021 | Toggle abstract · pdf

  We numerically study the possibility of many-body localization transition in a disordered quantum dimer model on the honeycomb lattice. By using the peculiar constraints of this model and state-of-the-art exact diagonalization and time evolution methods, we probe both eigenstates and dynamical properties and conclude on the existence of a localization transition, on the available time and length scales (system sizes of up to N=108 sites). We critically discuss these results and their implications.

  13 citations

• Tangent-space methods for truncating uniform MPS

  Bram Vanhecke, Maarten Van Damme, Jutho Haegeman, Laurens Vanderstraeten, Frank Verstraete

  SciPost Phys. Core 4, 004 (2021) · published 19 February 2021 | Toggle abstract · pdf

  A central primitive in quantum tensor network simulations is the problem of approximating a matrix product state with one of a lower bond dimension. This problem forms the central bottleneck in algorithms for time evolution and for contracting projected entangled pair states. We formulate a tangent-space based variational algorithm to achieve this for uniform (infinite) matrix product states. The algorithm exhibits a favourable scaling of the computational cost, and we demonstrate its usefulness by several examples involving the multiplication of a matrix product state with a matrix product operator.

  8 citations

• Riemannian optimization of isometric tensor networks

  Markus Hauru, Maarten Van Damme, Jutho Haegeman

  SciPost Phys. 10, 040 (2021) · published 17 February 2021 | Toggle abstract · pdf

  Several tensor networks are built of isometric tensors, i.e. tensors satisfying $W^\dagger W = \mathrm{I}$. Prominent examples include matrix product states (MPS) in canonical form, the multiscale entanglement renormalization ansatz (MERA), and quantum circuits in general, such as those needed in state preparation and quantum variational eigensolvers. We show how gradient-based optimization methods on Riemannian manifolds can be used to optimize tensor networks of isometries to represent e.g. ground states of 1D quantum Hamiltonians. We discuss the geometry of Grassmann and Stiefel manifolds, the Riemannian manifolds of isometric tensors, and review how state-of-the-art optimization methods like nonlinear conjugate gradient and quasi-Newton algorithms can be implemented in this context. We apply these methods in the context of infinite MPS and MERA, and show benchmark results in which they outperform the best previously-known optimization methods, which are tailor-made for those specific variational classes. We also provide open-source implementations of our algorithms.

  31 citations

• Bridging nano-optics and condensed matter formalisms in a unified description of inelastic scattering of relativistic electron beams

  Hugo Lourenço-Martins, Axel Lubk, Mathieu Kociak

  SciPost Phys. 10, 031 (2021) · published 10 February 2021 | Toggle abstract · pdf

  In the last decades, the blossoming of experimental breakthroughs in the domain of electron energy loss spectroscopy (EELS) has triggered a variety of theoretical developments. Those have to deal with completely different situations, from atomically resolved phonon mapping to electron circular dichroism passing by surface plasmon mapping. All of them rely on very different physical approximations and have not yet been reconciled, despite early attempts to do so. As an effort in that direction, we report on the development of a scalar relativistic quantum electrodynamic (QED) approach of the inelastic scattering of fast electrons. This theory can be adapted to describe all modern EELS experiments, and under the relevant approximations, can be reduced to any of the last EELS theories. In that aim, we present in this paper the state of the art and the basics of scalar relativistic QED relevant to the electron inelastic scattering. We then give a clear relation between the two once antagonist descriptions of the EELS, the retarded green Dyadic, usually applied to describe photonic excitations and the quasi-static mixed dynamic form factor (MDFF), more adapted to describe core electronic excitations of material. We then use this theory to establish two important EELS-related equations. The first one relates the spatially resolved EELS to the imaginary part of the photon propagator and the incoming and outgoing electron beam wavefunction, synthesizing the most common theories developed for analyzing spatially resolved EELS experiments. The second one shows that the evolution of the electron beam density matrix is proportional to the mutual coherence tensor, proving that quite universally, the electromagnetic correlations in the target are imprinted in the coherence properties of the probing electron beam.

  7 citations

• The multi-faceted inverted harmonic oscillator: Chaos and complexity

  Arpan Bhattacharyya, Wissam Chemissany, S. Shajidul Haque, Jeff Murugan, Bin Yan

  SciPost Phys. Core 4, 002 (2021) · published 8 February 2021 | Toggle abstract · pdf

  The harmonic oscillator is the paragon of physical models; conceptually and computationally simple, yet rich enough to teach us about physics on scales that span classical mechanics to quantum field theory. This multifaceted nature extends also to its inverted counterpart, in which the oscillator frequency is analytically continued to pure imaginary values. In this article we probe the inverted harmonic oscillator (IHO) with recently developed quantum chaos diagnostics such as the out-of-time-order correlator (OTOC) and the circuit complexity. In particular, we study the OTOC for the displacement operator of the IHO with and without a non-Gaussian cubic perturbation to explore genuine and quasi scrambling respectively. In addition, we compute the full quantum Lyapunov spectrum for the inverted oscillator, finding a paired structure among the Lyapunov exponents. We also use the Heisenberg group to compute the complexity for the time evolved displacement operator, which displays chaotic behaviour. Finally, we extended our analysis to N-inverted harmonic oscillators to study the behaviour of complexity at the different timescales encoded in dissipation, scrambling and asymptotic regimes.

  47 citations

• Shape of a sound wave in a weakly-perturbed Bose gas

  Oleksandr V. Marchukov, Artem G. Volosniev

  SciPost Phys. 10, 025 (2021) · published 3 February 2021 | Toggle abstract · pdf

  We employ the Gross-Pitaevskii equation to study acoustic emission generated in a uniform Bose gas by a static impurity. The impurity excites a sound-wave packet, which propagates through the gas. We calculate the shape of this wave packet in the limit of long wave lengths, and argue that it is possible to extract properties of the impurity by observing this shape. We illustrate here this possibility for a Bose gas with a trapped impurity atom -- an example of a relevant experimental setup. Presented results are general for all one-dimensional systems described by the nonlinear Schr\"odinger equation and can also be used in nonatomic systems, e.g., to analyze light propagation in nonlinear optical media. Finally, we calculate the shape of the sound-wave packet for a three-dimensional Bose gas assuming a spherically symmetric perturbation.

  3 citations