<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>