I spent some days improving the phase diagram of spin-1 cases but there has been no significant change.

Here is a diagram just as usual, where we calculate the correlation length (logarithm) for each phase point at two different bond dimensions 64 and 128.
The main part is still composed of gapped points that have no change in the correlation length, which is fine; but for certain points gathered in the center of two diagonal directions, they seemed gapless but also sometimes calculated inaccurately, which seems quite creepy.
Then, I look closer into the line Jx=Jy, do the DMRG in a series (iteratively using the previous psi as the input of the next DMRG engine), and here is an interesting figure for the correlation length (the original value, not log)

The two peak points correspond to Jx=Jy=.84 and Jx=Jy=1.40 respectively.
I spent some days improving the phase diagram of spin-1 cases but there has been no significant change.

Here is a diagram just as usual, where we calculate the correlation length (logarithm) for each phase point at two different bond dimensions 64 and 128.
The main part is still composed of gapped points that have no change in the correlation length, which is fine; but for certain points gathered in the center of two diagonal directions, they seemed gapless but also sometimes calculated inaccurately, which seems quite creepy.
Then, I look closer into the line

Jx=Jy, do the DMRG in a series (iteratively using the previouspsias the input of the next DMRG engine), and here is an interesting figure for the correlation length (the original value, not log)The two peak points correspond to
Jx=Jy=.84andJx=Jy=1.40respectively.