Explore Abinit and around
- ABINIT Coding
- Doc for new developers
- Organizing events
The format 6 for ABINIT pseudopotentials allows to use pseudopotentials generated from the FHI98PP code (Fritz-Haber-Institute, Berlin). This code is available at URL http://www.FHI-Berlin.MPG.DE/th/fhi98md/fhi98PP .
A few lines must be added to the file generated by the FHI98PP code, and are described in the present file. ABINITv1.9 is able to read format 6 pseudopotential files without core correction. ABINITv2.2 is able to read them with core correction (thanks to AF).
1) The original cpi file
We will suppose that the user has been able to generate a 'cpi' filen from the FHI98PP code. There is one such file in the ~abinit/tests/Psps_for_tests directory, with the name al_h.cpi (the FHI98PP code generated a file named “al:h.cpi”, but the name was changed, to allow portability under DOS/Windows) . It begins by:
0.30000000000000E+01 3 0.0000 0.0000 0.0000 0.0000 0.0000 .00e+00 .00e+00 0.0000 .00e+00 .00e+00 0.0000 .00e+00 .00e+00 0.0000 .00e+00 .00e+00 0.0000 .00e+00 .00e+00 0.0000 .00e+00 .00e+00 0.0000 .00e+00 .00e+00 0.0000 .00e+00 .00e+00 0.0000 .00e+00 .00e+00 493 0.10247000000000E+01 1 0.48076923076923E-03 0.91926957204792E-04 0.72762392428392E+00 2 0.49264423076923E-03 0.94197553047751E-04 0.72762434262851E+00 3 0.50481254326923E-03 0.96524232608030E-04 0.72762477651319E+00 ... ...
From these numbers, we will make direct use of:
2) The new file
For the purpose of generating the corresponding file readable by ABINIT (also found in the ~abinit/tests/Psps_for_tests directory, with the name 13al.981214.fhi), the following seven lines have been added at the beginning of the al_h.cpi file:
Aluminum, fhi98PP : Hamann-type, LDA CA PerdewWang, l=2 local 13.000 3.000 981214 zatom,zion,pspdat 6 7 2 2 493 0 pspcod,pspxc,lmax,lloc,mmax,r2well 0.0 0.0 0.0 rchrg, fchrg, qchrg 5--- These two lines are available for giving more information, later 6 7-Here follows the cpi file from the fhi98pp code-
Similar lines must be added at the beginning of other cpi files, in order to make them readable by ABINIT.
Line 1 is simply a header, that might include any information, and that will be printed without modification in the output of ABINIT.
Line 2 describes:
The two first information are crucial, the third one is not really important. The atomic number is 1 for Hydrogen, 8 for Oxygen, and so on. The ionic charge should be the same number as the first number mentioned in the cpi file (called 'number of valence electrons').
Line 3 describes:
Line 4 describes 3 parameters needed to read the model core-charge if core correction is needed.
rchrg can be inferred by the pp file. Indeed, if the pp has been generated with the core-correction option, at the end of 'cpi' file there will be:
do i=1,mmax rad(i), f(i), f1(i), f2(i) end do
that is, the radial grid, the core charge (f), and their first(f1) and second (f2) derivative on the radial mesh; rchrg can be equal to any value of rad(i) to which f,f1 and f2 are reasonably vanished.
Line 5-7 have no interest up to now.
A pseudopotential with format 6 will be treated by the routine psp6in.f, that calls psp5lo.f (local part) and psp5nl.f (non-local part). There is no psp6lo.f or psp6nl.f .
The integral of (V®+Zion/r) r^2 in psp5lo.f is performed from 0 to the highest allowed radius (usually about 100 a.u.), except that beyond 20 a.u. no value of abs(V®) larger than 2.0d-8 is tolerated. This will allow to cut off spurious behaviour of pseudopotential whose data file is written in single precision.