Convergence Problem for Quantum Impurity Systems
I am currently trying to calculate the entanglement entropy of a one-dimensional Spin chain with an impurity dot in the middle, i.e., \( J_{xx} = [J,...,J, \lambda J, \lambda J, J, ..., J].\). For testing purposes, I firstly chose \( J_z = 0 \) to compare the results with exact diagonalization methods (non-interacting tight-binding chain via Jordan-Wigner transformation). There appear to be large convergence issues especially for small \( \lambda \) and moderate/large \(N\). I have tried mixing, choosing different initial guesses, and played with the cutoff parameters. If I choose only one impurity bond at the edges, i.e., \( J_{xx} = [\lambda J,..., J]\), everything seems to be okay and I can reproduce results from the literature. My guess is that the initial state should be chosen more wisely than a simple Neel state or related states to account for long range entanglement.
Concretely, I am thinking about getting the ground states for a system with an impurity on the right (left), and combine both ground states to get an initial guess for the impurity in the middle. However, I am not sure how to properly implement the product of two ground states as an initial guess.
Any suggestions/experiences with impurity problems in DMRG would be greatly appreciated!