#! /bin/sh

# Calculation of the excitation energies of the complex TCNE-Benzene using the
# simplified Tamm-Dancoff (sTDA) method by Grimme, for a range-separated
# function. This method is meant for hybrid functionals or range-separated
# functionals during the SCF. In this example the range-separated functional 
# wB97 is used during the SCF. In the the calculation of the excitation
# energies the sTDA method is used in which the required integrals are
# approximated. For range-separated functional one needs to set the parameters
# manually. Note that the alpha and beta parameters for sTDA for range separated
# functionals in the Risthaus paper (dx.doi.org/10.1039/C3CP54517B) are mixed
# up. A TZP or TZ2P basis set is recommended for this method.

# First the atoms are calculated. To calculate this range-separated functional
# wB97 the RIHARTREEFOCK method is required, and LibXC is needed. For speed
# reasons the atoms are calculated separately, such that these calculations
# might run in parallel. If the BASIS key is used the atoms are calculated
# serially.

for n in H C N
do
$ADFBIN/adf <<eor
Create $n file=$ADFRESOURCES/DZ/$n
XC
 LibXC wB97
End
eor
mv TAPE21 $n.t21
rm logfile
done

# Next the excitation energies of the complex TCNE-Benzene are calculated, in
# which charge transfer (CT) states are important. Symmetry NOSYM is required.
# The subkey DESCRIPTORS is added to the key EXCITATIONS such that CT descriptors 
# are calculated. The the atomic distance criterion used for the calculation of 
# CT_AT was set to 3.0 Angstrom.
# The sTDA parameters should be set with the subkeys GRIMMEAEX, GRIMMEALPHA, and
# GRIMMEBETA of the key MODIFYEXCITATION.

$ADFBIN/adf <<eor
Atoms
C     -2.274205   -1.398789   -0.000000
C     -2.274206    1.398791   -0.000000
C     -2.274506   -0.699551    1.213150
C     -2.274505   -0.699551   -1.213150
C     -2.274506    0.699551    1.213150
C     -2.274505    0.699551   -1.213150
H     -2.278385   -2.491323   -0.000000
H     -2.278386    2.491321   -0.000000
C      1.368996   -0.686302    0.000000
C      1.368996    0.686301    0.000000
C      1.379968   -1.434052   -1.220900
C      1.379968   -1.434052    1.220895
C      1.379968    1.434051   -1.220900
C      1.379967    1.434051    1.220895
N      1.400955   -2.054479   -2.204820
N      1.400955   -2.054479    2.204824
N      1.400955    2.054481   -2.204820
N      1.400954    2.054481    2.204824
H     -2.276046   -1.246672    2.158733
H     -2.276045   -1.246672   -2.158730
H     -2.276045    1.246671   -2.158730
H     -2.276046    1.246671    2.158733
End
Fragments
 H H.t21
 C C.t21
 N N.t21
End
symmetry nosym
excitations
 stda
 onlysing
 descriptors
 descriptors_CT_AT_Rab 3.0
end
XC
 LibXC wB97
END
modifyexcitation
 GRIMMEPERTC
 GRIMMEAEX   0.61
 GRIMMEALPHA 4.41
 GRIMMEBETA  8.00
end
eor