Hi Luca,
> thanks again for your useful reply. > > If I understand correctly, then I am not sure > how to simulate the physics of disorder. > > Indeed, that task requires to create various > (in principle infinite in number) configurations, > measure some physical quantity on each of them, and then > to take the average of all these measurements on the created configurations. > > How can I create the different configurations ? > Starting e.g. from the piece of code Indeed, typically you'll want to take the average over different disorder configurations. This means using different values for the 'salt' parameter. The actual values used for 'salt' are not important, as long as they differ; 'kwant.digest.uniform' is a hash function (see https://en.wikipedia.org/wiki/Salt_(cryptography) ). You would do something like the following: syst = kitaev(L=...).finalized() params = dict(mu=..., W=..., t=..., Delta=...) # We (arbitrarily) choose the salts as the strings '0', '1', ... salts = map(str, range(100)) for params['salt'] in salts: ham = syst.hamiltonian_submatrix(params=params, sparse=False) evals = np.linalg.eigvalsh(ham) ... Happy Kwanting, Joe
signature.asc
Description: OpenPGP digital signature