I have (stupidly) tried to copy two slices in central region (one for each lead). After 3 hours, the memory blows up. My workstation has 128Gb DDR4 memory. (kwant is compiled against mumps 5.0.1 and multithreading is working)
I am now trying to use only one slice in central region and I will keep you informed whether this will work.
I am calculating a real system, with 32 orbitals per unit cell. Moreover, I need to use supercells so that I can simulate periodical boundary conditions. For example, I want to scatter along x direction. Also, I want z direction to be finite and y direction to be periodical. Since kwant does not allow hoppings beyond nearest unit cell in leads, I have to make a supercell that contains all unit cells in y direction!