[Wien] lcore stops

[Fred Nastos]

Hi,

We are doing some calculations on different Si surfaces.
One of the surfaces we are looking at is hydrogen-terminated.
We are having trouble with this surface.

We took the positions from a pseudo-potential MD code.
The Si-H bond is roughly 2.8 bohr.

We first set the RKM to 3 and, and the H muffin tin
to 0.6 bohr, and all the Si tin radii to 2.2.  The
run_lapw script ran ok.

We increased RKM to 4, but this increased gmin too
much (as expected), so we tried increasing the H
muffin tin radius to 0.8.  We lowered the nearest
Si radius to 1.97.  This lowered gmin and the wien
code seemed to work fine.

We then increased RKM to 5.  This gives a gmin of
12.  When we execute run_lapw, it runs one complete
cycle.  In the second cycle, when it reaches lcore,
it crashes.  The error given is:

################### begin cut and paste ############
Invalid decimal character M was detected (unit=20).
(the relative position causing an error in a record = 12 )
MBER OF SYMMETRY OPERATIONS
?
  Error occurs at or near line 95 of insld_
  Called from or near line 129 of MAIN__
0.030u 0.000s 0:00.02 150.0%    0+0k 0+0io 250pf+0w
################## end cut and paste ###############

In case.outputc the last line is NSTOP= 210 FOR THE 1S ORBITAL.

This last atom that it stops at is the Si atom that
binds to H atom, so it has a smaller RMT (of 1.97) than
the other Si atoms (2.2 bohr).  Output was given for
the H and the other Si atoms.

This according the usersguide, in the lstart section,
seems to mean that lcore could not find a solution for
the 1s state.  Is that right, or is something else
wrong with lcore? (besides the fact that it is written
with too many goto statements!)

Is there a way around this?

Why would the calculation run with RKM=4 though?

Thank you for your time