Finite Density Matrix Renormalisation Group Algorithm for Anyonic Systems

The numerical study of anyonic systems is known to be highly challenging due to their non-bosonic, non-fermionic particle exchange statistics, and with the exception of certain models for which analytical solutions exist, very little is known about their collective behaviour as a result. Meanwhile, the density matrix renormalisation group (DMRG) algorithm is an exceptionally powerful numerical technique for calculating the ground state of a low-dimensional lattice Hamiltonian, and has been applied to the study of bosonic, fermionic, and group-symmetric systems. The recent development of a tensor network formulation for anyonic systems opened up the possibility of studying these systems using algorithms such as DMRG, though this has proved challenging both in terms of programming complexity and computational cost. This paper presents the implementation of DMRG for finite anyonic systems, including a detailed scheme for the implementation of anyonic tensors with optimal scaling of computational cost. The anyonic DMRG algorithm is demonstrated by calculating the ground state energy of the Golden Chain, which has become the benchmark system for the numerical study of anyons, and is shown to produce results comparable to those of the anyonic TEBD algorithm and superior to the variationally optimised anyonic MERA, at far lesser computational cost.Comment: 24 pages, 37 figure files (25 floating figures). RevTeX 4.1. Minor changes for clarity in Figs. 9 & 11, matching published versio

Tensor network states and algorithms in the presence of a global U(1) symmetry

Tensor network decompositions offer an efficient description of certain many-body states of a lattice system and are the basis of a wealth of numerical simulation algorithms. In a recent paper [arXiv:0907.2994v1] we discussed how to incorporate a global internal symmetry, given by a compact, completely reducible group G, into tensor network decompositions and algorithms. Here we specialize to the case of Abelian groups and, for concreteness, to a U(1) symmetry, often associated with particle number conservation. We consider tensor networks made of tensors that are invariant (or covariant) under the symmetry, and explain how to decompose and manipulate such tensors in order to exploit their symmetry. In numerical calculations, the use of U(1) symmetric tensors allows selection of a specific number of particles, ensures the exact preservation of particle number, and significantly reduces computational costs. We illustrate all these points in the context of the multi-scale entanglement renormalization ansatz.Comment: 22 pages, 25 figures, RevTeX

Further Series Studies of the Spin-1/2 Heisenberg Antiferromagnet at T=0: Magnon Dispersion and Structure Factors

We have extended our previous series studies of quantum antiferromagnets at zero temperature by computing the one-magnon dispersion curves and various structure factors for the linear chain, square and simple cubic lattices. Many of these results are new; others are a substantial extension of previous work. These results are directly comparable with neutron scattering experiments and we make such comparisons where possible.Comment: 15 pages, 12 figures, revised versio

Controlling Mixing Inside a Droplet by Time Dependent Rigid-body Rotation

The use of microscopic discrete fluid volumes (i.e., droplets) as microreactors for digital microfluidic applications often requires mixing enhancement and control within droplets. In this work, we consider a translating spherical liquid droplet to which we impose a time periodic rigid-body rotation which we model using the superposition of a Hill vortex and an unsteady rigid body rotation. This perturbation in the form of a rotation not only creates a three-dimensional chaotic mixing region, which operates through the stretching and folding of material lines, but also offers the possibility of controlling both the size and the location of the mixing. Such a control is achieved by judiciously adjusting the three parameters that characterize the rotation, i.e., the rotation amplitude, frequency and orientation of the rotation. As the size of the mixing region is increased, complete mixing within the drop is obtained.Comment: 6 pages, 6 figure

Matrix product states for anyonic systems and efficient simulation of dynamics

Matrix product states (MPS) have proven to be a very successful tool to study lattice systems with local degrees of freedom such as spins or bosons. Topologically ordered systems can support anyonic particles which are labeled by conserved topological charges and collectively carry non-local degrees of freedom. In this paper we extend the formalism of MPS to lattice systems of anyons. The anyonic MPS is constructed from tensors that explicitly conserve topological charge. We describe how to adapt the time-evolving block decimation (TEBD) algorithm to the anyonic MPS in order to simulate dynamics under a local and charge-conserving Hamiltonian. To demonstrate the effectiveness of anyonic TEBD algorithm, we used it to simulate (i) the ground state (using imaginary time evolution) of an infinite 1D critical system of (a) Ising anyons and (b) Fibonacci anyons both of which are well studied, and (ii) the real time dynamics of an anyonic Hubbard-like model of a single Ising anyon hopping on a ladder geometry with an anyonic flux threading each island of the ladder. Our results pertaining to (ii) give insight into the transport properties of anyons. The anyonic MPS formalism can be readily adapted to study systems with conserved symmetry charges, as this is equivalent to a specialization of the more general anyonic case.Comment: 18 pages, 15 figue

Simulation of braiding anyons using Matrix Product States

Anyons exist as point like particles in two dimensions and carry braid statistics which enable interactions that are independent of the distance between the particles. Except for a relatively few number of models which are analytically tractable, much of the physics of anyons remain still unexplored. In this paper, we show how U(1)-symmetry can be combined with the previously proposed anyonic Matrix Product States to simulate ground states and dynamics of anyonic systems on a lattice at any rational particle number density. We provide proof of principle by studying itinerant anyons on a one dimensional chain where no natural notion of braiding arises and also on a two-leg ladder where the anyons hop between sites and possibly braid. We compare the result of the ground state energies of Fibonacci anyons against hardcore bosons and spinless fermions. In addition, we report the entanglement entropies of the ground states of interacting Fibonacci anyons on a fully filled two-leg ladder at different interaction strength, identifying gapped or gapless points in the parameter space. As an outlook, our approach can also prove useful in studying the time dynamics of a finite number of nonabelian anyons on a finite two-dimensional lattice.Comment: Revised version: 20 pages, 14 captioned figures, 2 new tables. We have moved a significant amount of material concerning symmetric tensors for anyons --- which can be found in prior works --- to Appendices in order to streamline our exposition of the modified Anyonic-U(1) ansat

Quantum mechanics without spacetime II : noncommutative geometry and the free point particle

In a recent paper we have suggested that a formulation of quantum mechanics should exist, which does not require the concept of time, and that the appropriate mathematical language for such a formulation is noncommutative differential geometry. In the present paper we discuss this formulation for the free point particle, by introducing a commutation relation for a set of noncommuting coordinates. The sought for background independent quantum mechanics is derived from this commutation relation for the coordinates. We propose that the basic equations are invariant under automorphisms which map one set of coordinates to another- this is a natural generalization of diffeomorphism invariance when one makes a transition to noncommutative geometry. The background independent description becomes equivalent to standard quantum mechanics if a spacetime manifold exists, because of the proposed automorphism invariance. The suggested basic equations also give a quantum gravitational description of the free particle.Comment: 8 page