/usr/share/psi4/samples/sapt1/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 | #! SAPT0 cc-pVDZ computation of the ethene-ethyne interaction energy, using the cc-pVDZ-JKFIT RI basis for SCF
#! and cc-pVDZ-RI for SAPT. Monomer geometries are specified using Cartesian coordinates.
Eref = [ 85.189064196429101, -0.00359915058, 0.00362911158, #TEST
-0.00083137117, -0.00150542374, -0.00230683391 ] #TEST
molecule ethene_ethyne {
0 1
C 0.000000 -0.667578 -2.124659
C 0.000000 0.667578 -2.124659
H 0.923621 -1.232253 -2.126185
H -0.923621 -1.232253 -2.126185
H -0.923621 1.232253 -2.126185
H 0.923621 1.232253 -2.126185
--
0 1
C 0.000000 0.000000 2.900503
C 0.000000 0.000000 1.693240
H 0.000000 0.000000 0.627352
H 0.000000 0.000000 3.963929
units angstrom
}
# this molecule will crash test if molecule passing broken
molecule barrier {
0 1
He
}
set {
basis cc-pvdz
guess sad
scf_type df
sad_print 2
d_convergence 11
puream true
print 1
}
energy('sapt0', molecule=ethene_ethyne)
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("SAPT0 TOTAL ENERGY")
compare_values(Eref[0], ethene_ethyne.nuclear_repulsion_energy(), 9, "Nuclear Repulsion Energy") #TEST
compare_values(Eref[1], Eelst, 6, "SAPT0 Eelst") #TEST
compare_values(Eref[2], Eexch, 6, "SAPT0 Eexch") #TEST
compare_values(Eref[3], Eind, 6, "SAPT0 Eind") #TEST
compare_values(Eref[4], Edisp, 6, "SAPT0 Edisp") #TEST
compare_values(Eref[5], ET, 6, "SAPT0 Etotal") #TEST
|