#!/bin/bash


# The trajectory is a concatenated xyz-file:

cat << eor > coords.xyz
C        -1.741968   -2.220880    0.000000
Pd       -2.137509   -0.237843    0.000000
H        -2.809570   -2.499547    0.000000
H        -1.265288   -2.629932    0.895677
H        -1.265288   -2.629932   -0.895677
H        -0.755099   -0.885698    0.00000

C        -1.741968   -2.220880    0.000000
Pd       -2.137509   -0.250000    0.000000
H        -2.809570   -2.499547    0.000000
H        -1.265288   -2.629932    0.895677
H        -1.265288   -2.629932   -0.895677
H        -0.755099   -0.885698    0.00000
eor

# Input options for ADF:

cat << eor > adfinputfile
NumericalQuality=VeryGood
XC.GGA=OPBE
Basis.core=None
Charge=0
Symmetry=AUTO
RELATIVISTIC=Scalar ZORA
eprint.sfo=eig ovl
eor

# Here we call the PyFrag script.

# PyFrag conducts a total of six single point calculations, one for each fragment and the 
# whole system for each of the two geometries in the trajectory, respectively. 
# The intermediate results for each geometry are printed during the calculation. 
# At the end of its run PyFrag summarizes all results in a table which is also written 
# to an output file PyFrag<SystemName>.txt

$ADFBIN/startpython $ADFHOME/scripting/standalone/pyfrag/PyFrag.py  \
--xyzpath coords.xyz \
--fragment 2 --fragment 1 3 4 5 6 --strain 0 --strain -554.1 \
--bondlength 1 6 1.093 \
--angle 2 1 6 180 \
--VDD 2 1 6 \
--hirshfeld frag1 \
--irrepOI AA \
--irrepOI AAA \
--orbitalenergy frag1 HOMO \
--orbitalenergy frag2 LUMO \
--orbitalenergy AA frag2 5 \
--population frag1 HOMO \
--population AA frag2 5 \
--overlap frag1 HOMO frag2 LUMO \
--overlap S frag1 5 AA frag2 4 \
--adfinput basis.type=DZP \
--adfinputfile adfinputfile