Pick the right basis for your calculation¶
The development of basis sets suitable for use in relativistic calculations reflects the relative lateness of the field’s development. Because the more consistent efforts in method development started at about the mid 1980’s, it wasn’t until well into the late 1990’s that the pioneering works of the early and mid 1990’s were substantially complemented and improved upon.
The availability of basis sets has dramatically improved in recent years, notably with the work of K. G. Dyall, and Dirac users are strongly advised to use Dyall’s basis sets whenever they are available. These sets follow roughly the “correlationconsistent” philosophy introduced by Dunning and coworkers [Dunning1989], so they already contain polarization and correlation functions, but the SCF sets are designed for an adequate SCF representation rather than to match correlating sets for the valence shells.
Dyall basis sets¶
We recommend that you use the Dyall basis set repository whenever the basis sets are available for the elements of interest. In order to make that usage as convenient as possible, the following files, containing all sets currently available at the URL above (published or to be published), are made available:
quality 
valence 
corevalence 
allelectron 
DFT 

doublezeta 
dyall.v2z 
dyall.cv2z 
dyall.ae2z 
dyall.2zp 
+diffuse functions 
dyall.av2z 
dyall.acv2z 
dyall.aae2z 

triplezeta 
dyall.v3z 
dyall.cv3z 
dyall.ae3z 
dyall.3zp 
+diffuse functions 
dyall.av3z 
dyall.acv3z 
dyall.aae3z 

quadruplezeta 
dyall.v4z 
dyall.cv4z 
dyall.ae4z 
dyall.4zp 
+diffuse functions 
dyall.av4z 
dyall.acv4z 
dyall.aae4z 
The basis sets with diffuse functions are available for the s, p, and d blocks. The diffuse functions consist of one additional function in each symmetry that has functions for valence correlation.
While the division into “valence” and “corevalence” can be at times not so clearcut as for lighter elements, the option was made to stick to the usual jargon of nonrelativistic theory, particularly in relation to the “correlationconsistent” family of basis sets.
The valence basis sets are defined differently for each block, and include functions for the correlation of the following shells:
block 
shells included 

s 
ns shell, (n1)s, p shells 
p 
ns, np shells 
d 
ns, np, nd shells, (n+1)s and p shells 
f 
ns, np, nd, nf, (n+1)s, p, d shells, (n+2) s, p shells 
The choice for the f block is necessary to cover correlation of the open f shell, which becomes a semicore shell towards the end of the row. The reason for including the outer core shells for the s and d blocks is that the correlation of these shells is usually necessary for accurate results.
The corevalence basis sets include the (n2) shell for the s elements, the (n1) shell for the p elements, the (n1) shell for the d elements, and nothing extra for the f elements.
The allelectron basis sets include correlating functions for all shells, down to the 1s for all elements. These are intended for use when correlating all electrons.
The basis sets also include functions for dipole polarization of the outer core shells, as this is important for many elements. For the s block these functions are included for the outer core (n1)s and p shells; for the p block, an f function is included to polarize the (n1)d shell; for the d block, polarizing f and higher functions are included for the nd shell; none are added for the f block as the f is very compact.
The DFT basis sets do not contain the correlating functions as these are not necessary for DFT (or HartreeFock) calculations, except for the outermost shells where the functions with one unit more of angular momentum are included from the correlating sets as polarization functions. The dipole polarization functions for the outer core are also included. These basis sets are the most economical choice for DFT calculations.
You are encouraged to look in the basis set files and in the original archives published in Theor. Chem. Acc. to get a feel for what is included in each case. The archive files are available on zenodo here.
With the recent addition of Dyall basis sets for the light elements, it is no longer necessary to use the standard nonrelativistic basis sets, such as the correlationconsistent sets of Dunning and coworkers. However, because the Dyall basis sets are quite a bit larger, you might want to continue using these basis sets. It is advisable that, in order to have a balanced description when light and heavy elements are present, that one uses either contracted or uncontracted sets thoughout. For 4component calculations it is strongly recommended to use the basis sets uncontracted (which is the case for the files listed above).
See the Dyall basis set repository for the latest updates and the appropriate basis set references. In case of errors or omissions on any of the files in this directory, users are kindly asked to contact the authors of DIRAC.
Other relativistic basis sets¶
Apart from Dyall’s sets, you can choose several different basis sets based upon geometric progressions of exponents. One such set is that of K. Faegri, also available in the basis set library, but with the drawback that you may need to extend it by adding polarization functions.
Nonrelativistic and scalarrelativistic basis sets¶
The DIRAC distribution shares a large library of standard nonrelativistic and scalarrelativistic basis sets with the Dalton program. These basis sets can be found in the directory basis_dalton of the DIRAC distribution.
These basis sets are not all suitable for relativistic calculations, especially not for the heavier elements. Basis sets developed for full relativistic calculations (including spinorbit coupling) can be found in the directory basis.