"Big problem with the spin orbit in PSP5NL"

temok · Post by **temok** » Fri Jul 25, 2014 5:55 am

Dear All,

I got the following error

--- !ERROR
message: |
BIG PROBLEM WITH THE SPIN ORBIT IN PSP5NL
src_file: psp5in.F90
src_line: 430
...

when I tried to compute the SCF ground state of GaAs using [Spin-Orbit] pseudo potentials obtained with the J.L.Martins code (atom.exe) followed by the martins2abinit.f90 code. This error seems to be trigger *only* by the Arsenic pseudo potential,
apparently right after reading it and during 'pspatm: epsatm' and the splines computation.

I would appreciate any comments on this, details are found below.
Thanks,
Temok

Details:

I downloaded the Martins2abinit tool
http://www.abinit.org/downloads/psp-links/psp-links/MARTINS2ABINIT.tgz
and followed the instructions by M. Torrent therein.

I successfully reproduced the test
MARTINS2ABINIT/_atom-5.695/atom.dat
and it matched the expected result
MARTINS2ABINIT/_atom-5.695/atom_f.out

Then following the template in
MARTINS2ABINIT/_atom-5.695/atom_table.dat

I used the next input for Gallium

Code: Select all

   pg      Gallium   From Kocinski   Ecut ~ 81Ry   l=0 (maybe 1) as local
        tm2
 n=Ga c=car
       0.0       0.0       0.0       0.0       0.0       0.0
    5    3 
    4    0     2.00      0.00
    4    1     1.00      0.00 
    3    2    10.00      0.00
   2.10     2.10     2.10

and the next input for Arsenic
[which ABINIT complained about]

Code: Select all

   pg      Arsenic    Guess   Ecut ~ 20Ry   l=2 (maybe 0,1) as local
        tm2
 n=As c=car
       0.0       0.0       0.0       0.0       0.0       0.0
    6    3 
    4    0     2.00      0.00
    4    1     3.00      0.00 
    4    2     0.00      0.00
   2.10     2.10     2.10

After running atom.exe on those input files, the results were:

for gallium [next one seems OK]

Code: Select all

 atom-LDAca20-JUL-14            Gallium   From Kocinski   Ecut ~ 81Ry   l=0 
*********************

 Ga pseudopotential generation
 -----------------------------

 r e l a t i v i s t i c ! !

 correlation = ca      spin-polarized

 nuclear charge             = 31.000000
 number of core orbitals    =  7
 number of valence orbitals =  5
 electronic charge          = 31.000000
 ionic charge               =  0.000000


 input data for orbitals

  i    n    l    s     j     occ

  1    1    0   0.5   0.5    2.0000
  2    2    0   0.5   0.5    2.0000
  3    2    1  -0.5   0.5    2.0000
  4    2    1   0.5   1.5    4.0000
  5    3    0   0.5   0.5    2.0000
  6    3    1  -0.5   0.5    2.0000
  7    3    1   0.5   1.5    4.0000
  8    4    0   0.5   0.5    2.0000
  9    4    1  -0.5   0.5    0.3333
 10    4    1   0.5   1.5    0.6667
 11    3    2  -0.5   1.5    4.0000
 12    3    2   0.5   2.5    6.0000


 radial grid parameters

 r(1) = .0 , r(2) =0.101E-05 , ... , r(1138) = 118.927
 a =  6.000  b =  80.000


 n = 1  l = 0  s = 0.5
        a extr         4.043
        b extr        -0.457
        r extr         0.032
        r zero    
        r 90/99 %      0.088   0.141

 n = 2  l = 0  s = 0.5
        a extr         1.115  -2.251
        b extr        -0.163   0.045
        r extr         0.024   0.183
        r zero         0.065
        r 90/99 %      0.338   0.492

 n = 2  l = 1  s =-0.5
        a extr         2.235
        b extr         0.055
        r extr         0.145
        r zero    
        r 90/99 %      0.310   0.468

 n = 2  l = 1  s = 0.5
        a extr         2.222
        b extr        -0.110
        r extr         0.148
        r zero    
        r 90/99 %      0.314   0.474

 n = 3  l = 0  s = 0.5
        a extr         0.428  -0.744   1.382
        b extr        -0.064   0.018  -0.009
        r extr         0.024   0.150   0.566
        r zero         0.063   0.267
        r 90/99 %      0.979   1.406

 n = 3  l = 1  s =-0.5
        a extr         0.700  -1.314
        b extr         0.023  -0.008
        r extr         0.112   0.562
        r zero         0.237
        r 90/99 %      1.042   1.535

 n = 3  l = 1  s = 0.5
        a extr         0.694  -1.304
        b extr        -0.044   0.017
        r extr         0.115   0.570
        r zero         0.241
        r 90/99 %      1.055   1.554

 n = 4  l = 0  s = 0.5
        a extr         0.114  -0.195   0.319  -0.715
        b extr        -0.017   0.005  -0.003   0.001
        r extr         0.024   0.147   0.488   1.885
        r zero         0.063   0.260   0.808
        r 90/99 %      3.547   5.291

 n = 4  l = 1  s =-0.5
        a extr         0.127  -0.214   0.568
        b extr         0.004  -0.002   0.001
        r extr         0.111   0.486   2.487
        r zero         0.232   0.860
        r 90/99 %      5.358   8.195

 n = 4  l = 1  s = 0.5
        a extr         0.126  -0.211   0.562
        b extr        -0.008   0.003  -0.002
        r extr         0.113   0.493   2.537
        r zero         0.236   0.875
        r 90/99 %      5.425   8.298

 n = 3  l = 2  s =-0.5
        a extr         1.098
        b extr         0.015
        r extr         0.529
        r zero    
        r 90/99 %      1.389   2.377

 n = 3  l = 2  s = 0.5
        a extr         1.093
        b extr        -0.022
        r extr         0.532
        r zero    
        r 90/99 %      1.406   2.377

Total number of iterations needed for electron screening potential is 22

 Ga output data for orbitals
 ----------------------------

 nl    s      occ         eigenvalue    kinetic energy      pot energy

 1s   0.5    2.0000    -747.93682004       0.00000000   -1934.31093355
 2s   0.5    2.0000     -92.36470510       0.00000000    -425.51311655
 2p  -0.5    2.0000     -81.93285893       0.00000000    -423.60254808
 2p   0.5    4.0000     -79.91418014       0.00000000    -412.80667019
 3s   0.5    2.0000     -10.80531543       0.00000000    -135.68432792
 3p  -0.5    2.0000      -7.41914706       0.00000000    -127.80939159
 3p   0.5    4.0000      -7.15275657       0.00000000    -125.16339730
 4s   0.5    2.0000      -0.67297295       0.00000000     -34.87750125
 4p  -0.5    0.3333      -0.20741728       0.00000000     -24.47891848
 4p   0.5    0.6667      -0.19848577       0.00000000     -23.96211732
 3d  -0.5    4.0000      -1.44854755       0.00000000    -102.73001473
 3d   0.5    6.0000      -1.41453160       0.00000000    -101.96111903


 total energies
 --------------

 sum of eigenvalues        =    -2245.01422863
 kinetic energy from ek    =        0.00000000
 el-ion interaction energy =    -9362.29706538
 el-el  interaction energy =     1677.62985754
 vxc    correction         =     -185.45936861
 virial correction         =       10.39071222
 exchange + corr energy    =     -141.69220451
 kinetic energy from ev    =     3947.48249029
 potential energy          =    -7826.35941236
 ---------------------------------------------
 total energy              =    -3878.87692207


 virial theorem(nonrelativistic)
 --------------

 kinetic energy  *  2      =     7894.96498058
 potential energy          =    -7826.35941236
 virial correction         =       10.39071222
 ---------------------------------------------
 virial sum                =       78.99628044


 Pseudopotential generation file number  1


 Ga pseudopotential generation using the Improved Troullier and Martins method
 ------------------------------------------------------------

 nl    s    eigenvalue      rc          cdrc       delta

 1s   0.5   -0.672973    2.097979    0.455005   -1.482336
 2p  -0.5   -0.207417    2.097979    0.201114   -1.603455
 2p   0.5   -0.198486    2.097979    0.188924   -1.656638
 3d  -0.5   -1.448548    2.097979    0.982087    2.234360
 3d   0.5   -1.414532    2.097979    0.981185    2.224849

 n = 1  l = 0  s = 0.5
        a extr         0.715
        r extr         1.891
        r zero    
        r 90/99 %      3.547   5.291

 n = 2  l = 1  s =-0.5
        a extr         0.568
        r extr         2.487
        r zero    
        r 90/99 %      5.358   8.195

 n = 2  l = 1  s = 0.5
        a extr         0.562
        r extr         2.537
        r zero    
        r 90/99 %      5.425   8.298

 n = 3  l = 2  s =-0.5
        a extr         1.193
        r extr         0.829
        r zero    
        r 90/99 %      1.424   2.377

 n = 3  l = 2  s = 0.5
        a extr         1.189
        r extr         0.832
        r zero    
        r 90/99 %      1.442   2.377


 test of eigenvalues

 rcut =  s   2.47  p   2.47  p   2.47  d   2.47  d   2.47
 eval =  -0.67296  -0.20740  -0.19848  -1.44855  -1.41453

 l    vps(0)    vpsmin      at r

 s   -22.825   -22.825     0.000
 p   -25.518   -25.518     0.000
 p   -25.470   -25.470     0.000
 d   -55.470   -55.470     0.000
 d   -55.322   -55.322     0.000


 Ga output data for orbitals
 ----------------------------

 nl    s      occ         eigenvalue    kinetic energy      pot energy

 1s   0.5    2.0000      -0.67296373       0.44583822     -12.24717724
 2p  -0.5    0.3333      -0.20740235       0.55165856      -9.25392284
 2p   0.5    0.6667      -0.19847548       0.53016941      -9.07825211
 3d  -0.5    4.0000      -1.44854723      12.73592217     -33.51368385
 3d   0.5    6.0000      -1.41453172      12.63352542     -33.33824704


 total energies
 --------------

 sum of eigenvalues        =      -15.82875781
 kinetic energy from ek    =      128.17385010
 el-ion interaction energy =     -367.71538111
 el-el  interaction energy =      123.32610578
 vxc    correction         =      -22.93931931
 virial correction         =        1.03454200
 exchange + corr energy    =      -17.46312498
 kinetic energy from ev    =      128.17373105
 potential energy          =     -261.85240031
 ---------------------------------------------
 total energy              =     -133.67866926


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 The total time for the calculation is      0.27355 seconds

and the result for Arsenic [which later turns problematic]

Code: Select all

 atom-LDAca20-JUL-14            Arsenic    Guess   Ecut ~ 20Ry   l=2 (maybe 
*********************

 As pseudopotential generation
 -----------------------------

 r e l a t i v i s t i c ! !

 correlation = ca      spin-polarized

 nuclear charge             = 33.000000
 number of core orbitals    =  9
 number of valence orbitals =  5
 electronic charge          = 33.000000
 ionic charge               =  0.000000


 input data for orbitals

  i    n    l    s     j     occ

  1    1    0   0.5   0.5    2.0000
  2    2    0   0.5   0.5    2.0000
  3    2    1  -0.5   0.5    2.0000
  4    2    1   0.5   1.5    4.0000
  5    3    0   0.5   0.5    2.0000
  6    3    1  -0.5   0.5    2.0000
  7    3    1   0.5   1.5    4.0000
  8    3    2  -0.5   1.5    4.0000
  9    3    2   0.5   2.5    6.0000
 10    4    0   0.5   0.5    2.0000
 11    4    1  -0.5   0.5    1.0000
 12    4    1   0.5   1.5    2.0000
 13    4    2  -0.5   1.5    0.0000
 14    4    2   0.5   2.5    0.0000


 radial grid parameters

 r(1) = .0 , r(2) =0.945E-06 , ... , r(1143) = 118.925
 a =  6.000  b =  80.000


 n = 1  l = 0  s = 0.5
        a extr         4.173
        b extr        -0.503
        r extr         0.030
        r zero    
        r 90/99 %      0.081   0.132

 n = 2  l = 0  s = 0.5
        a extr         1.161  -2.334
        b extr        -0.181   0.051
        r extr         0.023   0.171
        r zero         0.061
        r 90/99 %      0.314   0.456

 n = 2  l = 1  s =-0.5
        a extr         2.321
        b extr         0.061
        r extr         0.135
        r zero    
        r 90/99 %      0.287   0.434

 n = 2  l = 1  s = 0.5
        a extr         2.305
        b extr        -0.122
        r extr         0.137
        r zero    
        r 90/99 %      0.291   0.440

 n = 3  l = 0  s = 0.5
        a extr         0.453  -0.782   1.445
        b extr        -0.072   0.021  -0.010
        r extr         0.022   0.139   0.523
        r zero         0.059   0.248
        r 90/99 %      0.897   1.272

 n = 3  l = 1  s =-0.5
        a extr         0.743  -1.382
        b extr         0.026  -0.010
        r extr         0.104   0.516
        r zero         0.220
        r 90/99 %      0.943   1.372

 n = 3  l = 1  s = 0.5
        a extr         0.736  -1.371
        b extr        -0.050   0.019
        r extr         0.107   0.524
        r zero         0.224
        r 90/99 %      0.954   1.406

 n = 3  l = 2  s =-0.5
        a extr         1.204
        b extr         0.019
        r extr         0.474
        r zero    
        r 90/99 %      1.137   1.828

 n = 3  l = 2  s = 0.5
        a extr         1.200
        b extr        -0.027
        r extr         0.476
        r zero    
        r 90/99 %      1.137   1.851

 n = 4  l = 0  s = 0.5
        a extr         0.138  -0.236   0.381  -0.801
        b extr        -0.022   0.006  -0.003   0.002
        r extr         0.022   0.137   0.448   1.608
        r zero         0.059   0.241   0.735
        r 90/99 %      2.904   4.225

 n = 4  l = 1  s =-0.5
        a extr         0.175  -0.288   0.673
        b extr         0.006  -0.002   0.001
        r extr         0.102   0.442   1.946
        r zero         0.214   0.770
        r 90/99 %      3.920   5.995

 n = 4  l = 1  s = 0.5
        a extr         0.172  -0.284   0.666
        b extr        -0.012   0.005  -0.002
        r extr         0.105   0.449   1.986
        r zero         0.218   0.783
        r 90/99 %      4.019   6.071

 n = 4  l = 2  s =-0.5
        a extr         0.000  -0.204
        b extr         0.000  -0.000
        r extr         0.419 111.685
        r zero         0.974
        r 90/99 %    117.448 118.925
        WARNING:  This orbital is not bound!

 n = 4  l = 2  s = 0.5
        a extr         0.000  -0.204
        b extr        -0.000   0.000
        r extr         0.421 111.685
        r zero         0.981
        r 90/99 %    117.448 118.925
        WARNING:  This orbital is not bound!

Total number of iterations needed for electron screening potential is 25



 As output data for orbitals
 ----------------------------

 nl    s      occ         eigenvalue    kinetic energy      pot energy

 1s   0.5    2.0000    -856.76483983       0.00000000   -2201.59949732
 2s   0.5    2.0000    -108.78700069       0.00000000    -487.98343109
 2p  -0.5    2.0000     -97.38414173       0.00000000    -486.22424696
 2p   0.5    4.0000     -94.71578591       0.00000000    -472.13734644
 3s   0.5    2.0000     -13.88869473       0.00000000    -157.69811785
 3p  -0.5    2.0000     -10.03942327       0.00000000    -149.76218919
 3p   0.5    4.0000      -9.67099146       0.00000000    -146.33376884
 3d  -0.5    4.0000      -3.04019661       0.00000000    -126.47505350
 3d   0.5    6.0000      -2.98657277       0.00000000    -125.55349893
 4s   0.5    2.0000      -1.07824311       0.00000000     -45.21096031
 4p  -0.5    1.0000      -0.40706755       0.00000000     -34.74903601
 4p   0.5    2.0000      -0.38639266       0.00000000     -33.94651263
 4d  -0.5    0.0000       0.00000000       0.00000000      -0.65084801
 4d   0.5    0.0000       0.00000000       0.00000000      -0.65084800


 total energies
 --------------

 sum of eigenvalues        =    -2624.69187213
 kinetic energy from ek    =        0.00000000
 el-ion interaction energy =   -10892.70461543
 el-el  interaction energy =     1934.93911766
 vxc    correction         =     -204.91729016
 virial correction         =       12.09574791
 exchange + corr energy    =     -156.71190460
 kinetic energy from ev    =     4603.05179813
 potential energy          =    -9114.47740236
 ---------------------------------------------
 total energy              =    -4511.42560423


 virial theorem(nonrelativistic)
 --------------

 kinetic energy  *  2      =     9206.10359626
 potential energy          =    -9114.47740236
 virial correction         =       12.09574791
 ---------------------------------------------
 virial sum                =      103.72194180


 Pseudopotential generation file number  1


 As pseudopotential generation using the Improved Troullier and Martins method
 ------------------------------------------------------------

 nl    s    eigenvalue      rc          cdrc       delta

 1s   0.5   -1.078243    2.097942    0.644913   -0.965386
 2p  -0.5   -0.407068    2.097942    0.402207   -0.962379
 2p   0.5   -0.386393    2.097942    0.382631   -1.030067
 3d  -0.5    0.000000    2.097942    0.000000  -11.325842
 3d   0.5    0.000000    2.097942    0.000000  -11.297860

 n = 1  l = 0  s = 0.5
        a extr         0.793
        r extr         1.655
        r zero    
        r 90/99 %      2.904   4.225

 n = 2  l = 1  s =-0.5
        a extr         0.673
        r extr         1.947
        r zero    
        r 90/99 %      3.920   5.995

 n = 2  l = 1  s = 0.5
        a extr         0.666
        r extr         1.986
        r zero    
        r 90/99 %      4.019   6.071

 n = 3  l = 2  s =-0.5
        a extr         0.205
        r extr       111.091
        r zero    
        r 90/99 %    117.448 118.925
        WARNING:  This orbital is not bound!

 n = 3  l = 2  s = 0.5
        a extr         0.205
        r extr       111.091
        r zero    
        r 90/99 %    117.448 118.925
        WARNING:  This orbital is not bound!


 test of eigenvalues

 rcut =  s   3.33  p   3.33  p   3.33  d   3.33  d   3.33
 eval =  -1.07824  -0.40706  -0.38639   0.00000   0.00000



 l    vps(0)    vpsmin      at r

 s    -3.401    -6.914     1.460
 p    -5.759    -6.717     1.321
 p    -5.537    -6.622     1.355
 d   -16.826   -16.826     0.000
 d   -17.206   -17.206     0.000


 As output data for orbitals
 ----------------------------

 nl    s      occ         eigenvalue    kinetic energy      pot energy

 1s   0.5    2.0000      -1.07824364       0.62157728      -5.16879993
 2p  -0.5    1.0000      -0.40705615       0.93946229      -4.44241103
 2p   0.5    2.0000      -0.38638650       0.90064093      -4.34953262
 3d  -0.5    0.0000       0.00000000       0.00142615      -0.09882329
 3d   0.5    0.0000       0.00000000       0.00142615      -0.09882329


 total energies
 --------------

 sum of eigenvalues        =       -3.33631642
 kinetic energy from ek    =        3.98389871
 el-ion interaction energy =      -23.47907612
 el-el  interaction energy =        9.88037926
 vxc    correction         =       -3.60195322
 virial correction         =        0.24459470
 exchange + corr energy    =       -2.76261359
 kinetic energy from ev    =        3.98395439
 potential energy          =      -16.36131044
 ---------------------------------------------
 total energy              =      -12.37735605


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 The total time for the calculation is      0.27475 seconds

then I ran
Martins2Abinit.exe so
to produce Spin-Orbit pseudopotentials.

The header for ABINIT's Gallium pseudo was

Code: Select all

   Psp for Gallium (LDA-ca - fully relativ.) using Improved Troullier and Martins method (from JL Martins psp generator)
31.  13.00  140720                               : zatom,zion,pspdat
5   2  2  0   977  0.0                           : pspcod,pspxc,lmax,lloc,mmax,r2well
 5.908270609269090E-04  1.250000000000000E-02  2 : aa,bb,pspso
0   0.   0.  1  2.097979                         : l=1,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
1   0.   0.  2  2.097979                         : l=2,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
1   0.   0.  2  2.097979                         : l=2,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
2   0.   0.  2  2.097979                         : l=3,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
2   0.   0.  2  2.097979                         : l=3,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
0.0  0.0  0.0                                    : rchrg,fchrg,qchrg
  0 =l for pseudopotential
 -1.1412668751218E+01 -1.1412668750092E+01 -1.1412668748939E+01
 -1.1412668747755E+01 -1.1412668746541E+01 -1.1412668745297E+01
 -1.1412668744021E+01 -1.1412668742713E+01 -1.1412668741371E+01
.
.
.

and the header for ABINIT's Arsenic pseudo was
[which later turns problematic]

Code: Select all

Psp for Arsenic (LDA-ca - fully relativ.) using Improved Troullier and Martins method (from JL Martins psp generator)
33.   5.00  140720                               : zatom,zion,pspdat
5   2  2  0   982  0.0                           : pspcod,pspxc,lmax,lloc,mmax,r2well
 5.550193602646721E-04  1.250000000000000E-02  2 : aa,bb,pspso
0   0.   0.  1  2.097942                         : l=1,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
1   0.   0.  2  2.097942                         : l=2,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
1   0.   0.  2  2.097942                         : l=2,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
2   0.   0.  2  2.097942                         : l=3,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
2   0.   0.  2  2.097942                         : l=3,e99.0,e99.9,nproj,rcpsp
0.0  0.0  0.0  0.0                               : rms,ekb1,ekb2,epsatm
0.0  0.0  0.0                                    : rchrg,fchrg,qchrg
  0 =l for pseudopotential
 -1.7003192967764E+00 -1.7003192947281E+00 -1.7003192926279E+00
 -1.7003192904745E+00 -1.7003192882667E+00 -1.7003192860029E+00
 -1.7003192836818E+00 -1.7003192813019E+00 -1.7003192788618E+00
 -1.7003192763600E+00 -1.7003192737948E+00 -1.7003192711646E+00
 -1.7003192684679E+00 -1.7003192657029E+00 -1.7003192628680E+00
 -1.7003192599612E+00 -1.7003192569809E+00 -1.7003192539251E+00
.
.
.

In a moment I will post my ABINIT's input and log file in this thread,
apparently I have reach some limit in the message size.

temok · Post by **temok** » Fri Jul 25, 2014 6:03 am

.
.
.

continued ...

Finally I ran abinit,

my input file [whch runs fine with HGH pseudos]

Code: Select all

# Crystalline GaAs 
# J.L.Martins pseudos w  spin-orbit 
# Computation of the band structure.
# First, a SCF density computation, then a non-SCF band structure calculation.

ndtset 2
nspinor 2
prtvol  10 

#Dataset 1 : usual self-consistent calculation
kptopt1 1          # Option for the automatic generation of k points,
                   # taking into account the symmetry
nshiftk1 4
 shiftk1  0.5 0.5 0.5
         0.5 0.0 0.0
         0.0 0.5 0.0
         0.0 0.0 0.5

ngkpt1  6 6 6  
prtden1  1         # Print the density, for use by dataset 2
toldfe1  1.0d-12    # "is not worth to go lower"

#Dataset 2 : the band structure
istwfk2   *1   #
iscf2    -2
getden2  -1
kptopt2  -5
nband2   30   

ndivsm      20

kptbounds2  0.5  0.0  0.0 # L' point
            0.0  0.0  0.0 # Gamma point
            0.0  0.5  0.5 # X point
            1/4  1/2  3/4 # W point
            3/8  3/8  3/4 # K point
            1.0  1.0  1.0 # Gamma point in another cell.

# http://www.abinit.org/documentation/helpfiles/for-v6.6/input_variables/vargs.html#kptbounds
#L-Gamma-X-W-K, 1/2 0 0  0 0 0  0 1/2 1/2  1/4 1/2 3/4  3/8 3/8 3/4
#U-L-W-X-K,     1/2 1/2 1/2  1/4 1/2 3/4  1/2 1/2 1  
#U-Gamma     3/8 3/8 3/4  0 0 0

tolwfr2  1.0d-20       # choose between 1.0d-16 and 1.0d-22.
enunit2  1             # Will output the eigenenergies in eV 


#Definition of the unit cell
# acell   3*10.60 # unrelaxed
 acell    3*10.4510070380 # relaxed, with spin

 rprim   0.0    0.5    0.5  
         0.5    0.0    0.5  
         0.5    0.5    0.0


#Definition of the atom types
ntypat    2          # There are two types of atoms
znucl 31  33         
                         

#Definition of the atoms
natom 2           # There are two atoms
typat 1 2         # Atom 1 is Ga, atoms 2 is As.
   
 xred 0     0     0   # Triplet giving the REDUCED coordinate of atom 1.
      1/4  1/4  1/4   # Triplet giving the REDUCED coordinate of atom 2.


#Definition of the planewave basis set
ecut       35         # Maximal kinetic energy cut-off, in Hartree

#Definition of the SCF procedure
nstep 200          # Maximal number of SCF cycles

and got the next log file

Code: Select all

  ABINIT 
  
  Give name for formatted input file: 
gaas.in
  Give name for formatted output file:
gaas.out
  Give root name for generic input files:
gaasi
  Give root name for generic output files:
gaaso
  Give root name for generic temporary files:
gaas.tmp

.Version 7.8.1 of ABINIT 
.(sequential version, prepared for a x86_64_linux_gnu4.4 computer) 

.Copyright (C) 1998-2014 ABINIT group . 
 ABINIT comes with ABSOLUTELY NO WARRANTY.
 It is free software, and you are welcome to redistribute it
 under certain conditions (GNU General Public License,
 see ~abinit/COPYING or http://www.gnu.org/copyleft/gpl.txt).

 ABINIT is a project of the Universite Catholique de Louvain,
 Corning Inc. and other collaborators, see ~abinit/doc/developers/contributors.txt .
 Please read ~abinit/doc/users/acknowledgments.html for suggested
 acknowledgments of the ABINIT effort.
 For more information, see http://www.abinit.org .

.Starting date : Mon 21 Jul 2014.
- ( at 01h17 )
  

 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

 === Build Information === 
  Version       : 7.8.1
  Build target  : x86_64_linux_gnu4.4
  Build date    : 20140702

 === Compiler Suite === 
  C compiler       : gnu4.4
  C++ compiler     : gnu4.4
  Fortran compiler : gnu4.4
  CFLAGS           :  -g -O2 -mtune=native -march=native
  CXXFLAGS         :  -g -O2 -mtune=native -march=native
  FCFLAGS          :  -g -ffree-line-length-none
  FC_LDFLAGS       : 

 === Optimizations === 
  Debug level        : basic
  Optimization level : standard
  Architecture       : amd_opteron

 === Multicore === 
  Parallel build : no
  Parallel I/O   : no
  openMP support : no
  GPU support    : no

 === Connectors / Fallbacks === 
  Connectors on : yes
  Fallbacks on  : yes
  DFT flavor    : libxc-fallback+atompaw-fallback+bigdft-fallback+wannier90-fallback
  FFT flavor    : none
  LINALG flavor : netlib
  MATH flavor   : none
  TIMER flavor  : abinit
  TRIO flavor   : netcdf-fallback+etsf_io-fallback

 === Experimental features === 
  Bindings            : @enable_bindings@
  Exports             : no
  GW double-precision : no

 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++


 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 Default optimizations:
   -O2 -mtune=native -march=native


 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++


 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 CPP options activated during the build:

                    CC_GNU                   CXX_GNU                    FC_GNU
 
          HAVE_DFT_ATOMPAW           HAVE_DFT_BIGDFT            HAVE_DFT_LIBXC
 
        HAVE_DFT_WANNIER90 HAVE_FC_ALLOCATABLE_DT...  HAVE_FC_COMMAND_ARGUMENT
 
           HAVE_FC_CPUTIME              HAVE_FC_EXIT             HAVE_FC_FLUSH
 
             HAVE_FC_GAMMA            HAVE_FC_GETENV          HAVE_FC_INT_QUAD
 
             HAVE_FC_IOMSG     HAVE_FC_ISO_C_BINDING        HAVE_FC_LONG_LINES
 
           HAVE_FC_PRIVATE         HAVE_FC_PROTECTED         HAVE_FC_STREAM_IO
 
            HAVE_FC_SYSTEM          HAVE_FORTRAN2003        HAVE_LIBPAW_ABINIT
 
               HAVE_LINALG        HAVE_LINALG_SERIAL        HAVE_MEM_PROFILING
 
                HAVE_NUMPY             HAVE_OS_LINUX                HAVE_TIMER
 
         HAVE_TIMER_ABINIT         HAVE_TIMER_SERIAL         HAVE_TRIO_ETSF_IO
 
          HAVE_TRIO_NETCDF              USE_MACROAVE                            
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

- input  file    -> gaas.in
- output file    -> gaas.out
- root for input  files -> gaasi
- root for output files -> gaaso

-instrng :    84 lines of input have been read from file gaas.in


 iofn2 : Please give name of formatted atomic psp file
 iofn2 : for atom type   1 , psp file is /home/temok/AtomicPotentialFiles/martins2abinit/31_13_Ga_TM_soc_Kocinski

 iofn2 : Please give name of formatted atomic psp file
 iofn2 : for atom type   2 , psp file is /home/temok/AtomicPotentialFiles/martins2abinit/33_05_As_TM_soc
  read the values zionpsp= 13.0 , pspcod=   5 , lmax=   2
  read the values zionpsp=  5.0 , pspcod=   5 , lmax=   2

 inpspheads: deduce mpsang  =   3, n1xccc  =   0.

 invars1m : enter jdtset=1
 invars1 : treat image number     1

 symlatt : the Bravais lattice is cF (face-centered cubic)
  xred   is defined in input file
 ingeo : takes atomic coordinates from input array xred

 symlatt : the Bravais lattice is cF (face-centered cubic)

 symlatt : the Bravais lattice is cF (face-centered cubic)
 symspgr : spgroup= 216  F-4 3 m   (=Td^2)
       Simple Lattice Grid

 invars1m : enter jdtset=2
 invars1 : treat image number     1

 symlatt : the Bravais lattice is cF (face-centered cubic)
  xred   is defined in input file
 ingeo : takes atomic coordinates from input array xred

 symlatt : the Bravais lattice is cF (face-centered cubic)

 symlatt : the Bravais lattice is cF (face-centered cubic)
 symspgr : spgroup= 216  F-4 3 m   (=Td^2)
 Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
 R(1)=  0.0000000  5.2255035  5.2255035  G(1)= -0.0956846  0.0956846  0.0956846
 R(2)=  5.2255035  0.0000000  5.2255035  G(2)=  0.0956846 -0.0956846  0.0956846
 R(3)=  5.2255035  5.2255035  0.0000000  G(3)=  0.0956846  0.0956846 -0.0956846
 Unit cell volume ucvol=  2.8537402E+02 bohr^3
 Unit cell volume ucvol=  2.8537402E+02 bohr^3
 Angles (23,13,12)=  6.00000000E+01  6.00000000E+01  6.00000000E+01 degrees
 Angles (23,13,12)=  6.00000000E+01  6.00000000E+01  6.00000000E+01 degrees
  8.28652588822836600E-002  9.56845590443090033E-002  4.78422795221545016E-002  3.38296002775377441E-002  0.10148880083261325     

 mknormpath : total number of points on the path :      215

 Number of divisions for each segment of the normalized path :
 0.50000 0.00000 0.00000 ==>  0.00000 0.00000 0.00000 ( ndiv :    49 )
 0.00000 0.00000 0.00000 ==>  0.00000 0.50000 0.50000 ( ndiv :    57 )
 0.00000 0.50000 0.50000 ==>  0.25000 0.50000 0.75000 ( ndiv :    28 )
 0.25000 0.50000 0.75000 ==>  0.37500 0.37500 0.75000 ( ndiv :    20 )
 0.37500 0.37500 0.75000 ==>  1.00000 1.00000 1.00000 ( ndiv :    60 )


 invars2: take the default value of fband=  1.25000000E-01
       Simple Lattice Grid
  dtset%nelect=   18.000000000000000     
 chkneu : initialized the occupation numbers for occopt=    1
    spin-unpolarized case :
  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00
  1.00  1.00  1.00  1.00  1.00  1.00  0.00  0.00
 Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
 R(1)=  0.0000000  5.2255035  5.2255035  G(1)= -0.0956846  0.0956846  0.0956846
 R(2)=  5.2255035  0.0000000  5.2255035  G(2)=  0.0956846 -0.0956846  0.0956846
 R(3)=  5.2255035  5.2255035  0.0000000  G(3)=  0.0956846  0.0956846 -0.0956846
 Unit cell volume ucvol=  2.8537402E+02 bohr^3
 Unit cell volume ucvol=  2.8537402E+02 bohr^3
 Angles (23,13,12)=  6.00000000E+01  6.00000000E+01  6.00000000E+01 degrees
 Angles (23,13,12)=  6.00000000E+01  6.00000000E+01  6.00000000E+01 degrees
  8.28652588822836600E-002  9.56845590443090033E-002  4.78422795221545016E-002  3.38296002775377441E-002  0.10148880083261325     

 mknormpath : total number of points on the path :      215

 Number of divisions for each segment of the normalized path :
 0.50000 0.00000 0.00000 ==>  0.00000 0.00000 0.00000 ( ndiv :    49 )
 0.00000 0.00000 0.00000 ==>  0.00000 0.50000 0.50000 ( ndiv :    57 )
 0.00000 0.50000 0.50000 ==>  0.25000 0.50000 0.75000 ( ndiv :    28 )
 0.25000 0.50000 0.75000 ==>  0.37500 0.37500 0.75000 ( ndiv :    20 )
 0.37500 0.37500 0.75000 ==>  1.00000 1.00000 1.00000 ( ndiv :    60 )


  dtset%nelect=   18.000000000000000     
 chkneu : initialized the occupation numbers for occopt=    1
    spin-unpolarized case :
  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00
  1.00  1.00  1.00  1.00  1.00  1.00  0.00  0.00  0.00  0.00  0.00  0.00
  0.00  0.00  0.00  0.00  0.00  0.00
mpi_setup: mkmem  undefined in the input file.Use default mkmem  = nkpt
 mpi_setup: With nkpt_me=28 and mkmem  = 28, ground state wf handled in core.
mpi_setup: mkqmem undefined in the input file.Use default mkqmem = nkpt
 mpi_setup: With nkpt_me=28 and mkqmem = 28, ground state wf handled in core.
mpi_setup: mk1mem undefined in the input file.Use default mk1mem = nkpt
 mpi_setup: With nkpt_me=28 and mk1mem = 28, ground state wf handled in core.

--- !COMMENT
message: |
    The second and third dimension of the FFT grid: 0, 0
    were imposed to be multiple of the number of processors for the FFT: 1
src_file: getng.F90
src_line: 204
...

 For input ecut=  3.500000E+01 best grid ngfft=      40      40      40
       max ecut=  3.614460E+01

 ==== FFT mesh ====
  FFT mesh divisions ........................    40   40   40
  Augmented FFT divisions ...................    41   41   40
  FFT algorithm .............................   112
  FFT cache size ............................    16
 getmpw: optimal value of mpw=    2832
mpi_setup: mkmem  undefined in the input file.Use default mkmem  = nkpt
 mpi_setup: With nkpt_me=215 and mkmem  = 215, ground state wf handled in core.
mpi_setup: mkqmem undefined in the input file.Use default mkqmem = nkpt
 mpi_setup: With nkpt_me=215 and mkqmem = 215, ground state wf handled in core.
mpi_setup: mk1mem undefined in the input file.Use default mk1mem = nkpt
 mpi_setup: With nkpt_me=215 and mk1mem = 215, ground state wf handled in core.
 For input ecut=  3.500000E+01 best grid ngfft=      40      40      40
       max ecut=  3.614460E+01

 ==== FFT mesh ====
  FFT mesh divisions ........................    40   40   40
  Augmented FFT divisions ...................    41   41   40
  FFT algorithm .............................   112
  FFT cache size ............................    16
 getmpw: optimal value of mpw=    2842

 DATASET    1 : space group F-4 3 m (#216); Bravais cF (face-center cubic)

 getdim_nloc : deduce lmnmax  =  17, lnmax  =   5,
                      lmnmaxso=  33, lnmaxso=   9.
memory : analysis of memory needs
================================================================================
 Values of the parameters that define the memory need for DATASET  1.
   intxc =         0  ionmov =         0    iscf =         7 xclevel =         1
  lmnmax =         9   lnmax =         9   mband =        20  mffmem =         1
P  mgfft =        40   mkmem =        28 mpssoang=         5     mpw =      2832
  mqgrid =      3001   natom =         2    nfft =     64000    nkpt =        28
  nloalg =         4  nspden =         1 nspinor =         2  nsppol =         1
    nsym =        24  n1xccc =         0  ntypat =         2  occopt =         1
================================================================================
P This job should need less than                      69.668 Mbytes of memory.
  Rough estimation (10% accuracy) of disk space for files :
_ WF disk file :     24.201 Mbytes ; DEN or POT disk file :      0.490 Mbytes.
================================================================================

 Biggest array : cg(disk), with     24.2012 MBytes.
 memana : allocated an array of     24.201 Mbytes, for testing purposes.
 memana : allocated      69.668 Mbytes, for testing purposes. 
 The job will continue.

 DATASET    2 : space group F-4 3 m (#216); Bravais cF (face-center cubic)

 getdim_nloc : deduce lmnmax  =  17, lnmax  =   5,
                      lmnmaxso=  33, lnmaxso=   9.
memory : analysis of memory needs
================================================================================
 Values of the parameters that define the memory need for DATASET  2.
   intxc =         0  ionmov =         0    iscf =        -2 xclevel =         1
  lmnmax =         9   lnmax =         9   mband =        30  mffmem =         1
P  mgfft =        40   mkmem =       215 mpssoang=         5     mpw =      2842
  mqgrid =      3001   natom =         2    nfft =     64000    nkpt =       215
  nloalg =         4  nspden =         1 nspinor =         2  nsppol =         1
    nsym =        24  n1xccc =         0  ntypat =         2  occopt =         1
================================================================================
P This job should need less than                     579.876 Mbytes of memory.
  Rough estimation (10% accuracy) of disk space for files :
_ WF disk file :    279.709 Mbytes ; DEN or POT disk file :      0.490 Mbytes.
================================================================================

 Biggest array : cg(disk), with    279.7093 MBytes.
 memana : allocated an array of    279.709 Mbytes, for testing purposes.
 memana : allocated     579.876 Mbytes, for testing purposes. 
 The job will continue.
--------------------------------------------------------------------------------
------------- Echo of variables that govern the present computation ------------
--------------------------------------------------------------------------------
-
- outvars: echo of selected default values                                      
-   accesswff0 =  0 , fftalg0 =112 , wfoptalg0 =  0
-
- outvars: echo of global parameters not present in the input file              
-  max_nthreads =    0
-
 -outvars: echo values of preprocessed input variables --------
 
 These variables are accessible in NetCDF format (gaaso_OUT.nc)

            acell      1.0451007038E+01  1.0451007038E+01  1.0451007038E+01 Bohr
              amu      6.97230000E+01  7.49215900E+01
             ecut      3.50000000E+01 Hartree
           enunit1          0
           enunit2          1
-          fftalg         112
           getden1          0
           getden2         -1
             iscf1          7
             iscf2         -2
           istwfk2       1    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    1
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    1    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0    0    0    0    0    0
                         0    0    0    0    0
              ixc           2
           jdtset        1    2
              kpt1    -8.33333333E-02 -1.66666667E-01  0.00000000E+00
                      -8.33333333E-02 -3.33333333E-01  0.00000000E+00
                      -1.66666667E-01 -2.50000000E-01  0.00000000E+00
                      -8.33333333E-02 -2.50000000E-01  8.33333333E-02
                      -8.33333333E-02  5.00000000E-01  0.00000000E+00
                      -1.66666667E-01 -4.16666667E-01  0.00000000E+00
                      -8.33333333E-02 -4.16666667E-01  8.33333333E-02
                      -2.50000000E-01 -3.33333333E-01  0.00000000E+00
                      -1.66666667E-01 -3.33333333E-01  8.33333333E-02
                      -8.33333333E-02 -3.33333333E-01  1.66666667E-01
                      -8.33333333E-02  3.33333333E-01  0.00000000E+00
                      -1.66666667E-01  4.16666667E-01  0.00000000E+00
                      -2.50000000E-01  5.00000000E-01  0.00000000E+00
                      -1.66666667E-01  5.00000000E-01  8.33333333E-02
                      -3.33333333E-01 -4.16666667E-01  0.00000000E+00
                      -2.50000000E-01 -4.16666667E-01  8.33333333E-02
                      -1.66666667E-01 -4.16666667E-01  1.66666667E-01
                      -8.33333333E-02 -4.16666667E-01  2.50000000E-01
                      -8.33333333E-02  1.66666667E-01  0.00000000E+00
                      -1.66666667E-01  2.50000000E-01  0.00000000E+00
                      -2.50000000E-01  3.33333333E-01  0.00000000E+00
                      -3.33333333E-01  4.16666667E-01  0.00000000E+00
                      -4.16666667E-01  5.00000000E-01  0.00000000E+00
                      -3.33333333E-01  5.00000000E-01  8.33333333E-02
                      -2.50000000E-01  5.00000000E-01  1.66666667E-01
                      -8.33333333E-02  0.00000000E+00  0.00000000E+00
                      -2.50000000E-01  0.00000000E+00  0.00000000E+00
                      -4.16666667E-01  0.00000000E+00  0.00000000E+00
              kpt2     5.00000000E-01  0.00000000E+00  0.00000000E+00
                       4.89795918E-01  0.00000000E+00  0.00000000E+00
                       4.79591837E-01  0.00000000E+00  0.00000000E+00
                       4.69387755E-01  0.00000000E+00  0.00000000E+00
                       4.59183673E-01  0.00000000E+00  0.00000000E+00
                       4.48979592E-01  0.00000000E+00  0.00000000E+00
                       4.38775510E-01  0.00000000E+00  0.00000000E+00
                       4.28571429E-01  0.00000000E+00  0.00000000E+00
                       4.18367347E-01  0.00000000E+00  0.00000000E+00
                       4.08163265E-01  0.00000000E+00  0.00000000E+00
                       3.97959184E-01  0.00000000E+00  0.00000000E+00
                       3.87755102E-01  0.00000000E+00  0.00000000E+00
                       3.77551020E-01  0.00000000E+00  0.00000000E+00
                       3.67346939E-01  0.00000000E+00  0.00000000E+00
                       3.57142857E-01  0.00000000E+00  0.00000000E+00
                       3.46938776E-01  0.00000000E+00  0.00000000E+00
                       3.36734694E-01  0.00000000E+00  0.00000000E+00
                       3.26530612E-01  0.00000000E+00  0.00000000E+00
                       3.16326531E-01  0.00000000E+00  0.00000000E+00
                       3.06122449E-01  0.00000000E+00  0.00000000E+00
                       2.95918367E-01  0.00000000E+00  0.00000000E+00
                       2.85714286E-01  0.00000000E+00  0.00000000E+00
                       2.75510204E-01  0.00000000E+00  0.00000000E+00
                       2.65306122E-01  0.00000000E+00  0.00000000E+00
                       2.55102041E-01  0.00000000E+00  0.00000000E+00
                       2.44897959E-01  0.00000000E+00  0.00000000E+00
                       2.34693878E-01  0.00000000E+00  0.00000000E+00
                       2.24489796E-01  0.00000000E+00  0.00000000E+00
                       2.14285714E-01  0.00000000E+00  0.00000000E+00
                       2.04081633E-01  0.00000000E+00  0.00000000E+00
                       1.93877551E-01  0.00000000E+00  0.00000000E+00
                       1.83673469E-01  0.00000000E+00  0.00000000E+00
                       1.73469388E-01  0.00000000E+00  0.00000000E+00
                       1.63265306E-01  0.00000000E+00  0.00000000E+00
                       1.53061224E-01  0.00000000E+00  0.00000000E+00
                       1.42857143E-01  0.00000000E+00  0.00000000E+00
                       1.32653061E-01  0.00000000E+00  0.00000000E+00
                       1.22448980E-01  0.00000000E+00  0.00000000E+00
                       1.12244898E-01  0.00000000E+00  0.00000000E+00
                       1.02040816E-01  0.00000000E+00  0.00000000E+00
                       9.18367347E-02  0.00000000E+00  0.00000000E+00
                       8.16326531E-02  0.00000000E+00  0.00000000E+00
                       7.14285714E-02  0.00000000E+00  0.00000000E+00
                       6.12244898E-02  0.00000000E+00  0.00000000E+00
                       5.10204082E-02  0.00000000E+00  0.00000000E+00
                       4.08163265E-02  0.00000000E+00  0.00000000E+00
                       3.06122449E-02  0.00000000E+00  0.00000000E+00
                       2.04081633E-02  0.00000000E+00  0.00000000E+00
                       1.02040816E-02  0.00000000E+00  0.00000000E+00
                       0.00000000E+00  0.00000000E+00  0.00000000E+00
                       0.00000000E+00  8.77192982E-03  8.77192982E-03
                       0.00000000E+00  1.75438596E-02  1.75438596E-02
                       0.00000000E+00  2.63157895E-02  2.63157895E-02
                       0.00000000E+00  3.50877193E-02  3.50877193E-02
                       0.00000000E+00  4.38596491E-02  4.38596491E-02
                       0.00000000E+00  5.26315789E-02  5.26315789E-02
                       0.00000000E+00  6.14035088E-02  6.14035088E-02
                       0.00000000E+00  7.01754386E-02  7.01754386E-02
                       0.00000000E+00  7.89473684E-02  7.89473684E-02
                       0.00000000E+00  8.77192982E-02  8.77192982E-02
                       0.00000000E+00  9.64912281E-02  9.64912281E-02
                       0.00000000E+00  1.05263158E-01  1.05263158E-01
                       0.00000000E+00  1.14035088E-01  1.14035088E-01
                       0.00000000E+00  1.22807018E-01  1.22807018E-01
                       0.00000000E+00  1.31578947E-01  1.31578947E-01
                       0.00000000E+00  1.40350877E-01  1.40350877E-01
                       0.00000000E+00  1.49122807E-01  1.49122807E-01
                       0.00000000E+00  1.57894737E-01  1.57894737E-01
                       0.00000000E+00  1.66666667E-01  1.66666667E-01
                       0.00000000E+00  1.75438596E-01  1.75438596E-01
                       0.00000000E+00  1.84210526E-01  1.84210526E-01
                       0.00000000E+00  1.92982456E-01  1.92982456E-01
                       0.00000000E+00  2.01754386E-01  2.01754386E-01
                       0.00000000E+00  2.10526316E-01  2.10526316E-01
                       0.00000000E+00  2.19298246E-01  2.19298246E-01
                       0.00000000E+00  2.28070175E-01  2.28070175E-01
                       0.00000000E+00  2.36842105E-01  2.36842105E-01
                       0.00000000E+00  2.45614035E-01  2.45614035E-01
                       0.00000000E+00  2.54385965E-01  2.54385965E-01
                       0.00000000E+00  2.63157895E-01  2.63157895E-01
                       0.00000000E+00  2.71929825E-01  2.71929825E-01
                       0.00000000E+00  2.80701754E-01  2.80701754E-01
                       0.00000000E+00  2.89473684E-01  2.89473684E-01
                       0.00000000E+00  2.98245614E-01  2.98245614E-01
                       0.00000000E+00  3.07017544E-01  3.07017544E-01
                       0.00000000E+00  3.15789474E-01  3.15789474E-01
                       0.00000000E+00  3.24561404E-01  3.24561404E-01
                       0.00000000E+00  3.33333333E-01  3.33333333E-01
                       0.00000000E+00  3.42105263E-01  3.42105263E-01
                       0.00000000E+00  3.50877193E-01  3.50877193E-01
                       0.00000000E+00  3.59649123E-01  3.59649123E-01
                       0.00000000E+00  3.68421053E-01  3.68421053E-01
                       0.00000000E+00  3.77192982E-01  3.77192982E-01
                       0.00000000E+00  3.85964912E-01  3.85964912E-01
                       0.00000000E+00  3.94736842E-01  3.94736842E-01
                       0.00000000E+00  4.03508772E-01  4.03508772E-01
                       0.00000000E+00  4.12280702E-01  4.12280702E-01
                       0.00000000E+00  4.21052632E-01  4.21052632E-01
                       0.00000000E+00  4.29824561E-01  4.29824561E-01
                       0.00000000E+00  4.38596491E-01  4.38596491E-01
                       0.00000000E+00  4.47368421E-01  4.47368421E-01
                       0.00000000E+00  4.56140351E-01  4.56140351E-01
                       0.00000000E+00  4.64912281E-01  4.64912281E-01
                       0.00000000E+00  4.73684211E-01  4.73684211E-01
                       0.00000000E+00  4.82456140E-01  4.82456140E-01
                       0.00000000E+00  4.91228070E-01  4.91228070E-01
                       0.00000000E+00  5.00000000E-01  5.00000000E-01
                       8.92857143E-03  5.00000000E-01  5.08928571E-01
                       1.78571429E-02  5.00000000E-01  5.17857143E-01
                       2.67857143E-02  5.00000000E-01  5.26785714E-01
                       3.57142857E-02  5.00000000E-01  5.35714286E-01
                       4.46428571E-02  5.00000000E-01  5.44642857E-01
                       5.35714286E-02  5.00000000E-01  5.53571429E-01
                       6.25000000E-02  5.00000000E-01  5.62500000E-01
                       7.14285714E-02  5.00000000E-01  5.71428571E-01
                       8.03571429E-02  5.00000000E-01  5.80357143E-01
                       8.92857143E-02  5.00000000E-01  5.89285714E-01
                       9.82142857E-02  5.00000000E-01  5.98214286E-01
                       1.07142857E-01  5.00000000E-01  6.07142857E-01
                       1.16071429E-01  5.00000000E-01  6.16071429E-01
                       1.25000000E-01  5.00000000E-01  6.25000000E-01
                       1.33928571E-01  5.00000000E-01  6.33928571E-01
                       1.42857143E-01  5.00000000E-01  6.42857143E-01
                       1.51785714E-01  5.00000000E-01  6.51785714E-01
                       1.60714286E-01  5.00000000E-01  6.60714286E-01
                       1.69642857E-01  5.00000000E-01  6.69642857E-01
                       1.78571429E-01  5.00000000E-01  6.78571429E-01
                       1.87500000E-01  5.00000000E-01  6.87500000E-01
                       1.96428571E-01  5.00000000E-01  6.96428571E-01
                       2.05357143E-01  5.00000000E-01  7.05357143E-01
                       2.14285714E-01  5.00000000E-01  7.14285714E-01
                       2.23214286E-01  5.00000000E-01  7.23214286E-01
                       2.32142857E-01  5.00000000E-01  7.32142857E-01
                       2.41071429E-01  5.00000000E-01  7.41071429E-01
                       2.50000000E-01  5.00000000E-01  7.50000000E-01
                       2.56250000E-01  4.93750000E-01  7.50000000E-01
                       2.62500000E-01  4.87500000E-01  7.50000000E-01
                       2.68750000E-01  4.81250000E-01  7.50000000E-01
                       2.75000000E-01  4.75000000E-01  7.50000000E-01
                       2.81250000E-01  4.68750000E-01  7.50000000E-01
                       2.87500000E-01  4.62500000E-01  7.50000000E-01
                       2.93750000E-01  4.56250000E-01  7.50000000E-01
                       3.00000000E-01  4.50000000E-01  7.50000000E-01
                       3.06250000E-01  4.43750000E-01  7.50000000E-01
                       3.12500000E-01  4.37500000E-01  7.50000000E-01
                       3.18750000E-01  4.31250000E-01  7.50000000E-01
                       3.25000000E-01  4.25000000E-01  7.50000000E-01
                       3.31250000E-01  4.18750000E-01  7.50000000E-01
                       3.37500000E-01  4.12500000E-01  7.50000000E-01
                       3.43750000E-01  4.06250000E-01  7.50000000E-01
                       3.50000000E-01  4.00000000E-01  7.50000000E-01
                       3.56250000E-01  3.93750000E-01  7.50000000E-01
                       3.62500000E-01  3.87500000E-01  7.50000000E-01
                       3.68750000E-01  3.81250000E-01  7.50000000E-01
                       3.75000000E-01  3.75000000E-01  7.50000000E-01
                       3.85416667E-01  3.85416667E-01  7.54166667E-01
                       3.95833333E-01  3.95833333E-01  7.58333333E-01
                       4.06250000E-01  4.06250000E-01  7.62500000E-01
                       4.16666667E-01  4.16666667E-01  7.66666667E-01
                       4.27083333E-01  4.27083333E-01  7.70833333E-01
                       4.37500000E-01  4.37500000E-01  7.75000000E-01
                       4.47916667E-01  4.47916667E-01  7.79166667E-01
                       4.58333333E-01  4.58333333E-01  7.83333333E-01
                       4.68750000E-01  4.68750000E-01  7.87500000E-01
                       4.79166667E-01  4.79166667E-01  7.91666667E-01
                       4.89583333E-01  4.89583333E-01  7.95833333E-01
                       5.00000000E-01  5.00000000E-01  8.00000000E-01
                       5.10416667E-01  5.10416667E-01  8.04166667E-01
                       5.20833333E-01  5.20833333E-01  8.08333333E-01
                       5.31250000E-01  5.31250000E-01  8.12500000E-01
                       5.41666667E-01  5.41666667E-01  8.16666667E-01
                       5.52083333E-01  5.52083333E-01  8.20833333E-01
                       5.62500000E-01  5.62500000E-01  8.25000000E-01
                       5.72916667E-01  5.72916667E-01  8.29166667E-01
                       5.83333333E-01  5.83333333E-01  8.33333333E-01
                       5.93750000E-01  5.93750000E-01  8.37500000E-01
                       6.04166667E-01  6.04166667E-01  8.41666667E-01
                       6.14583333E-01  6.14583333E-01  8.45833333E-01
                       6.25000000E-01  6.25000000E-01  8.50000000E-01
                       6.35416667E-01  6.35416667E-01  8.54166667E-01
                       6.45833333E-01  6.45833333E-01  8.58333333E-01
                       6.56250000E-01  6.56250000E-01  8.62500000E-01
                       6.66666667E-01  6.66666667E-01  8.66666667E-01
                       6.77083333E-01  6.77083333E-01  8.70833333E-01
                       6.87500000E-01  6.87500000E-01  8.75000000E-01
                       6.97916667E-01  6.97916667E-01  8.79166667E-01
                       7.08333333E-01  7.08333333E-01  8.83333333E-01
                       7.18750000E-01  7.18750000E-01  8.87500000E-01
                       7.29166667E-01  7.29166667E-01  8.91666667E-01
                       7.39583333E-01  7.39583333E-01  8.95833333E-01
                       7.50000000E-01  7.50000000E-01  9.00000000E-01
                       7.60416667E-01  7.60416667E-01  9.04166667E-01
                       7.70833333E-01  7.70833333E-01  9.08333333E-01
                       7.81250000E-01  7.81250000E-01  9.12500000E-01
                       7.91666667E-01  7.91666667E-01  9.16666667E-01
                       8.02083333E-01  8.02083333E-01  9.20833333E-01
                       8.12500000E-01  8.12500000E-01  9.25000000E-01
                       8.22916667E-01  8.22916667E-01  9.29166667E-01
                       8.33333333E-01  8.33333333E-01  9.33333333E-01
                       8.43750000E-01  8.43750000E-01  9.37500000E-01
                       8.54166667E-01  8.54166667E-01  9.41666667E-01
                       8.64583333E-01  8.64583333E-01  9.45833333E-01
                       8.75000000E-01  8.75000000E-01  9.50000000E-01
                       8.85416667E-01  8.85416667E-01  9.54166667E-01
                       8.95833333E-01  8.95833333E-01  9.58333333E-01
                       9.06250000E-01  9.06250000E-01  9.62500000E-01
                       9.16666667E-01  9.16666667E-01  9.66666667E-01
                       9.27083333E-01  9.27083333E-01  9.70833333E-01
                       9.37500000E-01  9.37500000E-01  9.75000000E-01
                       9.47916667E-01  9.47916667E-01  9.79166667E-01
                       9.58333333E-01  9.58333333E-01  9.83333333E-01
                       9.68750000E-01  9.68750000E-01  9.87500000E-01
                       9.79166667E-01  9.79166667E-01  9.91666667E-01
                       9.89583333E-01  9.89583333E-01  9.95833333E-01
                       1.00000000E+00  1.00000000E+00  1.00000000E+00
           kptopt1          1
           kptopt2         -5
         kptrlatt        6   -6    6     -6    6    6     -6   -6    6
          kptrlen1     6.27060422E+01
          kptrlen2     3.00000000E+01
P           mkmem1         28
P           mkmem2        215
            natom           2
            nband1         20
            nband2         30
           nbdbuf1          0
           nbdbuf2          4
           ndtset           2
            ngfft          40      40      40
             nkpt1         28
             nkpt2        215
          nspinor           2
            nstep         200
             nsym          24
           ntypat           2
              occ      1.000000  1.000000  1.000000  1.000000  1.000000  1.000000
                       1.000000  1.000000  1.000000  1.000000  1.000000  1.000000
                       1.000000  1.000000  1.000000  1.000000  1.000000  1.000000
                       0.000000  0.000000
           prtvol          10
            rprim      0.0000000000E+00  5.0000000000E-01  5.0000000000E-01
                       5.0000000000E-01  0.0000000000E+00  5.0000000000E-01
                       5.0000000000E-01  5.0000000000E-01  0.0000000000E+00
           shiftk      5.00000000E-01  5.00000000E-01  5.00000000E-01
          spgroup         216
           symrel      1  0  0   0  1  0   0  0  1       0 -1  1   0 -1  0   1 -1  0
                      -1  0  0  -1  0  1  -1  1  0       0  1 -1   1  0 -1   0  0 -1
                      -1  0  0  -1  1  0  -1  0  1       0 -1  1   1 -1  0   0 -1  0
                       1  0  0   0  0  1   0  1  0       0  1 -1   0  0 -1   1  0 -1
                      -1  0  1  -1  1  0  -1  0  0       0 -1  0   1 -1  0   0 -1  1
                       1  0 -1   0  0 -1   0  1 -1       0  1  0   0  0  1   1  0  0
                       1  0 -1   0  1 -1   0  0 -1       0 -1  0   0 -1  1   1 -1  0
                      -1  0  1  -1  0  0  -1  1  0       0  1  0   1  0  0   0  0  1
                       0  0 -1   0  1 -1   1  0 -1       1 -1  0   0 -1  1   0 -1  0
                       0  0  1   1  0  0   0  1  0      -1  1  0  -1  0  0  -1  0  1
                       0  0  1   0  1  0   1  0  0       1 -1  0   0 -1  0   0 -1  1
                       0  0 -1   1  0 -1   0  1 -1      -1  1  0  -1  0  1  -1  0  0
           toldfe1     1.00000000E-12 Hartree
           toldfe2     0.00000000E+00 Hartree
           tolwfr1     0.00000000E+00
           tolwfr2     1.00000000E-20
            typat      1  2
              wtk1       0.02778    0.02778    0.02778    0.05556    0.02778    0.02778
                         0.05556    0.02778    0.05556    0.05556    0.02778    0.02778
                         0.02778    0.05556    0.02778    0.05556    0.05556    0.05556
                         0.02778    0.02778    0.02778    0.02778    0.02778    0.05556
                         0.05556    0.00926    0.00926    0.00926
              wtk2       1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000    1.00000
                         1.00000    1.00000    1.00000    1.00000    1.00000
           xangst      0.0000000000E+00  0.0000000000E+00  0.0000000000E+00
                       1.3826086828E+00  1.3826086828E+00  1.3826086828E+00
            xcart      0.0000000000E+00  0.0000000000E+00  0.0000000000E+00
                       2.6127517595E+00  2.6127517595E+00  2.6127517595E+00
             xred      0.0000000000E+00  0.0000000000E+00  0.0000000000E+00
                       2.5000000000E-01  2.5000000000E-01  2.5000000000E-01
            znucl       31.00000   33.00000

================================================================================

 chkinp: machine precision is   2.2204460492503131E-16

 chkinp: Checking input parameters for consistency, jdtset=   1.

 chkinp: Checking input parameters for consistency, jdtset=   2.
 DATA TYPE INFORMATION: 
 REAL:      Data type name: REAL(DP) 
            Kind value:      8
            Precision:      15
            Smallest nonnegligible quantity relative to 1:  0.22204460E-15
            Smallest positive number:                       0.22250739-307
            Largest representable number:                   0.17976931+309
 INTEGER:   Data type name: INTEGER(default) 
            Kind value: 4
            Bit size:   32
            Largest representable number: 2147483647
 LOGICAL:   Data type name: LOGICAL 
            Kind value: 4
 CHARACTER: Data type name: CHARACTER             Kind value: 1

  ==== OpenMP parallelism is OFF ====
 
  MPI-IO support is OFF

================================================================================
== DATASET  1 ==================================================================
-   nproc =    1


 getdim_nloc : deduce lmnmax  =  17, lnmax  =   5,
                      lmnmaxso=  33, lnmaxso=   9.
 Exchange-correlation functional for the present dataset will be:
  LDA: Perdew-Zunger-Ceperley-Alder - ixc=2
 Citation for XC functional:
  J.P.Perdew and A.Zunger, PRB 23, 5048 (1981)

 Unit cell volume ucvol=  2.8537402E+02 bohr^3
 Angles (23,13,12)=  6.00000000E+01  6.00000000E+01  6.00000000E+01 degrees

 getcut: wavevector=  0.0000  0.0000  0.0000  ngfft=  40  40  40
         ecut(hartree)=     35.000   => boxcut(ratio)=   2.03244
- pspini: atom type   1  psp file is /home/temok/AtomicPotentialFiles/martins2abinit/31_13_Ga_TM_soc_Kocinski
- pspatm: opening atomic psp file    /home/temok/AtomicPotentialFiles/martins2abinit/31_13_Ga_TM_soc_Kocinski
-    Psp for Gallium (LDA-ca - fully relativ.) using Improved Troullier and Martins method (from JL Martins psp generator)
- 31.00000  13.00000    140720                znucl, zion, pspdat
    5    2    2    0       977   0.00000      pspcod,pspxc,lmax,lloc,mmax,r2well
    5.908271E-04    1.250000E-02              r1 and al (Hamman grid)
 Pseudopotential is in spin-orbit format
    0   0.000   0.000    1   2.0979790        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    1   0.000   0.000    2   2.0979790        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    1   0.000   0.000    2   2.0979790        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    2   0.000   0.000    2   2.0979790        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    2   0.000   0.000    2   2.0979790        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    0.00000000000000    0.00000000000000    0.00000000000000   rchrg,fchrg,qchrg
 pspatm: epsatm=   53.18789583
         --- l  ekb(1:nproj) -->
             0    0.000000
             1    1.000000    0.500000
             2   -1.000000   -1.500000
  spin-orbit 1   -2.061553    2.061553
  spin-orbit 2    2.061553   -2.061553
 pspatm: atomic psp has been read  and splines computed

- pspini: atom type   2  psp file is /home/temok/AtomicPotentialFiles/martins2abinit/33_05_As_TM_soc
- pspatm: opening atomic psp file    /home/temok/AtomicPotentialFiles/martins2abinit/33_05_As_TM_soc
- Psp for Arsenic (LDA-ca - fully relativ.) using Improved Troullier and Martins method (from JL Martins psp generator)
- 33.00000   5.00000    140720                znucl, zion, pspdat
    5    2    2    0       982   0.00000      pspcod,pspxc,lmax,lloc,mmax,r2well
    5.550194E-04    1.250000E-02              r1 and al (Hamman grid)
 Pseudopotential is in spin-orbit format
    0   0.000   0.000    1   2.0979420        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    1   0.000   0.000    2   2.0979420        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    1   0.000   0.000    2   2.0979420        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    2   0.000   0.000    2   2.0979420        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    2   0.000   0.000    2   2.0979420        l,e99.0,e99.9,nproj,rcpsp
   0.00000000   0.00000000   0.00000000   0.00000000     rms, ekb1, ekb2, epsatm
    0.00000000000000    0.00000000000000    0.00000000000000   rchrg,fchrg,qchrg
 DATA TYPE INFORMATION: 
 REAL:      Data type name: REAL(DP) 
            Kind value:      8
            Precision:      15
            Smallest nonnegligible quantity relative to 1:  0.22204460E-15
            Smallest positive number:                       0.22250739-307
            Largest representable number:                   0.17976931+309
 INTEGER:   Data type name: INTEGER(default) 
            Kind value: 4
            Bit size:   32
            Largest representable number: 2147483647
 LOGICAL:   Data type name: LOGICAL 
            Kind value: 4
 CHARACTER: Data type name: CHARACTER             Kind value: 1
  MPI-IO support is OFF

 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

 === Build Information === 
  Version       : 7.8.1
  Build target  : x86_64_linux_gnu4.4
  Build date    : 20140702

 === Compiler Suite === 
  C compiler       : gnu4.4
  C++ compiler     : gnu4.4
  Fortran compiler : gnu4.4
  CFLAGS           :  -g -O2 -mtune=native -march=native
  CXXFLAGS         :  -g -O2 -mtune=native -march=native
  FCFLAGS          :  -g -ffree-line-length-none
  FC_LDFLAGS       : 

 === Optimizations === 
  Debug level        : basic
  Optimization level : standard
  Architecture       : amd_opteron

 === Multicore === 
  Parallel build : no
  Parallel I/O   : no
  openMP support : no
  GPU support    : no

 === Connectors / Fallbacks === 
  Connectors on : yes
  Fallbacks on  : yes
  DFT flavor    : libxc-fallback+atompaw-fallback+bigdft-fallback+wannier90-fallback
  FFT flavor    : none
  LINALG flavor : netlib
  MATH flavor   : none
  TIMER flavor  : abinit
  TRIO flavor   : netcdf-fallback+etsf_io-fallback

 === Experimental features === 
  Bindings            : @enable_bindings@
  Exports             : no
  GW double-precision : no

 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++


--- !ERROR
message: |
    BIG PROBLEM WITH THE SPIN ORBIT IN PSP5NL
src_file: psp5in.F90
src_line: 430
...


 leave_new : decision taken to exit ...

--- !ERROR
message: |
    BIG PROBLEM WITH THE SPIN ORBIT IN PSP5NL
src_file: psp5in.F90
src_line: 430
...

Thanks,
Temok

emilgojny · Post by **emilgojny** » Tue Oct 07, 2014 10:14 pm

Dear temok,

What a coincidence, I've just run into the same problem. I tried to calculate the bandstructure of GaAs, with the pseudopotentials generated with the APE package i got the same "BIG PROBLEM WITH THE SPIN ORBIT IN PSP5NL".

All other similar semiconductors (GaSb, InSb etc.) work just fine, it's the arsenic that is problematic.

Did you find a solution to this problem? I've tried different cutoffs while generating pseudos but that didn't seem to work.
I always generate my own psps, and I've always been successful.

Please let me know if you (or anyone) have a solution to this problem

Thanks in advance

Emil G.

ABINIT Discussion Forums

"Big problem with the spin orbit in PSP5NL"

"Big problem with the spin orbit in PSP5NL"

Re: "Big problem with the spin orbit in PSP5NL"

Re: "Big problem with the spin orbit in PSP5NL"