
Exchange correlation potential making use of GGA, the parametrization is given by the PBE functional
PBE exchange correlation functional
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | ir | |||
| integer, | intent(in) | :: | np | |||
| real(kind=dp), | intent(in) | :: | rv | |||
| real(kind=dp), | intent(in), | dimension(lmpotd, 2) | :: | rholm |
l,m decomposed charge density |
|
| real(kind=dp), | intent(inout), | dimension(ijd,2) | :: | vxcp |
XC-potential |
|
| real(kind=dp), | intent(inout), | dimension(ijd) | :: | excp |
XC-energy |
|
| real(kind=dp), | intent(in), | dimension(ijd, lmpotd) | :: | ylm |
real spherical harmonic to a given l,m |
|
| real(kind=dp), | intent(in), | dimension(ijd, lmpotd) | :: | dylmt1 | ||
| real(kind=dp), | intent(in), | dimension(ijd, lmpotd) | :: | dylmf1 | ||
| real(kind=dp), | intent(in), | dimension(ijd, lmpotd) | :: | dylmf2 | ||
| real(kind=dp), | intent(in), | dimension(ijd, lmpotd) | :: | dylmtf | ||
| real(kind=dp), | intent(in), | dimension(irmd, lmpotd) | :: | drrl | ||
| real(kind=dp), | intent(in), | dimension(irmd, lmpotd) | :: | ddrrl | ||
| real(kind=dp), | intent(in), | dimension(irmd, lmpotd) | :: | drrul | ||
| real(kind=dp), | intent(in), | dimension(irmd, lmpotd) | :: | ddrrul | ||
| integer, | intent(in) | :: | irmd |
Maximum number of radial points |
||
| integer, | intent(in) | :: | lmpotd |
(lpot+1)**2 |
||
| integer, | intent(in) | :: | lmmax0d |
(LMAX+1)^2 |
||
| logical, | intent(in) | :: | use_sol |
use_sol=0 -> PBE, use_sol=1 -> PBEsol |
Driver routine for PBE GGA subroutines
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp), | intent(in), | dimension(2) | :: | ro | ||
| real(kind=dp), | intent(in), | dimension(2) | :: | rol | ||
| real(kind=dp), | intent(in), | dimension(3,2) | :: | ro1 | ||
| real(kind=dp), | intent(in), | dimension(5,2) | :: | ro2 | ||
| real(kind=dp), | intent(in) | :: | xr | |||
| real(kind=dp), | intent(in) | :: | s | |||
| real(kind=dp), | intent(in) | :: | c | |||
| real(kind=dp), | intent(out) | :: | v1 |
xc-potential up |
||
| real(kind=dp), | intent(out) | :: | v2 |
xc-potential down |
||
| real(kind=dp), | intent(out) | :: | exc |
xc-energy |
||
| real(kind=dp), | intent(in) | :: | um | |||
| real(kind=dp), | intent(in) | :: | bet |
PBE exchange for a spin-unpolarized electronic system
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp), | intent(in) | :: | rho |
density |
||
| real(kind=dp), | intent(in) | :: | s |
with |
||
| real(kind=dp), | intent(in) | :: | u |
|
||
| real(kind=dp), | intent(in) | :: | v | |||
| real(kind=dp), | intent(out) | :: | ex |
exchange energy per electron |
||
| real(kind=dp), | intent(out) | :: | vx |
exchange potential |
||
| integer, | intent(in) | :: | llda | |||
| real(kind=dp), | intent(in) | :: | um |
This subroutine evaluates the correlation energy per particle for PBE xc-potential
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp), | intent(in) | :: | rs |
Wigner-Seitz radius |
||
| real(kind=dp), | intent(in) | :: | zeta |
relative spin polarization |
||
| real(kind=dp), | intent(in) | :: | t |
|
||
| real(kind=dp), | intent(in) | :: | uu |
|
||
| real(kind=dp), | intent(in) | :: | vv |
|
||
| real(kind=dp), | intent(in) | :: | ww |
|
||
| real(kind=dp), | intent(out) | :: | ec |
correlation energy per particle |
||
| real(kind=dp), | intent(out) | :: | vcup |
spin-dn correlation potential |
||
| real(kind=dp), | intent(out) | :: | vcdn |
spin-up correlation potential |
||
| integer, | intent(in) | :: | llda | |||
| real(kind=dp), | intent(in) | :: | bet |
coefficient in gradient expansion for correlation, a. |
Slimmed down version of GCOR used in PW91 routines, to interpolate LSD correlation energy
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp) | :: | a | ||||
| real(kind=dp) | :: | a1 | ||||
| real(kind=dp) | :: | b1 | ||||
| real(kind=dp) | :: | b2 | ||||
| real(kind=dp) | :: | b3 | ||||
| real(kind=dp) | :: | b4 | ||||
| real(kind=dp) | :: | rtrs | ||||
| real(kind=dp) | :: | gg | ||||
| real(kind=dp) | :: | ggrs |