Matrix product state algorithms for Gaussian fermionic states

While general quantum many-body systems require exponential resources to be simulated on a classical computer, systems of non-interacting fermions can be simulated exactly using polynomially scaling resources. Such systems may be of interest in their own right, but also occur as effective models in numerical methods for interacting systems, such as Hartree-Fock, density functional theory, and many others. Often it is desirable to solve systems of many thousand constituent particles, rendering these simulations computationally costly despite their polynomial scaling. We demonstrate how this scaling can be improved by adapting methods based on matrix product states, which have been enormously successful for low-dimensional interacting quantum systems, to the case of free fermions. Compared to the case of interacting systems, our methods achieve an exponential speedup in the entanglement entropy of the state. We demonstrate their use to solve systems of up to one million sites with an effective MPS bond dimension of 10^15.

[1]  P. Wiegmann,et al.  Thermodynamics of the s-d exchange model (Kondo problem) , 1981 .

[2]  Michael Wimmer,et al.  Algorithm 923: Efficient Numerical Computation of the Pfaffian for Dense and Banded Skew-Symmetric Matrices , 2011, TOMS.

[3]  Benni Reznik,et al.  BCS-like modewise entanglement of fermion Gaussian states , 2004 .

[4]  ScienceDirect Superlattices and microstructures , 1985 .

[5]  I. Peschel LETTER TO THE EDITOR: Calculation of reduced density matrices from correlation functions , 2002, cond-mat/0212631.

[6]  Roman Orus,et al.  A Practical Introduction to Tensor Networks: Matrix Product States and Projected Entangled Pair States , 2013, 1306.2164.

[7]  Román Orús,et al.  Advances on tensor network theory: symmetries, fermions, entanglement, and holography , 2014, 1407.6552.

[8]  Scott Kirkpatrick,et al.  Conductivity of the disordered linear chain , 1981 .

[9]  J. Cirac,et al.  Generalized Hartree–Fock theory for interacting fermions in lattices: numerical methods , 2010, 1005.5284.

[10]  Umesh Vazirani,et al.  Rigorous RG Algorithms and Area Laws for Low Energy Eigenstates in 1D , 2016, Communications in Mathematical Physics.

[11]  U. Schollwoeck The density-matrix renormalization group in the age of matrix product states , 2010, 1008.3477.

[12]  Eduardo R. Mucciolo,et al.  The recursive Green’s function method for graphene , 2013, Journal of Computational Electronics.

[13]  J. Eisert,et al.  Holography and criticality in matchgate tensor networks , 2017, Science Advances.

[14]  S. Datta Nanoscale device modeling: the Green’s function method , 2000 .

[15]  Christopher T. Chubb,et al.  Hand-waving and interpretive dance: an introductory course on tensor networks , 2016, 1603.03039.

[16]  U. Vazirani,et al.  A polynomial time algorithm for the ground state of one-dimensional gapped local Hamiltonians , 2015, Nature Physics.

[17]  Xavier Waintal,et al.  Numerical simulations of time-resolved quantum electronics , 2013, 1307.6419.

[18]  David Gosset,et al.  Complexity of Quantum Impurity Problems , 2016, 1609.00735.