#! /bin/sh
# General dispersion coefficients (beyond de dipole-dipole C6 interaction
# coefficient) are computed with the auxiliary program DISPER. It uses two
# output files from previous ADF Response calculations. In the example, the two
# ADF runs are one and the same and the relevant TENSOR output file is used
# twice.
$ADFBIN/adf <<EOR
title Van der Waals coefficients HF
atoms
H 0 0 -0.8708056087
F 0 0 0.04619439132
end
Basis
Type DZP
Core Small
End
symmetry C(8v)
RESPONSE
MAXWAALS 8
VANDERWAALS 7
ALLTENSOR
ALLCOMPONENTS
END
EOR
mv TAPE21 HF.t21
rm logfile
# Polarizabilities are computed at 7 (imaginary) frequencies between 0 and
# infinity. The program determines internally the actual frequency values in
# this range to use. The user only specifies the number of them, thereby
# determining the precision of, in fact, a numerical integration over the zero-
# infinity frequency range. A value of 7 is rather low.
# MaxWaals determines that not only the C6 but also C7 and C8 coefficients are
# computed. A value higher than 8 would not be recommended, because the
# available basis sets would be inadequate for higher coefficients.
# In DISPER calculations the preparatory Response calculation must use the
# AllTensor and AllComponents subkeys.
# The calculation produces a file TENSOR. The subsequent DISPER run uses two
# such files. In this example, both are taken from the same ADF run, copying the
# TENSOR file to, respectively, tensorA and tensorB. These names are prescribed
# for a DISPER calculation.
cp TENSOR tensorA
cp TENSOR tensorB
$ADFBIN/disper -n1 <<eor
eor
# The DISPER program needs no other input than just the files tensorA and
# tensorB, which must both be present as local files. Note the '-n1' flag: this
# enforces that a single-node (non-parallel) run is performed. The current
# implementation does not support parallelization of DISPER, because the kid
# processes may not have the (local to the master!) files tensorA and tensorB.