[Wien] Problem in volume optimization

Gavin Abo gsabo at crimson.ua.edu
Sat May 19 19:11:38 CEST 2018


Try a setrmt reduction of 5%.  I didn't run the entire calculation, but 
it seems to start off well with that:

username at computername:~/Desktop$ cd ~/wiendata
username at computername:~/wiendata$ mkdir NiO
username at computername:~/wiendata$ cd NiO
username at computername:~/wiendata/NiO$ makestruct_lapw
...
TITLE :NiO
...
Would you like to enter  Spacegroup or Lattice (S/L)(def=S)? S
SPACE GROUP: (type ENTER or give first LETTER for a list)
give SPACE GROUP as SYMBOL or NUMBER: 225
  Info:  space group is : 225 F Fm-3m -F4;2;3

Units of lattice parameters (Bohr/Angstrom) (b/A) (def=ANG):b
Lattice PARAMETERS as a b c (1 or 3 numbers - if you specify only 1 
number, a cubic system is assumed):
7.664157
ANGLES BETWEEN lattice vectors, as alpha beta gamma (def=90.0 90.0 
90.0):90 90 90
NUMBER INEQUEVALENT ATOMS :2
ATOM  1 (ELEMENT): Ni
POSITION OF ATOM Ni as X,Y,Z (def=0 0 0) :0 0 0
ATOM  2 (ELEMENT): O
POSITION OF ATOM O as X,Y,Z (def=0 0 0) :0.5 0.5 0.5
...
SETTING UP SPHERE RadII:
SPECIFY possible REDUCTION of SPHERE RADII in % (def=0)
5 <- Here is where I sent 5% for setrmt
...
atom  Z   RMT-max   RMT
  1  28.0  1.95   1.95
  2   8.0  1.67   1.67
file    init.struct_setrmt   generated
rerun setrmt ?(y,N) (def=N):
N
The file   init.struct   has been created

   for modifications of your input you can also edit file datastruct and run
   Tmaker / setrmt init -r X    individually

username at computername:~/wiendata/NiO$ cp init.struct NiO.struct
username at computername:~/wiendata/NiO$ init_lapw -b -sp
...
   init_lapw finished ok
username at computername:~/wiendata/NiO$ x optimize
  ********************************************
    GENERATES STRUCT-FILES AND optimize.job
  PLEASE CHOOSE ONE OF THE FOLLOWING FEATURES:

  [1]  VARY VOLUME with CONSTANT RATIO A:B:C
...
1

***************************************************
Using NiO_initial.struct as template.
***************************************************

  NUMBER OF STRUCTURE CHANGES ?
5
  PLEASE ENTER VALUE            1 (IN %)
-8
  PLEASE ENTER VALUE            2 (IN %)
-6
PLEASE ENTER VALUE            3 (IN %)
0
PLEASE ENTER VALUE            4 (IN %)
6
PLEASE ENTER VALUE            5 (IN %)
8
...
  Now run   optimize.job

username at computername:~/wiendata/NiO$ sed -i 's/run_lapw/runsp_lapw/g' 
optimize.job
username at computername:~/wiendata/NiO$ ./optimize.job
running dstart in single mode
DSTART ENDS
7.1u 0.0s 0:07.21 99.7% 0+0k 0+360io 0pf+0w
...
running dstart in single mode
DSTART ENDS
8.0u 0.0s 0:08.04 99.6% 0+0k 0+376io 0pf+0w
0.0u 0.0s 0:00.03 100.0% 0+0k 0+336io 0pf+0w
clmextrapol_lapw has generated a new NiO.clmdn
hup: Command not found.
  LAPW0 END
  LAPW1 END
  LAPW1 END
  LAPW2 END
  LAPW2 END
  CORE  END
  CORE  END
  MIXER END
in cycle 2    ETEST: 0   CTEST: 0
hup: Command not found.
  LAPW0 END
  LAPW1 END
  LAPW1 END
  LAPW2 END
  LAPW2 END
  CORE  END
  CORE  END
  MIXER END
...

I have seen a = 4.169 ang [ 
https://www.researchgate.net/post/How_can_I_calculate_NiO_band_structure_with_Materials_Studio 
].  So it might be that your lattice constant a = 4.056 = 7.664157 bohr 
could be a little too short resulting in the need for a higher setrmt 
reduction to avoid overlapping RMT spheres.

On 5/19/2018 3:48 AM, Arvind Kumar wrote:
> Dear Prof. Blaha and Wien2k users,
>
>                    I am trying to volume optimization of NiO and other 
> related compounds but it shows error during run scf
>
> ERROR status in NiO_vol___-8.00
> >   stop error
>
> NN - Error
>   LAPW0 END
> hup: Command not found.
> clmextrapol_lapw has generated a new NiO.clmdn
> 0.0u 0.0s 0:00.00 0.0% 0+0k 0+336io 0pf+0w
> 2.0u 0.0s 0:02.08 99.5% 0+0k 0+368io 0pf+0w
> DSTART ENDS
> running dstart in single mode
> clmextrapol_lapw has generated a new NiO.clmup
> 0.0u 0.0s 0:00.00 0.0% 0+0k 0+336io 0pf+0w
> 2.0u 0.0s 0:02.06 100.0% 0+0k 0+368io 0pf+0w
> DSTART ENDS
> running dstart in single mode
> clmextrapol_lapw has generated a new NiO.clmsum
> 0.0u 0.0s 0:00.00 0.0% 0+0k 0+336io 0pf+0w
> 2.0u 0.0s 0:02.07 99.5% 0+0k 0+368io 0pf+0w
> DSTART ENDS
> running dstart in single mode
> 1.7u 0.0s 0:01.79 99.4% 0+0k 0+352io 0pf+0w
> DSTART ENDS
> running dstart in single mode
> 1.7u 0.0s 0:01.78 100.0% 0+0k 0+352io 0pf+0w
> DSTART ENDS
> running dstart in single mode
> 1.7u 0.0s 0:01.79 100.0% 0+0k 0+360io 0pf+0w
> DSTART ENDS
> running dstart in single mode
>
>
>
> I am running these calculation by Wien2k 17.1 version compiled in HP 
> laptop having 8GB RAM and 1 TB harddisk with mkl library with ifort+ 
> cc compiler.
>
>
> NiO struct file as below:
> NiO
> F   LATTICE,NONEQUIV.ATOMS:  2 225_Fm-3m
> MODE OF CALC=RELA unit=ang
>   7.664157  7.664157  7.664157 90.000000 90.000000 90.000000
> ATOM   1: X=0.00000000 Y=0.00000000 Z=0.00000000
>           MULT= 1          ISPLIT= 2
> Ni         NPT=  781  R0=0.00005000 RMT=    2.1100   Z: 28.000
> LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
>                      0.0000000 1.0000000 0.0000000
>                      0.0000000 0.0000000 1.0000000
> ATOM   2: X=0.50000000 Y=0.50000000 Z=0.50000000
>           MULT= 1          ISPLIT= 2
> O          NPT=  781  R0=0.00010000 RMT=    1.8100   Z: 8.000
> LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
>                      0.0000000 1.0000000 0.0000000
>                      0.0000000 0.0000000 1.0000000
>   48      NUMBER OF SYMMETRY OPERATIONS
>
>
> -- 
> Dr. Arvind Kumar
> Assistant Professor
> Department of Physics
> Atma Ram Sanatan Dharma College
> (University of Delhi)
> Dhaula Kuan, New Delhi-110021

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20180519/2493daf5/attachment.html>


More information about the Wien mailing list