Eur. Phys. J. B (2012) 85: 239
https://doi.org/10.1140/epjb/e2012-30334-5

Regular Article

Calculation of core level shifts within DFT using pseudopotentials and localized basis sets

¹ Centre d’Investigació en Nanociència i Nanotecnologia – CIN2 (CSIC-ICN), Campus de la U.A.B., 08193 Bellaterra, Spain
² Centre d’Élaborationde Matériaux et d’ÉtudesStructurales (CNRS), 29 Rue Jeanne Marvig, 31055 Toulouse Cedex, France
³ Institut des Sciences Moléculaires d’Orsay (CNRS), Université Paris-Sud, Bâtiment 351, 91405 Orsay Cedex, France
⁴ Institut de Ciencia de Materials de Barcelona (ICMAB-CSIC), Campus de la U.A.B., 08193 Bellaterra, Spain

^a e-mail: pablo.ordejon@cin2.es

Received: 19 April 2012
Received in final form: 29 May 2012
Published online: 9 July 2012

Abstract

The calculation of core level shifts can be done in the context of density functional theory (DFT) using different approaches and physical approximations to the photoemission process. The initial state and the ΔSCF approximations are the most commonly used ones. Here, we describe the details of their implementation in the context of DFT using pseudopotentials and localized atomic orbitals as a basis set, and in particular as applied to the Siesta code. We give a full account of the technicalities involved in these calculations, including the details of the ionic pseudopotential generation, basis sets, charge states and reference potential. We test the method by computing the core level shifts of the Si 2p level for a series of molecules and the p(2 × 2) asymmetric-dimer reconstruction of the Si(001) surface.