Please use this identifier to cite or link to this item: https://doi.org/10.1103/PhysRevE.82.036702
Title: Dynamical simulations of classical stochastic systems using matrix product states
Authors: Johnson, T.H.
Clark, S.R. 
Jaksch, D. 
Issue Date: 16-Sep-2010
Citation: Johnson, T.H., Clark, S.R., Jaksch, D. (2010-09-16). Dynamical simulations of classical stochastic systems using matrix product states. Physical Review E - Statistical, Nonlinear, and Soft Matter Physics 82 (3) : -. ScholarBank@NUS Repository. https://doi.org/10.1103/PhysRevE.82.036702
Abstract: We adapt the time-evolving block decimation (TEBD) algorithm, originally devised to simulate the dynamics of one-dimensional quantum systems, to simulate the time evolution of nonequilibrium stochastic systems. We describe this method in detail; a system's probability distribution is represented by a matrix product state (MPS) of finite dimension and then its time evolution is efficiently simulated by repeatedly updating and approximately refactorizing this representation. We examine the use of MPS as an approximation method, looking at parallels between the interpretations of applying it to quantum state vectors and probability distributions. In the context of stochastic systems we consider two types of factorization for use in the TEBD algorithm: non-negative matrix factorization (NMF), which ensures that the approximate probability distribution is manifestly non-negative, and the singular value decomposition (SVD). Comparing these factorizations, we find the accuracy of the SVD to be substantially greater than current NMF algorithms. We then apply TEBD to simulate the totally asymmetric simple exclusion process (TASEP) for systems of up to hundreds of lattice sites in size. Using exact analytic results for the TASEP steady state, we find that TEBD reproduces this state such that the error in calculating expectation values can be made negligible even when severely compressing the description of the system by restricting the dimension of the MPS to be very small. Out of the steady state we show for specific observables that expectation values converge as the dimension of the MPS is increased to a moderate size. © 2010 The American Physical Society.
Source Title: Physical Review E - Statistical, Nonlinear, and Soft Matter Physics
URI: http://scholarbank.nus.edu.sg/handle/10635/115070
ISSN: 15393755
DOI: 10.1103/PhysRevE.82.036702
Appears in Collections:Staff Publications

Show full item record
Files in This Item:
There are no files associated with this item.

Google ScholarTM

Check

Altmetric


Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.