<div dir="ltr">Hi All,<br><br>I am trying to perform ab initio MD of pure sodium metal using CP2K. Attached is my input and output files for the same. Since for pure sodium I am using a semi core PP with nine valance electrons the required energy cut off is relatively higher (~ 1200 Ry). I am using diagonalization and fermi smearing at 300 K which is recommended for metallic systems. But still after a few fs of the run I find that the sodium atoms begin to lose its shape within the crystal. Also the fermi energy shown in the output while performing a geometry optimization is only 1/3 of that reported for pure sodium.<br><br>Kindly let me know what parameters I am setting wrong in doing these calculations? <br><br>Regards,<br>Aniruddha M Dive<br>PhD Candidate<br>Washington State University<br></div>