Comparative study of state-of-the-art matrix-product-state methods for lattice models with large local Hilbert spaces without U(1) symmetry

Lattice models consisting of high-dimensional local degrees of freedom without global particle-number conservation constitute an important problem class in the field of strongly correlated quantum many-body systems. For instance, they are realized in electron-phonon models, cavities, atom-molecule resonance models, or superconductors. In general, these systems elude a complete analytical treatment and need to be studied using numerical methods where matrix-product states (MPSs) provide a flexible and generic ansatz class. Typically, MPS algorithms scale at least quadratic in the dimension of the local Hilbert spaces. Hence, tailored methods, which truncate this dimension, are required to allow for efficient simulations. Here, we describe and compare three state-of-the-art MPS methods each of which exploits a different approach to tackle the computational complexity. We analyze the properties of these methods for the example of the Holstein model, performing high-precision calculations as well as a finite-size-scaling analysis of relevant ground-state observables. The calculations are performed at different points in the phase diagram yielding a comprehensive picture of the different approaches.