iDMRG shows poor convergence in the 1/3 FCI phase of Fermionic Haldane Model

[XiaohanYang]

Hi everyone,
I am trying to reproduce the result, especially the charge pumping result of arXiv:1407.6985. I directly use the module tenpy.models.haldane.FermionicHaldaneModel since it's just the Hamiltonian eq(1) in arXiv:1407.6985. To be concrete, let me show you the Hamiltonian:
\(H = -t_1\sum_{<ij>}c_i^\dagger c_j-t_2\sum_{<<ij>>}e^{i\phi_{ij}}c_i^\dagger c_j+m\sum_u (n_{u,A}-n_{u, B})+V\sum_{u}n_{u, A}n_{u, B}\)
The geometry is an infinite cylinder with circumference \(L_y\) which is the number of unit cells in the y direction(the finite periodic direction). The unit cell of the uMPS includes \(L_y\) lattice unit cells in a leg(\(2L_y\) sites). I have chosen two parameter dictionaries to reproduce Fig.2. of arXiv:1407.6985. They are:

Code: Select all

# Parameter for \nu=1 Chern Insulator
param_CI = {   "bc_MPS": "infinite",
              	"bc_x": "periodic",
            	"bc_y": "cylinder", 
            	"conserve": "N",
            	"Lx": 1,
           	"Ly": 3,
            	"order": "default",
            	"mu": 0.,
            	"V": 0.0,
            	"t1": -1., # leave t2 to be the default value
           	"phi_ext": phi_ext
        	}
# Parameter for \nu=1/3 Fractional Chern Insulator
param_FCI = {   "bc_MPS": "infinite",
              	"bc_x": "periodic",
            	"bc_y": "cylinder", 
            	"conserve": "N",
            	"Lx": 1,
           	"Ly": 6,
            	"order": "default",
            	"mu": 0.,
            	"V": 1.0,
            	"t1": -1., # leave t2 to be the default value
           	"phi_ext": phi_ext
        	}
        

The filling numbers are 1 fermion per lattice unit cell and 1 fermion per 3 lattice cells. These parameters are the same as the two cases in arXiv:1407.6985. The iDMRG for the Chern insulator case converges well with bond dimension 200 and the expected result has been obtained. For the FCI case, however, the iDMRG shows poor convergence in the suggested bond dimension 500. The errors of both energy and entropy fluctuate and cannot converge.

Since in both cases(CI and FCI) the systems are gapped, I think there is no gaplessness problem. I have tried to increase the interaction \(V\) in the FCI case to avoid competing with the metal phase but it doesn't work.

So, I want some suggestions from people who are familiar with iDMRG implementation or FCI systems.

Sincerely,
Xiaohan Yang

[XiaohanYang]

I give the Two site iDMRG such parameter dictionary:

Code: Select all

dmrg_params = {'mixer': True,
               'mixer_params': {
                    'amplitude': 5.e-5,
                    'decay': 1.1,
                    'disable_after': 100
                }, 
                'trunc_params': {
                    'svd_min': 1.e-10,
                },
                'max_E_err': 1.e-10,
                'max_S_err': 1.e-8,
                'min_sweeps': 10,
                'max_sweeps': 2000,
                'lanczos_params':{
                    'reortho': False, 
                    'N_min':5, 
                    'N_max':500
                },
                'chi_max': 500
            }

[XiaohanYang]

BTW, I found some insight from DHIMAN002's post:viewtopic.php?p=281&hilit=haldane#p281. In his post, Johannes replied that it may be the case that we need a large uMPS unit cell including several legs of the lattice. I have run the finite DMRG with \(L_y=3\) and \(L_x=30\) for both CI and FCI cases. The Fermion density profiles are:
CI:

CI.png (170.05 KiB) Viewed 226 times

FCI:

FCI.png (214.97 KiB) Viewed 226 times

You can see that CDW-like order exists in the case of FCI, so I wonder whether I need to use a large MPS unit cell in the infinite case. If the bond dimension is increased, will this CDW order be suppressed? When the system is in the FCI phase, does this CDW order always exist?