#! /bin/sh
# This sample illustrates the DOS property program to compute density-of-states
# data, for energy-dependent analysis.
# First, the Cu_4 CO molecule is calculated (ADF), using single-atom fragments.
$AMSBIN/ams <<eor
System
atoms [Bohr]
Cu 2.784 0.0 0.0
Cu -1.392 2.411014724135877 0.0
Cu -1.392 -2.411014724135877 0.0
Cu 0.0 0.0 -3.9371705576466964
C 0.0 0.0 2.65
O 0.0 0.0 4.91
end
end
Task SinglePoint
Engine ADF
title Cu4CO (3,1) from atoms
basis
core small
type TZP
end
print SFO
xc
gga PostSCF Becke Perdew
end
EndEngine
eor
# The PostSCF feature in the specification of the XC functional is used: the
# 'Becke-Perdew' GGA corrections are not included self-consistently but applied
# to the energy evaluation after the self-consistent LDA solution has been
# obtained.
# The utility program dos requires as input file a file named TAPE21 in the current directory,
# or if one uses the key ADFFILE a (relative) path to an adf.rkf file.
$AMSBIN/dos <<eor
adffile ams.results/adf.rkf
file dostxt
energyrange npoint=36 e-start=-25 e-end=10
tdos
! Cu 3d partial DOS
gpdos
a1 14:22 32:34
a2 5:10
e1:1 18:32 37:42
e1:2 18:32 37:42
end
! The same but using BAS
gpdos
BAS 17:34 57:74 97:114 137:154
end
! The same as above, but using much less complicated input
gpdos
ATYPE Cu d
end
! Overlap PDOS Cu 3d and CO 2p
opdos
ATYPE Cu 3d
SUBEND
ATOM 5 2p
ATOM 6 2p
end
eor
echo "DOS result file"
echo "==============="
cat dostxt
echo "END of DOS result file"
# Here, the total density of states, as well as various partial densities of
# states, are computed. You may feed the results found in the dostxt file into a
# plotting program such as gnuplot. The result is not displayed here. See the
# ADF manual for more detailed info about the dos program.