<div dir="ltr">Hi,<div><br></div><div>it might be best to post a bit of one of your pdos files to make it easier to help.</div><div><br></div><div>You say about a half the electrons seem to be missing?  Is it possible you running an unrestricted calculation and the other half is in the other spin channel - then there would be two pdos files, one for each spin? Without more output/input it is impossible to tell.Do you get the correct number of electrons overall when you do the integration?</div><div><br></div><div>Matt <br><br>On Thursday, April 3, 2014 2:59:18 PM UTC+1, Michael Hartmann wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear all,<div><br></div><div>I have just started using cp2k and am trying to calculate projected density of states for 13 atom platinum clusters.  My input is as follows:</div><div><br></div><div>&PRINT<br></div><div><div>     &PDOS</div><div>           NLUMO -1</div><div>           &LDOS</div><div>                COMPONENTS .TRUE.</div><div>                LIST 1</div><div>           &END LDOS</div><div>     &END PDOS</div><div>&END PRINT</div></div><div><br></div><div>This gives me all of the information I am looking for, but the results do not make any sense to me.  When I integrate the output to the fermi level that is printed in the file that lists each angular momentum component, I get an average d occupation for each platinum atom of about 4.6 electrons, about half of what is expected.  I have done similar calculations on gold clusters with other programs (turbomole, quantum espresso, gpaw) and have not had problems.  I don't understand what could be wrong, am I missing something in the input or miss interpreting the output?</div><div><br></div><div>Thank you</div></div></blockquote></div></div>