1 if { [ catch {exec which qchem } ] } then {

2     set msg "Unsupported test: Requires Q-Chem set up in path)"

3     puts stdout $msg

4     return $msg

5 }

6 #

7 # Default (STO-3G, HF calculation)

8 #

9 z_create zmatrix=z {

10 zmatrix angstrom

11 o

12 h 1 1

13 h 1 1 2 100

14 }

15 z_to_c zmatrix=z coords=c

16

17 matrix e

18 energy coords=c theory=qchem : { basis=VDZ** } energy=e

19 set final_energy  [ get_matrix_element matrix=e indices= {0 0 } ]

20 delete_object e

21

22 # Test result

23 # set etest  -74.96565701

24

25 # Replace with SV value until matching sto3g basis can be fixed up

26 set etest  -75.91167156

27

28 if { abs( $final_energy - $etest) > 1.0e-6 } {

29    chemerr "Failed .. E= $final_energy, should be $etest"

30 } else {

31    return "OK"

32 }