Hello CP2K Community,<div><br /><div>I have a system consisting of one pair of [EMIM]+ and [TFSI]- ions for which I want to perform an energy convergence test. I use the all-electron GAPW method along with 6-311G** basis set and BLYP exchange-correlation functional with DFT-D3 correction and no smoothing. I have set EPS_DEFAULT and EPS_SCF to 1E-12 and 1E-8, respectively. Here are the total energies obtained:</div><div><br /></div><div><span cellspacing="0" border="0"><span height="17" align="right" valign="bottom">cutoff(Ry)     Total Energy(a.u.)<br /><font face="Liberation Serif">200                 </font></span><span align="right" valign="bottom"><font face="Liberation Serif">-2172.105218573</font></span>
        
        
                <span height="17" align="right" valign="bottom"><br /><font face="Liberation Serif">400                 </font></span><span align="right" valign="bottom"><font face="Liberation Serif">-2172.1056897835</font></span>
        
        
                <span height="17" align="right" valign="bottom"><br /><font face="Liberation Serif">600                 </font></span><span align="right" valign="bottom"><font face="Liberation Serif">-2172.1056642192</font></span>
        
        
                <span height="17" align="right" valign="bottom"><br /><font face="Liberation Serif">800                 </font></span><span align="right" valign="bottom"><font face="Liberation Serif">-2172.105677308</font></span>
        
        
                <span height="17" align="right" valign="bottom"><br /><font face="Liberation Serif">1000               </font></span><span align="right" valign="bottom"><font face="Liberation Serif">-2172.1056768808</font></span>
        
        
                <span height="17" align="right" valign="bottom"><br /><font face="Liberation Serif">1200               </font></span><span align="right" valign="bottom"><font face="Liberation Serif">-2172.1056755387</font></span>
        
        
                <span height="17" align="right" valign="bottom"><br /><font face="Liberation Serif">1400               </font></span><span align="right" valign="bottom"><font face="Liberation Serif">-2172.1056755104</font></span>
        
        
                <span height="17" align="right" valign="bottom"><br /><font face="Liberation Serif">1600               </font></span><span align="right" valign="bottom"><font face="Liberation Serif">-2172.1056755094</font></span>
        
</span><br /></div><div><span cellspacing="0" border="0"><span align="right" valign="bottom"><font face="Liberation Serif"><br /></font></span></span></div><div><span cellspacing="0" border="0"><span align="right" valign="bottom"><span style="text-align: start;">Although similar studies (studying ionic liquid and using similar setting) use a cutoff energy around 280-300 Ry, my system seems to converge at much higher cutoff values. My question is if there is any way to improve this result. You can find my input and data files attached.</span><font face="Liberation Serif"><br /></font></span></span></div><div><span cellspacing="0" border="0"><span align="right" valign="bottom"><span style="text-align: start;"><br /></span></span></span></div><div><span cellspacing="0" border="0"><span align="right" valign="bottom"><span style="text-align: start;">Thanks,</span></span></span></div><div><span cellspacing="0" border="0"><span align="right" valign="bottom"><span style="text-align: start;">Hana</span></span></span></div><div><br /></div><div><br /></div></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/2e0435d8-63e3-42a4-8494-99d8ef0a59b1n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/2e0435d8-63e3-42a4-8494-99d8ef0a59b1n%40googlegroups.com</a>.<br />