Group  1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  
1A  2A  3B  4B  5B  6B  7B  8B  1B  2B  3A  4A  5A  6A  7A  8A  
Period  
1  1 H  2 He 

2  3 Li  4 Be  5 B  6 C  7 N  8 O  9 F  10 Ne 

3  11 Na  12 Mg  13 Al  14 Si  15 P  16 S  17 Cl  18 Ar 

4  19 K  20 Ca  21 Sc  22 Ti  23 V  24 Cr  25 Mn  26 Fe  27 Co  28 Ni  29 Cu  30 Zn  31 Ga  32 Ge  33 As  34 Se  35 Br  36 Kr 

5  37 Rb  38 Sr  39 Y  40 Zr  41 Nb  42 Mo  43 Tc  44 Ru  45 Rh  46 Pd  47 Ag  48 Cd  49 In  50 Sn  51 Sb  52 Te  53 I  54 Xe 

6  55 Cs  56 Ba  71 Lu  72 Hf  73 Ta  74 W  75 Re  76 Os  77 Ir  78 Pt  79 Au  80 Hg  81 Tl  82 Pb  83 Bi  84 Po  85 At  86 Rn 

7  87 Fr  88 Ra  103 Lr  104 Rf  105 Db  106 Sg  107 Bh  108 Hs  109 Mt  110 UUn  111 UUu  112 UUb  113 Uut  114 Uuq  115 Uup  116 Uuh  117 Uus  118 Uuo 

lanthanides  57 La  58 Ce  59 Pr  60 Nd  61 Pm  62 Sm  63 Eu  64 Gd  65 Tb  66 Dy  67 Ho  68 Er  69 Tm  70 Yb  
actinides  89 Ac  90 Th  91 Pa  92 U  93 Np  94 Pu  95 Am  96 Cm  97 Bk  98 Cf  99 Es  100 Fm  101 Md  102 No 
The periodic table above gives HartreeFock and DiracFock Average Relativistic Effective (AREP) pseudopotentials (effective core potentials) for the atoms H to Ba and Lu to Hg. These pseudopotentials are finite at the origin, which makes them particularly suitable for use in QMC calculations, although they are also suitable for many other methods, including quantum chemistry calculations. These pseudopotentials were generated by John Trail and, if you use them, please cite this website and both of the following papers:
1. “Smooth relativistic HartreeFock pseudopotentials for H to Ba and Lu to Hg”
J.R. Trail and R.J. Needs, J. Chem. Phys. 122, 174109 (2005)[link]
2. “Normconserving HartreeFock pseudopotentials and their asymptotic behaviour”
J.R. Trail and R.J. Needs, J. Chem. Phys. 122, 014112 (2005)[link]
A reasonable compact reference would be:
J.R. Trail and R.J. Needs, J. Chem. Phys. 122, 174109 (2005); J.R. Trail and R.J. Needs,
J. Chem. Phys. 122, 014112 (2005), see also www.vallico.net/casinoqmc/pplib/.
More details
There are currently three types of pseudopotential available for each of the atoms H to Ba and Lu to Hg: a DiracFock Average Relativistic Effective Potential (AREP) with small core radii^{*}, a HartreeFock pseudopotential generated with small core radii, and an AREP generated with larger core radii. The pseudopotentials are designed for use with methods in which the nonrelativistic Schrödinger equation is solved. Solving the Schrödinger equation with DiracFock AREP pseudopotentials will result in the inclusion of scalar relativistic effects (spinorbit potentials are also included in the table, and when these are included the calculations will also include spinorbit effects). For most purposes DiracFock AREPs are to be preferred because they contain important relativistic effects, but we also provide HartreeFock pseudopotentials, which may be useful in some circumstances.
The pseudopotentials have s, p, and d angular momentum channels. The pseudopotentials are finite at the origin, which is very important for QMC applications and may also be advantageous in other methods. We envisage the small core DiracFock and HartreeFock pseudopotentials being used with localized basis sets such as Gaussian functions. These pseudopotentials are given tabulated on a grid, and as fits to Gaussian basis sets for use with various quantum chemistry packages. The large core (or “softer”) DiracFock pseudopotentials are designed for use with planewave basis sets. The larger core radii improve convergence with the size of the planewave basis, but the region over which the pseudopotential is nonlocal is then slightly larger, so they are more costly to use within QMC calculations.
The table also gives plots of each pseudopotential (for Gaussian fits the original tabulated representation is also plotted in the same figure for comparison), a summary of the properties of each pseudopotential, and atomic orbitals for some typical configurations. Each of the pseudopotentials has been tested in atomic calculations.
^{* }The “core radii” are the radii outside of which we demand that the atomic ground state valence orbitals resulting from allelectron and pseudopotential calculations agree.
Description of Gaussian fits
The method of fitting the pseudopotentials to a Gaussian basis set is described in Ref. 1 above [link]. Fits are reported in formats suitable for the CRYSTAL and GAUSSIAN codes, and for the GAMESS code. It is important to realise that the CRYSTAL and GAUSSIAN fits are identical, and only the format of the files differs. Unfortunately not all quantum chemistry codes can deal with the large powers of r used to multiply the Gaussian functions. GAMESS is one of the codes which does not currently support all of the powers we would like to use and we have had to fit the pseudopotentials with a more restricted basis set for use in GAMESS. This leads to some loss of quality and the CRYSTAL/GAUSSIAN fits should be used whenever possible.
Spin orbit potentials are also provided, in the standard Gaussian parameterization.
For some of the atoms parameters for core polarization potentials are also listed, from
“Manybody corevalence partitioning” E.L. Shirley and R.M. Martin, Phys. Rev. B 47, 15413 (1993).
Description of the summary file
The core radii are the radii outside of which we demand that the pseudoorbitals and allelectron orbitals agree for the initial pseudopotential. Pseudopotentials generated by inversion of the Schrödinger or DiracFock equations have long ranged nonlocal tails which must be cutoff (as described in Ref. 1 [link] and Ref. 2[link] above). The cutoff procedure forces the pseudopotentials to be local outside of the localisation radius r_{loc} whilst correcting the pseudopotentials inside of r_{loc} to conserve ground state properties^{*}. This pseudopotential is the “tabulated” pseudopotential given in the periodic table. Parameterising such a pseudopotential in terms of a Gaussian basis set normally increases r_{loc}.
The allelectron and pseudoHF LScoupling eigenvalues should agree to high precision for a HartreeFock pseudopotential, but not necessarily for an AREP pseudopotential that is generated within DiracFock theory.
The terms in the energy for a pseudoHF calculation of the ground state of the atom are given. These are useful for checking purposes, for example these energies should be reproduced by a variational Monte Carlo calculation using the pseudopotential and the atomic orbitals resulting from a HartreeFock calculation using the same
pseudopotential.
Data for a few excitation energies are also included in the summary. It should be noted that a pseudoHF calculation with an AREP pseudopotential will not reproduce the same relativistic effects as those present in the original allelectron DF atom, even if scalar relativistic corrections and spinorbit coupling are included. This implies that excitation energies for the pseudoHF atom would differ from those of the allelectron DF atom even if the pseudopotential were able to reproduce the influence of the core electrons exactly.
Please note that there is a separate summary file for each representation of each type of pseudopotential. These give the energies resulting from pseudoHF calculations carried out using the tabulated representation, the CRYSTAL/GAUSSIAN Gaussian expansion, and the GAMESS Gaussian expansion of the pseudopotential. Small differences between the energies resulting from the tabulated and parameterised representations reflect the accuracy of the Gaussian fits.
In addition please note that the atomic orbitals provided in the periodic table have been calculated from the tabulated pseudopotentials, not the parameterised pseudopotentials.
^{*} Without this procedure V_{l}(r) – V_{l’}(r) remains finite as r→∞ for some l and l’, and for all atoms except a few special cases.
Tests
Results for atomic ground and excited states obtained using these pseudopotentials are also available, together with an error analysis:
¤ Ground state eigenvalues resulting from DiracFock AE calculations (spin averaged for each l) compared with HartreeFock calculations using the AREP pseudopotentials. These calculations are not strictly speaking HartreeFock, but are the result of taking the c→∞ limit of DiracFock theory (see Ref. 1 above [link]).
¤ Excitation energies resulting from AE HartreeFock calculations (with scalar relativistic corrections), compared with those resulting from HartreeFock calculations with AREP pseudopotentials.
¤ Ground state eigenvalues resulting from AE HartreeFock calculations, compared with those resulting from HartreeFock calculations with HF pseudopotentials.
¤ Excitation energies resulting from AE HartreeFock calculations, compared with those resulting from HartreeFock calculations with HF pseudopotentials.
In the data given above the “parameterised pseudopotential” refers to the GAUSSIAN/CRYSTAL Gaussian fit for each pseudopotential. Equivalent analyses of the data resulting from the GAMESS parameterisations are given in:
¤ Ground state eigenvalues from AE DiracFock calculations, compared with those resulting from HartreeFock calculations using the AREP pseudopotentials.
¤ Excitation energies from AE HartreeFock calculations, compared with those resulting from HartreeFock calculations using the AREP pseudopotentials.
¤ Ground state eigenvalues from AE HartreeFock calculations, compared with those resulting from HartreeFock calculations using the HF pseudopotentials.
¤ Excitation energies from AE HartreeFock calculations, compared with those resulting from HartreeFock calculations using the HF pseudopotential.
QMC notes
The pseudopotentials fitted to Gaussians are meant for use in quantum chemistry codes, which are often used to generate orbitals for QMC trial wave functions. The tabulated pseudopotentials can be used in other codes, such as planewave DFT codes. CASINO uses tabulated pseudopotentials, and the periodic table provides input files in the correct format. Note that the tabulated pseudopotentials are those which were generated on radial grids which have not been fitted to Gaussians. We recommend that these tabulated pseudopotentials be used in QMC calculations even when the orbitals were generated using the (very similar, but not exactly the same) pseudopotentials fitted to Gaussians. The reasons for this are:
(1) The tabulated pseudopotentials reproduce the orbitals and eigenvalues they were calculated from slightly more accurately than the Gaussian fits.
(2) The tabulated pseudopotentials have smaller localisation radii (the radius outside of which the components of a pseudopotential differ by less than some tolerance). The cost of evaluating the nonlocal energy in QMC calculations is lower if the localisation
radius is smaller.
The orbitals from HartreeFock calculations using the tabulated pseudopotentials are also given. These are in the correct format for CASINO, and are given for a number of configurations, including the ground state.
Look in the utils/pseudo_converters directory in the CASINO distribution for utilities that convert pseudopotentials formatted for other codes into the correct format for CASINO. Utilities are currently available for ABINIT, CASTEP, CHAMP, PWSCF, and GP.