<div dir="ltr">Hey Matt,<div><br></div><div>Thank you for the answer. I am working on running the calculations with bigger cell size at the moment. In the meantime, I am wondering if there is a way to turn off the periodicity in CP2K. I tried using the keyword PERIODIC NONE but got the exact same results.</div><div><br></div><div>Any help will be much appreciated it.</div><div><br></div><div>Thank you,</div><div>Bobby<br><br>On Monday, January 23, 2017 at 4:59:51 PM UTC-5, 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>it looks like you cell is very small. Unlike Gaussian CP2K tends to run with atoms periodically repeated inside a box. Your's is about 5x5x5 Angstrom^3. Probably you have some steric clashes between images of the molecule, and likely these are different between the conformers.</div><div><br></div><div>Matt<br><br>On Monday, January 23, 2017 at 8:05:24 PM UTC, Babgen Manookian wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello all,<div><br></div><div>I am new to CP2K and am working on making the transfer from Gaussian09. I am trying to match energy calculations between the two programs but am having trouble so I hope someone can help me.</div><div><br></div><div>System: Just for testing purposes, I want to calculate the energy difference between chair and boat conformations for a cyclohexane molecule.</div><div><br></div><div>Gaussian09: for my Gaussian calculations I am using DFT with BLYP functional and a 6-31G* basis set. Using this setup I obtain an energy difference of ~34 kJ/mol which is in fair agreement with a textbook value of ~30 kJ/mol.</div><div><br></div><div>CP2K: For CP2K I am using DZVP-GTH-BLYP basis sets and GTH-BLYP pseudopotentials for both carbon and hydrogen atoms. Pretty much all other input parameters are set to default values or in accordance to the ones given in the Si test files found on the CP2K website. The input file is attached for further detail if needed. With this set-up I am obtaining an energy difference of ~288 kJ/mol. Clearly this value is way off the expected value and my Gaussian-obtained value and I am totally lost as to why.</div><div><br></div><div>I have done these calculations with different values for cutoffs and different basis sets and in all cases I get a similar result where the difference between the two methods is roughly ten-fold. I am wondering if there may be something trivial that I am doing that is causing this.</div><div><br></div><div>Any sort of insight will be truly appreciated.</div><div><br></div><div>Thank you very much in advance.</div><div><br></div><div>Babgen Manookian</div></div></blockquote></div></div></blockquote></div></div>