Search

We link linear prediction of Chebyshev and Fourier expansions to analytic continuation. We push the resolution in the Chebyshev-based computation of T = 0 many-body spectral functions to a much higher precision by deriving a modified Chebyshev series expansion that allows to reduce the expansion order by a factor ~1/6. We show that in a certain limit the Chebyshev technique becomes equivalent to computing spectral functions via time evolution and subsequent Fourier transform. This introduces a novel recursive time-evolution algorithm that instead of the group operator e-iHt only involves the action of the generator H. For quantum impurity problems, we introduce an adapted discretization scheme for the bath spectral function. We discuss the relevance of these results for matrix product state (MPS) based DMRG-type algorithms, and their use within the dynamical mean-field theory (DMFT). We present strong evidence that the Chebyshev recursion extracts less spectral information from H than time evolution algorithms when fixing a given amount of created entanglement. [ABSTRACT FROM AUTHOR]

We solve nonequilibrium dynamical mean-field theory (DMFT) using matrix product states (MPS). This allows us to treat much larger bath sizes and by that reach substantially longer times (factor ∼2–3) than with exact diagonalization. We show that the star geometry of the underlying impurity problem can have substantially better entanglement properties than the previously favored chain geometry. This has immense consequences for the efficiency of an MPS-based description of general impurity problems: in the case of equilibrium DMFT, it leads to an orders-of-magnitude speedup. We introduce an approximation for the two-time hybridization function that uses time-translational invariance, which can be observed after a certain relaxation time after a quench to a time-independent Hamiltonian. [ABSTRACT FROM AUTHOR]

We compute the spectral functions for the two-site dynamical cluster theory and for the two-orbital dynamical mean-field theory in the density-matrix renormalization group (DMRG) framework using Chebyshev expansions represented with matrix product states (MPS). We obtain quantitatively precise results at modest computational effort through technical improvements regarding the truncation scheme and the Chebyshev rescaling procedure. We furthermore establish the relation of the Chebyshev iteration to real-time evolution and discuss technical aspects as computation time and implementation in detail. [ABSTRACT FROM AUTHOR]

A quantum many-body system that is prepared in the ground state of an integrable Hamiltonian does not directly thermalize after a sudden small parameter quench away from integrability. Rather, it will be trapped in a prethermalized state and can thermalize only at a later stage. We discuss several examples for which this prethermalized state shares some properties with the nonthermal steady state that emerges in the corresponding integrable system. These examples support the notion that nonthermal steady states in integrable systems may be viewed as prethermalized states that never decay further. Furthermore, we show that prethermalization plateaus are under certain conditions correctly predicted by generalized Gibbs ensembles, which are the appropriate extension of standard statistical mechanics in the presence of many constants of motion. This establishes that the relaxation behaviors of integrable and nearly integrable systems are continuously connected and described by the same statistical theory. [ABSTRACT FROM AUTHOR]

Many numerical techniques for the description of quantum systems that are coupled to a continuous bath require the discretization of the latter. To this end, a wealth of methods has been developed in the literature, which we classify as (i) direct discretization, (ii) orthogonal polynomial, and (iii) numerical optimization strategies. We recapitulate strategies (i) and (ii) to clarify their relation. For quadratic Hamiltonians, we show that (ii) is the best strategy in the sense that it gives the numerically exact time evolution up to a maximum time tmax, for which we give a simple expression. For nonquadratic Hamiltonians, we show that no such best strategy exists. We present numerical examples relevant to open quantum systems and strongly correlated systems, as treated by dynamical mean-field theory (DMFT). [ABSTRACT FROM AUTHOR]