The 4f electrons of lanthanides, because of their strong localization in the region around the nucleus, are traditionally included in a pseudopotential core. This approximation is scrutinized by optimizing the structures and calculating the interaction energies of Gd3+(H2O) and Gd3+(NH3) microsolvation complexes within plane wave Perdew-Burke-Ernzerhof calculations using ultrasoft pseudopotentials where the 4f electrons are included either in the core or in the valence space. Upon comparison to quantum chemical MP2 and CCSD(T) reference calculations it is found that the explicit treatment of the 4f electrons in the valence shell yields quite accurate results including the required small spin polarization due to ligand charge transfer with only modest computational overhead.
