gkc1000 / potato Goto Github PK

Hi, I use the code calculate Single-site DMFT with 9 bath sites system in RCCSD(I hope to reproduce the results in the paper of PHYSICAL REVIEW B 100, 115154 (2019)). And i use:
nx=50
U=2
nw=9
mydmft, w, delta = dmft.hub_1d (nx, U, nw)
freqs = _get_linear_freqs(-5., 5., 128)[0]
ldos = mydmft.get_ldos(freqs, delta)
But the ldos can't have a good match with the paper results in the paper of FIG. 1.

Where did I go wrong？