Example: Asymptotically correct XC potentials: CO

Download CO_model.run

For property calculations, xc potentials with asymptotically correct (-1/r) behavior outside the molecule, the results tend to be superior to regular LDA or GGA calculations. This is especially true for small molecules and for properties that depend heavily on the proper description of the outer region of the molecule. In the example, all-electron basis sets are used. This is mandatory for the SAOP potential.

In the first example, excitation energies are calculated with the GRACLB potential. This potential requires one number as argument: the experimental ionization potential in atomic units. This number can be either based on an experimental value, or on previous GGA total energy calculations.

$ADFBIN/adf <<EOR
title CO excitations grac potential

BeckeGrid
 Quality good
End

XC
 Model GRACLB 0.515
End

Atoms
O  0             0   0
C  1.128205364   0   0
end

Excitation
 Lowest 10
 Onlysing
End

Basis
 Type TZ2P
 Core None
End

end input
EOR

rm TAPE21 logfile

The same calculation with the SAOP xc potential would differ in the XC block only:

XC
 Model SAOP
End

SAOP depends on the orbitals which makes it more expensive to evaluate than GRAC for large molecules, but is much easier to use, since it does not require an ionization potential parameter as input.