<p>Dear CP2K users</p><p> I would like to smear my DOS using gaussians profiles. In the webpage </p><p> <a href="https://www.cp2k.org/exercises:2017_uzh_cmest:pdos" style="text-decoration: underline;">https://www.cp2k.org/exercises:2017_uzh_cmest:pdos</a></p><p> there are two scripts suggested for that: get-smearing-pdos.py from WPI and cp2k_pdos.py by Tiziano/Marcella. </p><p> I tried both scripts on my pdos output (attached at link below). </p><p>> python cp2k_pdos.py --output tiziano FE2-ALPHA_k1-1.pdos FE2-BETA_k1-1.pdos</p><p>>get-smearing-pdos.pyFE2-ALPHA_k1-1.pdos FE2-BETA_k1-1.pdos</p><p>I see that the normalization of the curves is different between the two codes but what worries me is the position of the peaks. In Tiziano script, I see a double peak at  ~-5eV -3 eV that I do not see using the other one.<br></p><p>I tried to change the sigma parameter used in cp2k_pdos.py  to match the profiles with no success</p><p>Moreover, both codes shift the level by the fermi energy (i.e., Fermi level at 0 eV). The Fermi level is reported in the first line of the pdos file and the Fermi level is different for ALPHA and BETA channels. Which is the one used for the shift?</p><p>I attach at the following link the PDOS file:</p><p>https://drive.google.com/drive/folders/17-Hget0l8ezblMmjO5FqBDQB8UJgE36n?usp=sharing</p><p>Thank you very much for any possible help<br></p><p>Best</p><p>Ivan</p>