Example: NMR Spin-spin coupling constants, finite nucleus: PbH4

Download PbH4_finitenuc.run

#!/bin/sh

# Example for a finite nucleus calculation and the calculation of NMR spin-spin
# coupling constants.

# One of the quality-determining factors for the calculation of NMR coupling
# constants is the chosen basis set, especially one needs enough tight s
# functions. If one has a large enough basis set in the core region one can see
# an effect of using a finite size of the nucleus instead of a point nucleus,
# especially for heavy nuclei. Such large basis sets can be found for some
# elements in $ADFRESOURCES/ZORA/jcpl, which are basis sets especially designed
# for NMR spin-spin coupling calculations. In this example first a basis set for
# Pb is made which has many tight s functions.


cat <<eor > $PWD/Pb
Pb Basis TZ2P - all steep, plus high exponent fcts. added: iteration 3

BASIS
 1s    39358.2303582207
 1s    20080.7297746024
 1s    10245.2702931645
 1s    5227.17872100229
 1s    2666.92791887872
 1s    1360.67750963200
 1s    694.223219200000
 1s    354.195520000000
 1s    180.712000000000
 2p    39358.2303582207
 2p    20080.7297746024
 2p    10245.2702931645
 2p    5227.17872100229
 2p    2666.92791887872
 2p    1360.67750963200
 2p    694.223219200000
 2p    354.195520000000
 2p    180.712000000000
        1S   92.200
        2S  106.000
        2S   40.200
        3S   51.150
        3S   21.850
        4S   13.400
        4S    9.550
        5S    7.650
        5S    5.200
        6S    3.700
        6S    2.400
        6S    1.580
        2P  196.000
        2P   70.200
        2P   38.300
        3P   21.750
        3P   15.900
        4P   14.500
        4P    9.900
        5P    6.600
        5P    4.200
        6P    2.700
        6P    1.600
        6P     .960
        3D   29.400
        3D   18.800
        4D   13.000
        4D    8.450
        5D    6.000
        5D    3.650
        5D    2.200
        4F   17.560
        4F    9.950
        4F    5.750
        6D    1.600
        5F    2.500
END

CORE   0 0 0 0
END

DESCRIPTION

END

FIT
 1S   78716.46
 1S   47664.51
 1S   28861.89
 1S   17476.49
 1S   10582.39
 1S    6407.87
 1S    3880.10
 1S    2349.49
 1S    1422.66
 1S     861.45
 1S     521.63
 1S     315.86
 1S     191.26
 1S     115.81
 1S      70.13
 2S      78.58
 2S      53.32
 3S      53.08
 3S      38.26
 4S      36.45
 4S      27.31
 5S      25.48
 5S      19.63
 6S      18.11
 6S      14.24
 7S      13.06
 7S      10.45
 8S       9.55
 8S       7.74
 9S       7.06
 9S       5.79
10S       5.28
10S       4.37
10S       3.62
11S       3.30
11S       2.75
11S       2.30
11S       1.92
 2P 130.50
 3P 105.00
 4P  83.37
 5P  66.14
 6P  52.66
 6P  35.96
 7P  28.77
 7P  20.17
 8P  16.23
 8P  11.63
 9P   9.41
 9P   6.86
10P   5.58
10P   4.14
11P   3.38
11P   2.54
 3D 111.00
 4D  81.70
 5D  60.58
 6D  45.35
 7D  34.28
 7D  22.74
 8D  17.35
 8D  11.80
 9D   9.08
 9D   6.31
10D   4.89
10D   3.46
 4F 102.15
 5F  63.70
 6F  40.72
 7F  26.59
 8F  17.70
 9F  11.98
 9F   7.34
10F   5.04
10F   3.16
 5G  48.25
 6G  32.57
 7G  22.37
 8G  15.60
 8G   9.73
 9G   6.88
 9G   4.40
END

eor


# This large basis set for Pb is used in ADF calculations on PbH4 and the
# calculation of the NMR spin-spin coupling constants.


$ADFBIN/adf <<eor
UNITS
    length Angstrom
    angle Degree
END

NuclearModel Gaussian
Print Nuclei

ATOMS
 Pb      0.000000      0.000000      0.000000
  H     -1.023703      1.023703      1.023703
  H      1.023703     -1.023703      1.023703
  H     -1.023703     -1.023703     -1.023703
  H      1.023703      1.023703     -1.023703
END

BASIS
type DZP
core None
Pb $PWD/Pb
END

XC
 gga PBE
END
usespcode

SCF
iterations 100
converge 1.0e-8 1.0e-6
END

BeckeGrid
 quality verygood
End

relativistic scalar zora

SAVE TAPE10

eor

## end scalar run, now do the coupling constant

$ADFBIN/cpl <<eor
gga
Print Nuclei
nmrcoupling
scf converge=1e-5 iterations=25
 nuclei 1 2
end
eor