<div dir="ltr">Also for ngrid 1 the output is:<div>shouldn't we not see any changes less than 10e-7 in energy with changing cutoff?</div><div>I relay appreciate any help here<br><div><br><div> 50.0000000000 -230.3665655634 0</div><div> 100.0000000000 -231.5201721491 0</div><div> 150.0000000000 -231.6333696011 0</div><div> 200.0000000000 -231.6954662742 0</div><div> 250.0000000000 -231.7211472374 0</div><div> 300.0000000000 -231.7441677975 0</div><div> 350.0000000000 -231.7546966069 0</div><div> 400.0000000000 -231.7608341845 0</div><div> 450.0000000000 -231.7712553463 0</div><div> 500.0000000000 -231.7712553463 0</div><div> 550.0000000000 -231.7762288308 0</div><div> 600.0000000000 -231.7833659880 0</div><div><br></div>On Monday, April 23, 2018 at 6:55:33 PM UTC+2, simin pahlavi wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div><br></div><div>Dear Matteu,<br></div><div>Thank you very much for your quick answer,</div><div>I am not sure what did you mean by " Your set up is converged at the lowest setting you report",</div><div> because after first few step with convergence I had a different energy in next cutoffs. </div><div>but the result for lower cutoff is as the following.</div><div>Can I ask where is the suitable cutoff is now? and why you say so?</div><div><br></div><div><div> Grid cutoff vs total energy</div><div># Date: Mon Apr 23 18:44:42 CEST 2018</div><div># PWD: /u/nsaf/JOBS/Ad-Ml.1/Cutoff-<wbr>test</div><div># REL_CUTOFF = 60</div><div># Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on grid 3 | NG on grid 4</div><div> 50.0 -230.3665656028 88679 13322 48 0</div><div> 100.0 -231.5201721496 69162 25725 7162 0</div><div> 150.0 231.6333696411 66211 22468 13322 48</div><div> 200.0 -231.6954662441 48271 28520 25210 48</div><div> 250.0 -231.7211472373 45554 27891 28476 128</div><div> 300.0 -231.7441677949 43072 26090 25725 7162</div><div> 350.0 -231.7546966111 43072 23982 22834 12161</div><div> 400.0 -231.7608341848 40587 25752 23501 12209</div><div> 450.0 -231.7712553838 36594 29617 22468 13370</div><div> 500.0 -231.7712553419 31682 21399 35598 13370</div><div> 550.0 -231.7762288293 24651 26716 28756 21926</div><div> 600.0 -231.7833659579 24651 23620 28520 25258</div></div><div><br></div><div>with best regards</div><div>Simin</div><br>On Monday, April 23, 2018 at 2:57:15 PM UTC+2, 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"><div>To see really clean convergence you need to turn off multigrids (set ngrids 1).<br></div><div><br></div><div>Your set up is converged at the lowest setting you report.</div><div><br></div><div>We need a new tutorial on how to get happy pw settings.</div><br>On Monday, April 23, 2018 at 11:38:03 AM UTC+1, simin pahlavi wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear aal,<div>Hello. </div><div>I am trying to find a suitable cutoff for my DFT calculation on hydrocarbons according to the cp2k manual </div><div><a href="https://www.cp2k.org/howto:converging_cutoff" rel="nofollow" target="_blank" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFHeI16YFpXApR5fGmB8YOPSHi7nQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFHeI16YFpXApR5fGmB8YOPSHi7nQ';return true;">https://www.cp2k.org/howto:<wbr>converging_cutoff</a>. I appriciate any help in advance.<br></div><div>The results never converge and the calculated output for REL_CUTOFF = 60, is as following:</div><div> </div><div><div># Grid cutoff vs total energy</div><div># Date: Fri Apr 20 20:39:04 CEST 2018</div><div># PWD: /u/nsaf/JOBS/Ad-Ml.1/Cutoff-<wbr>test</div><div># REL_CUTOFF = 60</div><div># Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on grid 3 | NG on grid 4</div><div> 800.0 -231.7994673423 73445 28476 128 0</div><div> 850.0 -231.7994673423 73445 28476 128 0</div><div> 900.0 -231.8023753152 69162 25725 7162 0</div><div> 950.0 -231.8023753152 69162 25725 7162 0</div><div>1000.0 -231.8065368280 69162 25725 7162 0</div><div>1050.0 -231.8065368280 69162 25725 7162 0</div><div>1100.0 -231.8065368280 67054 27833 7162 0</div><div>1150.0 -231.8129239556 67054 22834 12161 0</div><div>1200.0 -231.8129239556 67054 22834 12161 0</div><div>1250.0 -231.8129239556 67054 22834 12161 0</div><div>1300.0 -231.8129239556 66339 23501 12209 0</div><div>1350.0 -231.8129239556 66339 23501 12161 48</div><div>1400.0 -231.8148206637 66259 22420 13322 48</div><div>1450.0 -231.8148206637 66211 22468 13322 48</div><div>1500.0 -231.8148206637 66211 22468 13322 48</div><div>1550.0 -231.8158559068 63323 25356 13322 48</div><div>1600.0 -231.8175292623 57217 31462 13322 48</div></div><div>******************************<wbr>******************************<wbr>********************</div><div>Also my input file is :</div><div>******************************<wbr>******************************<wbr>******************</div><div><div>&GLOBAL</div><div> PRINT_LEVEL MEDIUM</div><div># WALLTIME 36:00:00 # changed</div><div> PROJECT cuttoff-test</div><div> RUN_TYPE ENERGY</div><div> SEED 300</div><div>&END GLOBAL</div><div><br></div><div><br></div><div>&FORCE_EVAL</div><div> METHOD Quickstep</div><div> &DFT</div><div> BASIS_SET_FILE_NAME ../data/BASIS_MOLOPT</div><div> POTENTIAL_FILE_NAME ../data/GTH_POTENTIALS</div><div> CHARGE 0</div><div> MULTIPLICITY 1</div><div> &MGRID</div><div> NGRIDS 4</div><div> CUTOFF LT_cutoff</div><div> REL_CUTOFF LT_rel_cutoff</div><div> &END</div><div> &QS</div><div> EPS_DEFAULT 1.0E-10</div><div> WF_INTERPOLATION PS</div><div> EXTRAPOLATION_ORDER 3</div><div> METHOD GPW</div><div> &END</div><div><br></div><div> &SCF</div><div> SCF_GUESS ATOMIC</div><div> EPS_SCF 5.e-7</div><div> MAX_SCF 1</div><div> &OT</div><div> MINIMIZER DIIS</div><div> &END OT</div><div> &OUTER_SCF</div><div> TYPE NONE</div><div> OPTIMIZER NONE</div><div> EPS_SCF 5.0E-7 </div><div> MAX_SCF 1</div><div> &END OUTER_SCF </div><div> &END SCF</div><div><br></div><div> &XC</div><div> &XC_FUNCTIONAL PBE</div><div> &END XC_FUNCTIONAL</div><div> &XC_GRID</div><div> XC_SMOOTH_RHO NN50</div><div> XC_DERIV NN50_SMOOTH</div><div> &END XC_GRID</div><div> &VDW_POTENTIAL </div><div> DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div> &PAIR_POTENTIAL</div><div> TYPE DFTD3</div><div> PARAMETER_FILE_NAME ../data/dftd3.dat</div><div> REFERENCE_FUNCTIONAL PBE</div><div> CALCULATE_C9_TERM TRUE</div><div> REFERENCE_C9_TERM TRUE</div><div> LONG_RANGE_CORRECTION TRUE</div><div> VERBOSE_OUTPUT TRUE</div><div> R_CUTOFF 15.0</div><div> &END PAIR_POTENTIAL </div><div> &END VDW_POTENTIAL</div><div> &END XC</div><div> &END DFT</div><div><br></div><div> &SUBSYS</div><div> &CELL</div><div> ABC 16.0 16.0 16.0</div><div> PERIODIC XYZ</div><div> &END CELL</div><div><br></div><div> &KIND H</div><div> BASIS_SET DZVP-MOLOPT-SR-GTH-q1</div><div> POTENTIAL GTH-PBE-q1</div><div> &END</div><div><br></div><div> &KIND C</div><div> BASIS_SET DZVP-MOLOPT-SR-GTH-q4</div><div> POTENTIAL GTH-PBE-q4</div><div> &END KIND</div><div><br></div><div><br></div><div> &TOPOLOGY</div><div> COORDINATE XYZ</div><div> COORD_FILE_NAME ../data/structure </div><div> CONNECTIVITY OFF</div><div> &END TOPOLOGY</div><div> &END SUBSYS</div><div>&END FORCE_EVAL</div></div><div><br></div><div><br></div><div><br></div></div></blockquote></div></blockquote></div></blockquote></div></div></div>