<div dir="ltr">Hi Matt,<br><br>Yeah my wavefunction is converging at each step. I shall try with NVT and check.<br><br>Thanks,<br>Aniruddha M Dive<br><br>On Thursday, March 23, 2017 at 2:12:38 AM UTC-7, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<div><br></div><div>I'd suggest checking you can run NVT before trying NPT. Your cutoff might not be large enough for the stress tensor (yes I know the values are getting ridiculous).</div><div><br></div><div>Is your wavefunction converging at each step?<br><br>Matt<br><br>On Wednesday, March 22, 2017 at 5:54:36 PM UTC, Aniruddha Dive wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<br><br>I am performing ab - initio MD for pure sodium at 300K which is room temperature. After 1 ps run, I find that the crystal structure of the pure sodium does not remain purely crystalline. I am not sure as to why this is happening. There shall be certain vibrations of the sodium ions but the rdf plot becomes more of an amorphous type after 1 ps of MD run. The sodium uses semi core basis sets which has 9 valence electrons and hence I need to use a higher value of cutoff (~ 950 Ry) which is equivalent to 12925 eV. Is this high value of cutoff correct? Or am I just using too high a cutoff and as a result of which the sodium atoms are behaving as liquid sodium?<br><br>I have attached my input file and trajectory files hereby. <br><br>Kindly let me know what am I doing wrong.<br><br>Thanks and Regards,<br>Aniruddha M Dive<br>PhD Student<br>Washington State University<br></div></blockquote></div></div></blockquote></div>