We report smooth relativistic Hartree-Fock pseudopotentials (also known as averaged relativistic effective potentials or AREPs) and spin-orbit operators for the atoms H to Ba and Lu to Hg. We remove the unphysical extremely non-local behaviour resulting from the exchange interaction in a controlled manner, and represent the resulting pseudopotentials in an analytic form suitable for use within standard quantum chemistry codes. These pseudopotentials are suitable for use within Hartree-Fock and correlated wave function methods, including diffusion quantum Monte Carlo calculations.PACS numbers: 71.15. Dx, 02.70.Ss Pseudopotentials or effective core potentials (ECPs) are commonly used within electronic structure calculations to replace the chemically inert core electrons. The influence of the core on the valence electrons is then described by an angular-momentum-dependent effective potential, leading to greatly improved computational efficiency in ab initio calculations for heavy atoms. The use of pseudopotentials is well established within HartreeFock (HF) and Density Functional Theory (DFT), and in correlated wave function calculations.Our main interest is in diffusion quantum Monte Carlo (DMC) calculations. 1,2 This technique provides an accurate solution of the interacting electron problem for which the computational effort scales with the number of electrons, N , as approximately N 3 , which is better than other correlated wave function approaches. Unfortunately, scaling with atomic number, Z, is approximately 3,4 Z 5−6.5 . The use of a pseudopotential reduces the effective value of Z, making DMC calculations feasible for heavy atoms.There is evidence that HF pseudopotentials give better results within DMC than DFT pseudopotentials. 5 It appears that the complete neglect of core-valence correlation within HF theory leads to better pseudopotentials than the description of core-valence correlation provided by DFT. Moreover, core-valence correlation can be included within correlated wave function calculations performed with HF pseudopotentials by using core polarization potentials. 6,7,8 Core polarization potentials mimic the effects of dynamical polarization of the core by the valence electrons, as well as static polarization effects due to the other ions. We would therefore like to use HF pseudopotentials in our DMC calculations, preferably constructed from Dirac-Fock (DF) theory in order to include the relativistic effects which are significant for heavy atoms.Standard quantum chemistry packages are convenient for generating the "guiding wave functions" required in DMC calculations. We would therefore like our pseudopotentials to be available in the standard parameterized form of a sum of Gaussian functions multiplied by powers of the electron-nucleus separation.Extensive sets of parameterized pseudopotentials are available in the literature, but they have generally been constructed with different goals to ours. Relativistic pseudopotentials 9 generated within DFT and the local density approximation have be...