<html>
  <head>

    <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
  </head>
  <body text="#000000" bgcolor="#FFFFFF">
    <div class="moz-text-flowed" style="font-family: -moz-fixed;
      font-size: 14px;" lang="x-western">Hello everyone,
      <br>
      <br>
      I'm trying to generate a series of (110) surface slabs of varying
      thickness for a zincblende crystal. I'm doing this using the
      structeditor functions in an octave script (and have done this
      with both wien2k (11), 12 and 13).
      <br>
      <br>
      At certain thicknesses of the slab a problem occurs, which I've
      seen in the past for other systems, where it appears to make a
      diagonal cut at the bottom of the slab along a {001} direction.
      Does anyone by any chance have a patch to resolve this? I know
      ways around it myself but being able to use the provided functions
      directly would be nice.
      <br>
      <br>
      I've added the bulk struct and octave script needed to reproduce
      the problem below.
      <br>
      <br>
      <br>
      Thanks,
      <br>
      <br>
      Michael Sluydts
      <br>
      <br>
      --------------------bulk.struct------------------------
      <br>
      blebleble
      <br>
      F&nbsp;&nbsp; LATTICE,NONEQUIV.ATOMS&nbsp;&nbsp; 2&nbsp; 216 F-43m
      <br>
      MODE OF CALC=RELA unit=bohr
      <br>
      &nbsp;10.847217 10.847217 10.847217 90.000000 90.000000 90.000000
      <br>
      ATOM&nbsp; -1: X=0.00000000 Y=0.00000000 Z=0.00000000
      <br>
      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT=15
      <br>
      Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp; R0=.000050000 RMT=&nbsp;&nbsp; 2.00000&nbsp;&nbsp; Z:&nbsp; 30.00000
      <br>
      LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 1.0000000 0.0000000 0.0000000
      <br>
      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 1.0000000 0.0000000
      <br>
      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 0.0000000 1.0000000
      <br>
      ATOM&nbsp; -2: X=0.25000000 Y=0.25000000 Z=0.25000000
      <br>
      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT=15
      <br>
      Se&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp; R0=.000050000 RMT=&nbsp;&nbsp; 2.00000&nbsp;&nbsp; Z:&nbsp; 34.00000
      <br>
      LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 1.0000000 0.0000000 0.0000000
      <br>
      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 1.0000000 0.0000000
      <br>
      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 0.0000000 1.0000000
      <br>
      &nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NUMBER OF SYMMETRY OPERATIONS
      <br>
      <br>
      ------------------------script.m----------------------------
      <br>
      bulk = makeconventional(loadstruct('bulk.struct'));
      <br>
      surf = makesurface(bulk,[1 1
      0],2,(bulk.a(3)/2*sin(pi/4))*10+0.01,37.8);
      <br>
      savestruct(makeprimitive(surf),['bug.struct'],0);
      <br>
    </div>
  </body>
</html>