Thanks again:)

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.


Chong Wang

From: Joseph Weston <>
Sent: Wednesday, May 11, 2016 8:44 PM
To: Chong Wang
Subject: Re: [Kwant] scatter without a center area

Hello again,

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!

Ok, I think we've got mixed in our nomenclature somewhere. My basic question
still stands, though. Have you tried doing as Anton and I suggested and just
putting a single "slice" of one of the leads as a central region, and seeing how
long it takes to compute whatever you are interested in? This solution is the easiest,
and my guess is that the cost will not be inhibitive.