File: input.dat

package info (click to toggle)
psi4 1%3A1.2.1-2
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 229,808 kB
  • sloc: cpp: 416,201; python: 207,680; perl: 4,328; makefile: 384; sh: 158; csh: 6
file content (40 lines) | stat: -rw-r--r-- 1,078 bytes parent folder | download
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
#! DFT Functional Smoke Test

test_dict = {                         #TEST
          "B2PLYP":   -76.2924549107, #TEST
       "CAM-B3LYP":   -76.3552628688, #TEST
        "DLDF+D09":   -76.8410958136, #TEST
 "DSD-PBEP86-D3BJ":   -76.2032218018, #TEST
             "PBE":   -76.2994541311, #TEST
}                                     #TEST

molecule h2o {
0 1
O
H 1 1.0
H 1 1.0 2 104.5
}


set basis 6-31g
set dft_spherical_points 230
set dft_radial_points 75
set scf_type df

for k, v in test_dict.items():
    if k == 'DSD-PBEP86-D3BJ' and 'dftd3' not in psi4.addons():
        continue
    ret = energy(k)
    compare_values(v, ret, 5, "RKS 0 1  %15s Energy" % k)

pbe = {
    "name": "PBE-TEST",
    "x_functionals": { "GGA_X_PBE": {} },
    "c_functionals": { "GGA_C_PBE": {} },
}

ret = energy("scf", dft_functional="PBE")
compare_values(test_dict["PBE"], ret, 5, "RKS 0 1              PBE Energy, dft_functional='PBE'")  #TEST

ret = energy("scf", dft_functional=pbe)
compare_values(test_dict["PBE"], ret, 5, "RKS 0 1              PBE Energy, dft_functional=dict()") #TEST