Hello Kwant community,
I write since I am trying to simulate a Kitaev chain with disorder,
here distributed uniformly in the range [1,-1],
and I am experiencing troubles to implement the disorder itself.
In the absence of disorder,
the function that defines the tight-binding (open) Hamiltonian can be written
e.g. as follows:
def kitaev(L ,t, mu, Delta) :
syst = kwant.Builder()
lat = kwant.lattice.chain(norbs=2)
syst[lat(0)] = - (mu/2)*pauli3
for i in range(1,L):
syst[lat(i)] = -(mu/2)*pauli3
syst[lat(i), lat(i - 1)] = -(t/2)*pauli3 + 1j*(Delta/2)*pauli2
return syst
where the Pauli matrices have been previously defined, t is the hopping,
mu is the chemical potential, and Delta is the gap.
Now, I want to add disorder, proportional to W and uniformly
distributed in [-1,1]. The tutorial online for Kwant states that the
function kwant.digest.uniform(input, salt='') creates a random offset
in the interval [0,1]. Therefore I thought that the Hamiltonian could
be changed e.g. as follows:
def kitaev(L ,t, mu, Delta, W) :
syst = kwant.Builder()
lat = kwant.lattice.chain(norbs=2)
syst[lat(0)] = -(1/2)*(mu + W*(2*onsite - 1))*pauli3
for i in range(1,L):
syst[lat(i)] = -(1/2)*(mu + W*(2*onsite - 1))*pauli3
syst[lat(i), lat(i - 1)] = -(t/2)*pauli3 + 1j*(Delta/2)*pauli2
return syst
after that the onsite function has been defined as
def onsite(site, phi, salt):
return uniform(repr(site), salt)
(I am still trying to follow the tutorial).
However, when trying to calculate eigenvalues and eigenvectors,
via the code
L = 50
W = 0
mub = 3
systb = kitaev(L ,1, mub, 1, W)
systb = systb.finalized()
hamatb = systb.hamiltonian_submatrix(sparse = False)
evalsb, evecsb = la.eigh(hamatb)
I encountered the problem
-----------------------------------------------------------------------------------------
TypeError Traceback (most recent call last)
<ipython-input-10-12a574eefcd0> in <module>
3
4 mub = 3
----> 5 systb = kitaev(L ,1, mub, 1, W)
6 systb = systb.finalized()
7 hamatb = systb.hamiltonian_submatrix(sparse = False)
<ipython-input-9-661d5ddb539a> in kitaev(L, t, mu, Delta, W)
4 lat = kwant.lattice.chain(norbs=2)
5
----> 6 syst[lat(0)] = -(1/2)*(mu + W*(2*onsite - 1))*pauli3
7
8 for i in range(1,L):
TypeError: unsupported operand type(s) for *: 'int' and 'function'
—————————————————————————————------
that instead does no occur in the absence of disorder (then when the
first definition above for the Kitaev chain is assumed).
Probably I did not understand properly how to use the
function uniform(repr(site), salt).
Can you help me ?
Thank you very much in advance and best regards
L. L.
