[gpaw-users] NEGF-STM
Lars Smykalla
lars.smykalla at physik.tu-chemnitz.de
Wed Dec 28 13:52:34 CET 2011
Hi all,
I really can't get the stm module to work, not even the tutorial with H on
a Al slab. I hope someone can help me because I have no idea what to try
further.
When I do stm.scan() and stm.plot(), I get this UnicodeDecodeError:
Exception in Tkinter callback
Traceback (most recent call last):
File "/usr/lib/python2.7/lib-tk/Tkinter.py", line 1413, in __call__
return self.func(*args)
File
"/usr/lib/pymodules/python2.7/matplotlib/backends/backend_tkagg.py", line
245, in resize
self.show()
File
"/usr/lib/pymodules/python2.7/matplotlib/backends/backend_tkagg.py", line
248, in draw
FigureCanvasAgg.draw(self)
File "/usr/lib/pymodules/python2.7/matplotlib/backends/backend_agg.py",
line 394, in draw
self.figure.draw(self.renderer)
File "/usr/lib/pymodules/python2.7/matplotlib/artist.py", line 55, in
draw_wrapper
draw(artist, renderer, *args, **kwargs)
File "/usr/lib/pymodules/python2.7/matplotlib/figure.py", line 798, in
draw
func(*args)
File "/usr/lib/pymodules/python2.7/matplotlib/artist.py", line 55, in
draw_wrapper
draw(artist, renderer, *args, **kwargs)
File "/usr/lib/pymodules/python2.7/matplotlib/axes.py", line 1946, in
draw
a.draw(renderer)
File "/usr/lib/pymodules/python2.7/matplotlib/artist.py", line 55, in
draw_wrapper
draw(artist, renderer, *args, **kwargs)
File "/usr/lib/pymodules/python2.7/matplotlib/axis.py", line 1032, in
draw
self.label.draw(renderer)
File "/usr/lib/pymodules/python2.7/matplotlib/artist.py", line 55, in
draw_wrapper
draw(artist, renderer, *args, **kwargs)
File "/usr/lib/pymodules/python2.7/matplotlib/text.py", line 526, in draw
bbox, info = self._get_layout(renderer)
File "/usr/lib/pymodules/python2.7/matplotlib/text.py", line 305, in
_get_layout
clean_line, ismath = self.is_math_text(line)
File "/usr/lib/pymodules/python2.7/matplotlib/text.py", line 989, in
is_math_text
if cbook.is_math_text(s):
File "/usr/lib/pymodules/python2.7/matplotlib/cbook.py", line 1793, in
is_math_text
s = unicode(s)
File "/usr/lib/python2.7/encodings/utf_8.py", line 16, in decode
return codecs.utf_8_decode(input, errors, True)
UnicodeDecodeError: 'utf8' codec can't decode byte 0xc5 in position 0:
unexpected end of data
GPAW CLEANUP for serial binary: <type 'exceptions.UnicodeDecodeError'>
occured. Calling sys.exit()
I also tried to calculate a STM image of Si(111)-(7x7), but because of the
hexagonal cell, I get:
File "/usr/local/lib/python2.7/dist-packages/gpaw/transport/jstm.py", line
1667, in smallestbox
new_ct = np.dot(R(theta), ct.T).T
NameError: global name 'R' is not defined
GPAW CLEANUP for serial binary: <type 'exceptions.NameError'> occured.
Calling sys.exit()
So I suppose, in the code of jstm.py it should call rotate:
new_ct = np.dot(rotate(theta), ct.T).T
But after I changed this, I get a new error:
File "/usr/local/lib/python2.7/dist-packages/gpaw/transport/jstm.py", line
1667, in smallestbox
new_ct = np.dot(rotate(theta), ct.T).T
ValueError: matrices are not aligned
Greetings,
Lars
Am 02.12.2011, 18:58 Uhr, schrieb Lars Smykalla
<lars.smykalla at physik.tu-chemnitz.de>:
> Hi,
>
> I read a localized basis is necessary for STM simulation, so from the
> tutorial and documentation I it seems only possible for lcao mode.
> But now I read it also works when I create Wannier functions from a plane
> wave or grid based calculation.
>
> What python code do I have to use exactly to do this and to dump in this
> case the overlap matrix and the Hamiltonian matrix?
>
>
>
> Further, I don't understand in the NEGF-STM tutorial:
>
> # Here the keyword 'cvl' refers to the number of basis functions in the
> # in the convergence layer, i.e. for the present system four atomic
> layers
> are used.
> dump_hs(calc, 'srf', region='surface', cvl=4*4*9)
>
>
> There a 16 Ag atoms in 4 layers, but where does this 9 come from?
>
>
> Greetings,
> Lars
> _______________________________________________
> gpaw-users mailing list
> gpaw-users at listserv.fysik.dtu.dk
> https://listserv.fysik.dtu.dk/mailman/listinfo/gpaw-users
More information about the gpaw-users
mailing list