#!/bin/sh

# A calculation of NMR nuclear spin-spin coupling constants (NSCCs) for the
# hybrid PBE0.

# The hybrid PBE0 is chosen as exchange-correlation potential in the ADF
# calculation. The key 'usespcode' is required for consistency reasons of the
# PBE0 implementation in ADF and the kernel that is used in the 'CPL' program,
# that calculates NMR spin-spin coupling constants. Symmetry should be NOSYM.
# The basis sets used are specially optimized all-electron basis sets for NMR
# spin-spin coupling calculations (in the directory
# $ADFHOME/atomicdata/ZORA/jcpl), which have extra tight functions, compared to
# a default ADF basis set. The integration accuracy is extra high (Quality
# VeryGood).


$ADFBIN/adf <<eor
   ! experimental bond length
   ATOMS
      F 0.0000 0.0000 0.0000
      H 0.0000 0.0000 0.9170
   END

   BASIS
      Type ZORA/jcpl
      Core None
   END

   usespcode

   XC
      hybrid PBE0
   END

   SYMMETRY nosym 

   BeckeGrid
      Quality verygood
   End

   scf
      converge 1e-8 1e-7
   End

eor


# The first call to cpl is as follows.
# The key 'gga' is included to use the first-order GGA potential instead of the
# first-order VWN potential. The Hartree-Fock part of the kernel is included
# automatically if a hybrid potential is used in the ADF calculation.


$ADFBIN/cpl <<eor
   gga
   nmrcoupling
   dso
   pso
   scf converge=1e-6 iterations=20 
      nuclei 1 2 
   End
eor


# The second CPL run also includes the spin-dipole (SD) term, through the SD
# subkey, which is much more time-consuming.


$ADFBIN/cpl <<eor
   gga
   nmrcoupling
   dso
   pso
   sd
   scf converge=1e-6 iterations=20
      nuclei 1 2
   End
eor