Hello everyone, i would appreciate any suggestion you could give me to
solve this
Im trying to plot the conductance for a two dimensional array of 3 basis
atoms, i have defined the atoms positions, the onsite energies and hoppings
matrices and the leads but it is showing me this error:
(i search in the mailing list but theres no info about it)
ValueError Traceback (most recent call
last)<ipython-input-1-ee484ce5ea50> in <module>() 111 112 if
__name__ == '__main__':--> 113 main()
<ipython-input-1-ee484ce5ea50> in main() 107
site_color=family_colors) 108 syst = syst.finalized()--> 109
plot_conductance(syst, energies=[0.01 * i - 0.3 for i in range(100)])
110 111
<ipython-input-1-ee484ce5ea50> in plot_conductance(syst, energies)
85 data = [] 86 for energy in energies:---> 87
smatrix = kwant.smatrix(syst, energy) 88
data.append(smatrix.transmission(1, 0)) 89 pyplot.figure()
~/anaconda3/lib/python3.6/site-packages/kwant/solvers/common.py in
smatrix(self, sys, energy, args, out_leads, in_leads,
check_hermiticity, params) 370 linsys, lead_info =
self._make_linear_sys(syst, in_leads, energy, args, 371
check_hermiticity, False,--> 372
params=params)
373 374 kept_vars = np.concatenate([coords for i, coords
in
~/anaconda3/lib/python3.6/site-packages/kwant/solvers/common.py in
_make_linear_sys(self, sys, in_leads, energy, args, check_hermiticity,
realspace, params) 162 lhs, norb =
syst.hamiltonian_submatrix(args, sparse=True, 163
return_norb=True,--> 164
params=params)[:2] 165 lhs =
getattr(lhs, 'to' + self.lhsformat)() 166 lhs = lhs -
energy * sp.identity(lhs.shape[0], format=self.lhsformat)
kwant/_system.pyx in kwant._system.hamiltonian_submatrix()
kwant/_system.pyx in kwant._system.make_sparse_full()
ValueError: Hopping from site 407 to site 814 does not match the
dimensions of onsite Hamiltonians of these sites.
Here is the code so far:
import kwant
import tinyarray
import numpy
from matplotlib import pyplot
%matplotlib notebook
VecPrim = [(3.176, 0, 0), (1.588, 2.7504, 0), (0, 0, 17.49)]
base = [(0, 0, 0), (1.588, 0.9168, 0.8745), (0, 0, 1.749)]
lat = kwant.lattice.general(prim_vecs=VecPrim, basis=base)
c1, c2, c3 = lat.sublattices
c1_ons = numpy.array([[-1.888842, -0.014064-0.000409j, 0.026908+0.003422j],
[-0.014064+0.000409j, -1.784936, -0.031384+0.021565j],
[0.026908-0.003422j, -0.031384-0.021565j, 0.710443]])
c2_ons = numpy.array([[7.168134, -0.005189-0.047376j, 0.001352+0.000103j,
-0.001189-0.00095j, -0.048166-0.007663j, 0.429882+0.00052j],
[-0.005189+0.047376j, 0.245358, 0.113162-0.072937j,
0.660849-0.06348j, -0.017571+0.035587j, 0.020961-0.017679j],
[0.001352-0.000103j, 0.113162+0.072937j, 1.338059,
-0.003126-0.016285j, -0.624702-0.065652j, 0.002353+0.001639j],
[-0.001189+0.00095j, 0.660849+0.06348j,
-0.003126+0.016285j, 1.325337, 0.110694+0.01548j, 0.009041+0.004483j],
[-0.048166+0.007663j, -0.017571-0.035587j,
-0.624702+0.065652j, 0.110694-0.01548j, 0.362102, -0.030934+0.002401j],
[0.429882-0.00052j, 0.020961+0.017679j,
0.002353-0.001639j, 0.009041-0.004483j, -0.030934-0.002401j, -0.650839]])
c3_ons = numpy.array([[-3.398186, 0.001839+0.001631j, -0.117366+0.00498j],
[0.001839-0.001631j, -3.401815, 0.023339+0.009067j],
[-0.117366-0.00498j, 0.023339-0.009067j, -0.637207]])
c1_c3_hop = numpy.array([[-1.852931-0.147524j, -0.274338+0.015064j,
-0.103127-0.001438j],
[0.274586+0.034451j, -1.796449+0.055331j,
0.004816+0.00237j],
[-0.015051+0.001677j, -0.034396-0.019407j,
3.410475-0.001248j]])
c1_c2_hop = numpy.array([[0.003644-0.00088j, -0.046131-0.04229j,
0.23188-0.001924j],
[-0.305706+0.050514j, -1.963829+0.280977j,
-0.597505+0.078772j],
[0.638425+0.083597j, -1.171301+0.088146j,
-0.146976+0.012943j],
[-1.105645-0.086829j, -1.112788+0.056971j,
-0.210584-0.002248j],
[1.463462+0.012994j, 0.174769-0.106872j,
0.177022-0.031485j],
[-0.01338+0.00264j, -0.058159+0.007987j,
0.953363+0.001855j]])
c2_c3_hop = numpy.array([[0.184152+0.002778j, 0.321743-0.010571j,
0.183538-0.00067j],
[0.275145+0.007691j, -0.54144+0.012815j,
0.045847+0.015379j],
[0.47922+0.015775j, 1.256012+0.007326j,
-0.523174+0.00585j],
[1.255561-0.002255j, 1.903063+0.010647j,
-0.876392+0.010246j],
[1.282837-0.091022j, -0.54898+0.001161j,
0.057341-0.004095j],
[0.02765-0.001793j, -0.002902-0.001866j,
-0.746426-0.000748j]])
def make_cuboid(t=1.0, a=15, b=10, c=5):
def cuboid_shape(pos):
x, y, z = pos
return 0 <= x < a and 0 <= y < b and 0 <= z < c
syst = kwant.Builder()
syst[c1.shape(cuboid_shape, (0, 0, 0))] = c1_ons
syst[c2.shape(cuboid_shape, (0, 1, 0))] = c2_ons
syst[c3.shape(cuboid_shape, (0, 0, 1))] = c3_ons
syst[kwant.builder.HoppingKind((0, 0, 0), c2, c1)] = c1_c2_hop
syst[kwant.builder.HoppingKind((-1, 0, 0), c2, c1)] = c1_c2_hop
syst[kwant.builder.HoppingKind((0, -1, 0), c2, c1)] = c1_c2_hop
syst[kwant.builder.HoppingKind((0, 0, 0), c3, c2)] = c2_c3_hop
syst[kwant.builder.HoppingKind((1, 0, 0), c3, c2)] = c2_c3_hop
syst[kwant.builder.HoppingKind((0, 1, 0), c3, c2)] = c2_c3_hop
lead = kwant.Builder(kwant.TranslationalSymmetry((-3.176, 0, 0)))
lead[c1.shape(cuboid_shape, (0, 0, 0))] = c1_ons
lead[c2.shape(cuboid_shape, (0, 1, 0))] = c2_ons
lead[c3.shape(cuboid_shape, (0, 0, 1))] = c3_ons
lead[kwant.builder.HoppingKind((0, 0, 0), c2, c1)] = c1_c2_hop
lead[kwant.builder.HoppingKind((-1, 0, 0), c2, c1)] = c1_c2_hop
lead[kwant.builder.HoppingKind((0, -1, 0), c2, c1)] = c1_c2_hop
lead[kwant.builder.HoppingKind((0, 0, 0), c3, c2)] = c2_c3_hop
lead[kwant.builder.HoppingKind((1, 0, 0), c3, c2)] = c2_c3_hop
lead[kwant.builder.HoppingKind((0, 1, 0), c3, c2)] = c2_c3_hop
syst.attach_lead(lead)
syst.attach_lead(lead.reversed())
return syst
def plot_conductance(syst, energies):
data = []
for energy in energies:
smatrix = kwant.smatrix(syst, energy)
data.append(smatrix.transmission(1, 0))
pyplot.figure()
pyplot.plot(energies, data)
pyplot.xlabel("energy [t]")
pyplot.ylabel("conductance [e^2/h]")
pyplot.show()
def main():
syst = make_cuboid()
kwant.plot(syst)
syst = make_cuboid(a=100, b=28, c=4)
def family_colors(site):
if site.family == c1:
return 'yellow'
elif site.family == c2:
return 'gray'
else:
return 'blue'
kwant.plot(syst, site_size=0.05, site_lw=0.01, hop_lw=0.01,
site_color=family_colors)
syst = syst.finalized()
plot_conductance(syst, energies=[0.01 * i - 0.3 for i in range(100)])
if __name__ == '__main__':
main()