/usr/share/psi4/samples/sapt8/test.in 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 | #! SAPT0(ROHF) open-shell computation of CN - Ne interaction energy
#! First with jun-cc-pVDZ and density fitted integrals with ROHF
#! Then with cc-pVDZ and direct integrals, except for dispersion
#! that is computed with cc-pVDZ-ri density fitting with ROHF.
memory 1 GB
molecule {
0 2
C 0.000000 0.000000 0.000000
N 0.000000 0.000000 2.21
--
0 1
Ne 0.000000 0.000000 -6.410228688
units bohr
symmetry c1
no_reorient
no_com
}
set {
reference rohf
scf_type df
guess sad
basis jun-cc-pVDZ
df_basis_scf aug-cc-pVDZ-jkfit
df_basis_sapt aug-cc-pVDZ-ri
}
energy('sapt0')
Eelst = psi4.get_variable("SAPT ELST ENERGY")
Eexch = psi4.get_variable("SAPT EXCH ENERGY")
Eind = psi4.get_variable("SAPT IND ENERGY")
Edisp = psi4.get_variable("SAPT DISP ENERGY")
ET = psi4.get_variable("SAPT ENERGY")
Eref = [ -0.00005818602, #TEST
0.00019907113, #TEST
-0.00002399403, #TEST
-0.00010370304, #TEST
0.00001318804 ] #TEST
compare_values(Eref[0], Eelst, 7, "SAPT0 Eelst") #TEST
compare_values(Eref[1], Eexch, 7, "SAPT0 Eexch") #TEST
compare_values(Eref[2], Eind, 7, "SAPT0 Eind") #TEST
compare_values(Eref[3], Edisp, 7, "SAPT0 Edisp") #TEST
compare_values(Eref[4], ET, 7, "SAPT0 Etotal") #TEST
set {
guess sad
basis cc-pVDZ
df_basis_sapt cc-pVDZ-ri
scf_type direct
}
energy('sapt0')
Eelst = psi4.get_variable("SAPT ELST ENERGY")
Eexch = psi4.get_variable("SAPT EXCH ENERGY")
Eind = psi4.get_variable("SAPT IND ENERGY")
Edisp = psi4.get_variable("SAPT DISP ENERGY")
ET = psi4.get_variable("SAPT ENERGY")
Eref = [ -0.00012024341, #TEST
0.00018948726, #TEST
-0.00001750737, #TEST
-0.00007559844, #TEST
-0.00002386197 ] #TEST
compare_values(Eref[0], Eelst, 7, "SAPT0 Eelst") #TEST
compare_values(Eref[1], Eexch, 7, "SAPT0 Eexch") #TEST
compare_values(Eref[2], Eind, 7, "SAPT0 Eind") #TEST
compare_values(Eref[3], Edisp, 7, "SAPT0 Edisp") #TEST
compare_values(Eref[4], ET, 7, "SAPT0 Etotal") #TEST
|