|
|
|
Use numerical orbitals (from the DIRAC internal subprogram) and add STOs for increased flexibility. This turns out to be more (although not by much) rapidly converging to the basis set limit than when you use STOs only, in particular for relatively small basis sets.