/usr/share/doc/psi4/examples/pywrap-cbs1.dat is in psi4-data 1:1.1-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 | #! Various basis set extrapolation tests
molecule {
Li
}
set {
reference rohf
scf_type pk
df_scf_guess false
}
e_cbs = energy(cbs,
scf_wfn='scf',
scf_basis='aug-cc-pV[DTQ]Z',
scf_scheme=scf_xtpl_helgaker_3)
clean()
# Here's another example
molecule h2 {
H
H 1 R
R = 1
}
set {
reference rhf
guess core
mp2_type conv
df_scf_guess true
}
e_cbs = energy(cbs,
corl_wfn='mp2',
corl_basis='cc-pV[TQ]Z',
corl_scheme=driver_cbs.corl_xtpl_helgaker_2,
delta_wfn='ccsd',
delta_basis='cc-pV[DT]Z',
delta_scheme=driver_cbs.corl_xtpl_helgaker_2)
clean()
# One final example
molecule {
He
}
set {
reference rhf
df_scf_guess false
}
e_cbs = energy(cbs,
scf_basis='aug-cc-pV[DTQ]Z',
scf_scheme=driver_cbs.scf_xtpl_helgaker_3,
corl_wfn='ccsd',
corl_basis='aug-cc-pV[DT]Z',
corl_scheme=driver_cbs.corl_xtpl_helgaker_2)
clean()
# Example with default extrapolation schemes
e_cbs = energy(cbs,
corl_wfn='mp2',
corl_basis='cc-pv[dt]z',
delta_wfn='mp5',
delta_basis='cc-pVDZ')
clean()
def myownmy_2(functionname, zLO, valueLO, zHI, valueHI, verbose=True):
"""Note that the _2 is necessary to signal this is a fn involving 2 zeta-levels.
Yes, this fn is just the corl_xtpl_helgaker_2 w/o error checking and printing.
"""
if isinstance(valueLO, float):
value = (valueHI * zHI ** 3 - valueLO * zLO ** 3) / (zHI ** 3 - zLO ** 3)
beta = (valueHI - valueLO) / (zHI ** (-3) - zLO ** (-3))
return value
elif isinstance(valueLO, (core.Matrix, core.Vector)):
beta = valueHI.clone()
beta.subtract(valueLO)
beta.scale(1 / (zHI ** (-3) - zLO ** (-3)))
beta.name = 'Helgaker SCF (%s, %s) beta' % (zLO, zHI)
value = valueHI.clone()
value.scale(zHI ** 3)
tmp = valueLO.clone()
tmp.scale(zLO ** 3)
value.subtract(tmp)
value.scale(1 / (zHI ** 3 - zLO ** 3))
value.name = 'Helgaker Corr (%s, %s) data' % (zLO, zHI)
return value
else:
raise ValidationError("scf_xtpl_helgaker_2: datatype is not recognized '%s'." % type(valueLO))
h2.R = 0.747953788665
# 5a is equivalent specification to 5b
#e_cbs = energy(cbs,
# molecule=h2,
# corl_wfn='mp2',
# corl_basis='aug-cc-pv[t,Q]z',
# delta_scheme=myownmy_2,
# delta_wfn='ccsd(t)',
# delta_basis='cc-pV[dt]Z')
#compare_values(-1.17505568495, e_cbs, 6, "[5a] user-def xtpl cbs()")
e_cbs = energy('MP2/aug-cc-pV[T,Q]Z + D:CCSD(T)/cc-pv[dt]z', delta_scheme=myownmy_2, molecule=h2)
|