<div dir="ltr">Hi Matt,<br><br>thanks for your answer. The issue was caused by wrongly formatted psf file - last two indices for torsion angles were pointing to the same atoms. After correction calcs start running OK but after while they crash. I need to make more tries but it seems that there are no problems if I use default units. All the parameters are taken from Pavese and Parker's potential  which I used with DL_POLY, but maybe it will be better to move them into default units   <br><br>cheers<br><br>Wojtek<br><br>On Tuesday, November 25, 2014 5:35:57 PM UTC+1, 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 Wojtek,<div><font size="2"><span style="white-space:nowrap"><br></span></font></div><div>not sure exactly, but:</div><div><span style="font-size:small;white-space:nowrap">(i) your bend potential is ludicrously strong - think you've probably misconverted units</span></div><div><span style="font-size:small;white-space:nowrap">(ii) not sure the ordering of atoms in your torsion is what you want - or if you want a pure torsion for a planar trigonal molecule... </span></div><div><font size="2"><span style="white-space:nowrap"><br></span></font></div><div><font size="2"><span style="white-space:nowrap">HTH you track down the problem</span></font></div><div><font size="2"><span style="white-space:nowrap"><br></span></font></div><div><font size="2"><span style="white-space:nowrap">Matt</span></font><div><br>On Wednesday, November 19, 2014 10:35:37 PM UTC, Wojtek 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>I'm trying to run calcite force field calculations. I've started from the simplest configuration with just one CaCO3 unit. The first energy evaluation looks OK, but in the next step I'm getting NaN. Also coordinates are NaNs.<br>Interestingly when I exclude angles and dihedrals, the simulation runs without problems. <br>What I'm doing wrong? <br>I'm attaching the input and output files for both cases. <br><br>Thanks <br></div></blockquote></div></div></div></blockquote></div>