Hello,<div>I am trying to reproduce the results of the following exercise:</div><div><br></div><div>https://www.cp2k.org/exercises:2019_conexs_newcastle:ex4<br></div><div><br></div><div>using the attached input .inp and compare it to the data given here:</div><div>https://www.cp2k.org/_media/exercises:2019_conexs_newcastle:h2o_pes.pdf<br></div><div><br></div><div>I used the python code that the exercise suggests to draw the plot</div><div>but "s states" have a much lower value than in the data above, the plot is almost flat. (64_water_pdos_fwhm_0.5pdos.png), same after MD (BOMD_pdos_FWHM_0.4.png).</div><div>Why is this happening?</div><div><br></div><div>To get a better look at the shape of s states I multiplied s*20 in python code(BOMD_pdos__20s_FWHM_0.4.png)</div><div>In the paper is mentioned "the 2s contributions were obtained from a spectrum computed using an r^2 transition operator" how could i get the value of this operator?</div><div><br></div><div>Thank you in advance</div><div><br></div><div><br></div><div><br></div>