Dear All,<br><br>Is it possible to perform a calculation in Wien2k for the charge density of a very heavy atom eg. #116 or 120? Also, would the answers be sensible?<br><br>Eg. for #116 I tried this case.inst:<br><br>Up<br>
Rn 5  <br>5, 3,3.0  N<br>5, 3,3.0  N<br>5,-4,4.0  N<br>5,-4,4.0  N<br>6, 2,1.0  N<br>6, 2,0.0  N<br>7,-1,1.0  N<br>7,-1,1.0  N<br>7, 1,1.0  N<br>7, 1,1.0  N<br>****<br>****         END of input (instgen_lapw)<br><br>And case.struct<br>
Title<br>P   LATTICE,NONEQUIV.ATOMS:  1<br>MODE OF CALC=RELA unit=ang<br> 18.897269 18.897269 18.897269 90.000000 90.000000 90.000000<br>ATOM   1: X=0.00000000 Y=0.00000000 Z=0.00000000<br>          MULT= 1          ISPLIT= 2<br>
Up         NPT=  781  R0=0.00000100 RMT=    3.5000   Z:116.0<br>LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000<br>                     0.0000000 1.0000000 0.0000000<br>                     0.0000000 0.0000000 1.0000000<br>
<br>But I receive errors after trying to run lstart:<br>error: Elemant120.inst not consistent with Z<br>edit Elemant120.inst and rerun lstart afterwards or change Z in StructGen!<br><br>If I try to put 5g states into case.inst (for element #120) lstart stops with a format error.<br>
<br>Thank you,<br>David Tompsett.<br>