Example: Compound Fragments: Ni(CO)4

Download Frags_NiCO4.run

#! /bin/sh

# An illustration of the fragment feature of ADF.

# A transition metal complex is built from a Nickel atom and four CO fragments.
# The outcomes allows for an analysis (of molecular orbitals and the Bonding
# energy) in terms of the fragment properties. It is a Single Point calculation.
# Geometry optimization would not have been possible in this set-up because an
# optimization requires that only single-atom fragments are used.

# The three atoms are created first: C, O, and Ni. For Carbon and Oxygen a type-
# DZ basis set is used (double-zeta) using the Basis key, while Ni gets a type-
# TZP basis (triple-zeta plus polarization).


# == CO ==


# The CO molecule, to serve as a fragment template in Ni(CO)4 , is computed from
# the atomic fragments C and O. The coordinate values (atoms) are in bohr,
# rather than in Angstrom because the unit-of-length is redefined by the key
# units with subkey length.

# The key scf is used to specify a somewhat tighter convergence criterion than
# the default, just to illustrate how to do this (normal settings are quite
# adequate).

# The TAPE21 result file is renamed t21.CO.


$ADFBIN/adf -n1 <<eor
  create Ni  $ADFRESOURCES/TZP/Ni.2p
eor
mv TAPE21 t21.Ni
rm logfile


$ADFBIN/adf -n1 <<eor
  title CO (as fragment for NiCO4)
  SCF
    converge  1e-8
  end
  EPRINT
    SFO eig ovl
  END
  units
    length  bohr
  end
  atoms
    C   0  0  0
    O   0  0  2.15617844
  end
  Basis
    Type DZ
    Core Small
  End
eor
mv TAPE21 CO.t21
rm logfile


# One needs to include the subkey SFO of the key EPRINT with arguments eig and
# ovl in order to get the SFO MO coefficients and SFO overlap matrix printed on
# standard output.


# == Main calculation ==


# Apart from the title, the input contains comment. This does not specify
# computational parameters but is only echoed in the output header, similar to
# the title. Contrary to the title, however, such comments are not preserved,
# apart from their echo in output and they are not written to TAPE21 or any
# other result file.

# The atomic coordinates (atoms) are given in bohr (Units). To supply the
# numerical values use is made of user-defined constants (define): 2.0053211 and
# 3.2501913. This is convenient and it prevents typing errors when several
# coordinate values are identical, in particular when they carry a lot of
# decimal places.

# The Atoms records contain also a specification of the fragments to which the
# respective atoms belong: four different CO fragments. No fragment is specified
# for the Ni atom, which implies that it is a fragment on its own.

# The numbers at the very left of the records (1 through 9, with (optionally) a
# period after them), have no relevance. You can set them for ease of reference
# or counting.


$ADFBIN/adf <<eor
  title Ni(CO)4,  from fragments Ni and CO
  COMMENT
    No geometry optimization possible, because not all fragments
    are single atoms
  END
  units
    length  bohr
  end
  EPRINT
    SFO eig ovl
  END
  atoms
    1. Ni   0       0        0  
    2. C    2.0053211    2.0053211     2.0053211     f=CO/1
    3. C   -2.0053211   -2.0053211     2.0053211     f=CO/2
    4. C    2.0053211   -2.0053211    -2.0053211     f=CO/3
    5. C   -2.0053211    2.0053211    -2.0053211     f=CO/4
    6. O    3.2501913   3.2501913    3.2501913     f=CO/1
    7. O   -3.2501913  -3.2501913    3.2501913     f=CO/2
    8. O    3.2501913  -3.2501913   -3.2501913     f=CO/3
    9. O   -3.2501913   3.2501913   -3.2501913     f=CO/4
  end
  fragments
    CO  CO.t21
    Ni  t21.Ni
  end
eor
mv TAPE21 NiCO4.t21
rm logfile