Example: Basis set and integration accuracy convergence test: Methane

Download ConvergenceTestCH4.run

In this example you will find how to use adfprep to test convergence of the bonding energy with respect to basis set and integration accuracy. adfreport is used to collect the resulting bonding energies.

rm -f runset
for b in SZ DZ DZP TZP TZ2P QZ4P
do
    "$ADFBIN/adfprep" -t "$ADFHOME/examples/adf/ConvergenceTestCH4/methane.adf" -b $b -j methane.$b >> runset
done

chmod +x runset
./runset

echo Results
echo Basis set convergence of Bonding Energy, SZ DZ DZP TZP TZ2P QZ4P
for b in SZ DZ DZP TZP TZ2P QZ4P
do
    "$ADFBIN/adfreport" "methane.$b.t21" BondingEnergy
done

rm -f runset
for i in 2 3 4 5
do
    "$ADFBIN/adfprep" -t "$ADFHOME/examples/adf/ConvergenceTestCH4/methane.adf" -b DZP -i $i -j methane.$i >> runset
done
"$ADFBIN/adfprep" -t "$ADFHOME/examples/adf/ConvergenceTestCH4/methane.adf" -b DZP -i Basic -j methane.bb >> runset
"$ADFBIN/adfprep" -t "$ADFHOME/examples/adf/ConvergenceTestCH4/methane.adf" -b DZP -i Normal -j methane.bn >> runset
"$ADFBIN/adfprep" -t "$ADFHOME/examples/adf/ConvergenceTestCH4/methane.adf" -b DZP -i Good -j methane.bg >> runset
"$ADFBIN/adfprep" -t "$ADFHOME/examples/adf/ConvergenceTestCH4/methane.adf" -b DZP -i VeryGood -j methane.bv >> runset
"$ADFBIN/adfprep" -t "$ADFHOME/examples/adf/ConvergenceTestCH4/methane.adf" -b DZP -i Excellent -j methane.be >> runset

chmod +x runset
./runset

echo Integration convergence of Bonding Energy, 2 3 4 5
for i in 2 3 4 5
do
    "$ADFBIN/adfreport" "methane.$i.t21" BondingEnergy
done

echo Integration Becke convergence of Bonding Energy, Basic, Normal, Good, VeryGood, Excellent
"$ADFBIN/adfreport" "methane.bb.t21" BondingEnergy
"$ADFBIN/adfreport" "methane.bn.t21" BondingEnergy
"$ADFBIN/adfreport" "methane.bg.t21" BondingEnergy
"$ADFBIN/adfreport" "methane.bv.t21" BondingEnergy
"$ADFBIN/adfreport" "methane.be.t21" BondingEnergy

echo Ready