<div dir="ltr"><div>Hey Xin Wu,</div><div><br></div><div>My name is Julien and I'm a colleague of Salome's (hence I am answering on her behalf).</div><div><br></div><div>Thank you for your suggestions -- your first suggestion is correct and is the source of the discrepancy in the results for our runs and those from the Gaussian tutorial. Both of us come from a classical MD background and hadn't considered *not* running with periodic boundary conditions...</div><div><br></div><div>For good measure, I also tried your other suggestions -- they had little to no effect on the energy differences.</div><div><br></div><div>Thanks again for solving this,</div><div dir="ltr">Julien</div></div>