Example: Asymptotically correct XC potentials: CO

Download CO_model.run

#!/bin/sh

# 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
  Excitations
   Lowest 10
   Onlysing
  End
  Basis
   Type TZ2P
   Core None
  End
EOR

rm TAPE21 logfile


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


$ADFBIN/adf <<EOR
  title CO excitations SAOP 
  BeckeGrid
   Quality good
  End
  XC
   Model SAOP
  End
  Atoms
    O  0             0   0
    C  1.128205364   0   0
  end
  Excitations
   Lowest 10
   Onlysing
  End
  Basis
   Type TZ2P
   Core None
  End
EOR

mv TAPE21 CO.t21


# 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.