<div dir="auto">Thanks a lot for your suggestions. They worked.<div dir="auto">eftrsd</div></div><div class="gmail_extra"><br><div class="gmail_quote">On 29 Oct 2017 22:08, "Matt W" <<a href="mailto:mattwa...@gmail.com">mattwa...@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>one possible problem is that you are running a periodic array of endo fullerenes, rather than an single isolated one (do I understand correctly you just want to run the individual fullerence really?). There could be some attractive electostatic / vdW interactions from that (though I wouldn't expect 0.8 eV!).</div><div><br></div><div>If you set PERIODIC NONE in the CELL section, and add </div><div><br></div><div>&POISSON</div><div>   PSOLVER WAVELET</div><div>   PERIODIC NONE<br>&END</div><div><br></div><div>in the DFT section you will run an isolated molecule. Your molecule needs to be centred in the box if you use this poisson solver. You can add </div><div><br></div><div>&TOPOLOGY</div><div>    &CENTER_COORDINATES</div><div><br></div><div>to do this automatically.</div><div><br></div><div>Otherwise at a quick look your setup looks to have very (meaning over the top) accurate settings.</div><div><br></div><div>Matt</div><div><br><br>On Thursday, October 26, 2017 at 1:52:46 PM UTC+1, eft rsd wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear all,<br>I want to calculate binding energy of an endofullerene. The routine to do so is clear to me but I am getting energy 20 kcal/mol higher than the reported values. I used the highest basis set available in CP2K with the vdW corrected PBE. The script I am using for geometry optimization is given at the end. The optimized geometries of the individual and supra molecules I get using these inputs are very much similar to the ones given by other workers. It is just the high binding energy that I am not getting correct even after correcting for BSSE. I cannot figure out what the reason could be.. Could someone please suggest where the problem is?<br>&GLOBAL<br>  PROJECT      ABC<br>  RUN_TYPE     GEO_OPT<br>  PRINT_LEVEL  MEDIUM<br>  WALLTIME     1000000<br>&END GLOBAL<br><br>@SET RESTART      FALSE<br>@SET BASISFILE    ./BASIS_SET<br>@SET PSEUDOFILE   ./GTH_POTENTIALS-05.07.13<br>!@SET WAVEFILE ----------------<br>@SET CUTOFF       1300<br>@SET GRIDS        6<br>@SET SCF_NCYCLES  800<br>@SET SCF_OCYCLES  200<br>@SET SCF_CONV     1E-8<br>@IF ( ${RESTART} == TRUE )<br>  @SET SCF_GUESS RESTART<br>@ENDIF<br>@IF ( ${RESTART} == FALSE )<br>  @SET SCF_GUESS ATOMIC<br>@ENDIF<br>@SET SCF_MINI     CG<br>@SET FUNCTIONAL   PBE<br>@SET OUT_FORM     XYZ<br>@SET OUT_UNIT     angstrom<br>@SET OUT_STEPS    1<br>&FORCE_EVAL<br>  &PRINT<br>    &FORCES<br>    &END<br>  &END<br>  METHOD Quickstep<br>  &DFT<br>    BASIS_SET_FILE_NAME ${BASISFILE}<br>    POTENTIAL_FILE_NAME ${PSEUDOFILE}<br>@IF ( ${RESTART} == TRUE )<br>    WFN_RESTART_FILE_NAME ${WAVEFILE}<br>@ENDIF<br>    &MGRID<br>      CUTOFF ${CUTOFF}          <br>      NGRIDS ${GRIDS}           <br>      REL_CUTOFF 60<br>    &END MGRID<br>    &QS<br>     METHOD GPW               <br>     EPS_DEFAULT 1.0E-10        <br>     MAP_CONSISTENT TRUE        <br>     EXTRAPOLATION ASPC         <br>     EXTRAPOLATION_ORDER 3     <br>    &END QS<br>    &SCF<br>      SCF_GUESS ${SCF_GUESS}<br>      EPS_SCF 1E-8<br>      MAX_SCF 800<br>      ADDED_MOS 900<br>      CHOLESKY INVERSE<br>      &SMEAR ON<br>        METHOD FERMI_DIRAC<br>        ELECTRONIC_TEMPERATURE [K] 300<br>      &END SMEAR<br>      &DIAGONALIZATION<br>        ALGORITHM STANDARD<br>      &END DIAGONALIZATION<br>      &MIXING<br>        METHOD BROYDEN_MIXING<br>        ALPHA 0.1<br>        BETA 1.5<br>        NBROYDEN 8<br>      &END MIXING<br>      &OUTER_SCF<br>        EPS_SCF 5.0E-9<br>        MAX_SCF 50<br>      &END OUTER_SCF<br>    &END SCF<br>    &XC<br>      &XC_FUNCTIONAL ${FUNCTIONAL}<br>      &END XC_FUNCTIONAL<br>      &XC_GRID<br>        XC_DERIV SPLINE2<br>      &END<br>      &VDW_POTENTIAL<br>        POTENTIAL_TYPE PAIR_POTENTIAL<br>        &PAIR_POTENTIAL<br>          TYPE DFTD3<br>          PARAMETER_FILE_NAME ./dftd3.dat<br>          REFERENCE_FUNCTIONAL PBE<br>          D3_SCALING 1.0 1.217 0.722<br>        &END PAIR_POTENTIAL<br>      &END VDW_POTENTIAL<br>    &END XC<br>   &END DFT<br>&SUBSYS<br>    &KIND C<br>        ELEMENT C<br>        BASIS_SET TZV2PX-MOLOPT-GTH-q4<br>        POTENTIAL GTH-PBE-q4<br>    &END KIND<br>    &KIND O<br>        ELEMENT O<br>        BASIS_SET TZV2PX-MOLOPT-GTH-q6<br>        POTENTIAL GTH-PBE-q6<br>    &END KIND<br>    &KIND H<br>        ELEMENT H<br>        BASIS_SET TZV2PX-MOLOPT-GTH-q1<br>        POTENTIAL GTH-PBE-q1<br>    &END KIND<br>  &CELL<br>A    20.000000000    0.000000000    0.000000000<br>B     0.000000000   20.000000000    0.000000000<br>C     0.000000000    0.000000000   20.000000000<br>  &END CELL<br>  &COORD<br> C    8.727500000   10.302462366   12.600119727<br> C    8.000000000   11.479582093   10.695500000<br> C    7.272500000   10.302462366   12.600119727<br> C    9.177119727   11.029962366   11.423000000<br> C    6.822880273   11.029962366   11.423000000<br> C    8.727500000    5.697537634   12.600119727<br> C    8.000000000    4.520417907   10.695500000<br> C    7.272500000    5.697537634   12.600119727<br> C    6.822880273    4.970037634   11.423000000<br> C    9.177119727    4.970037634   11.423000000<br> C    9.423000000    9.177119727   13.029962366<br> C    9.423000000    6.822880273   13.029962366<br> C    8.695500000    8.000000000   13.479582093<br> C   10.600119727    8.727500000   12.302462366<br> C   10.600119727    7.272500000   12.302462366<br> C    9.423000000    6.822880273    6.970037634<br> C    9.423000000    9.177119727    6.970037634<br> C    8.695500000    8.000000000    6.520417907<br> C   10.600119727    8.727500000    7.697537634<br> C   10.600119727    7.272500000    7.697537634<br> C    8.000000000    4.520417907    9.304500000<br> C    8.727500000    5.697537634    7.399880273<br> C    7.272500000    5.697537634    7.399880273<br> C    6.822880273    4.970037634    8.577000000<br> C    9.177119727    4.970037634    8.577000000<br> C    8.000000000   11.479582093    9.304500000<br> C    8.727500000   10.302462366    7.399880273<br> C    7.272500000   10.302462366    7.399880273<br> C    9.177119727   11.029962366    8.577000000<br> C    6.822880273   11.029962366    8.577000000<br> C    7.304500000    8.000000000    6.520417907<br> C    6.577000000    6.822880273    6.970037634<br> C    6.577000000    9.177119727    6.970037634<br> C    5.399880273    8.727500000    7.697537634<br> C    5.399880273    7.272500000    7.697537634<br> C    6.577000000    9.177119727   13.029962366<br> C    6.577000000    6.822880273   13.029962366<br> C    7.304500000    8.000000000   13.479582093<br> C    5.399880273    8.727500000   12.302462366<br> C    5.399880273    7.272500000   12.302462366<br> C   10.302462366   10.600119727   10.727500000<br> C   11.029962366    9.423000000   11.177119727<br> C   11.029962366    9.423000000    8.822880273<br> C   10.302462366   10.600119727    9.272500000<br> C   11.479582093    8.695500000   10.000000000<br> C    5.697537634   10.600119727   10.727500000<br> C    5.697537634   10.600119727    9.272500000<br> C    4.970037634    9.423000000    8.822880273<br> C    4.970037634    9.423000000   11.177119727<br> C    4.520417907    8.695500000   10.000000000<br> C    5.697537634    5.399880273   10.727500000<br> C    5.697537634    5.399880273    9.272500000<br> C    4.970037634    6.577000000    8.822880273<br> C    4.970037634    6.577000000   11.177119727<br> C    4.520417907    7.304500000   10.000000000<br> C   10.302462366    5.399880273   10.727500000<br> C   11.029962366    6.577000000   11.177119272<br> C   11.029962366    6.577000000    8.822880273<br> C   10.302462366    5.399880273    9.272500000<br> C   11.479582093    7.304500000   10.000000000<br> O    8.000000000    8.000000000   10.000000000<br> H    8.757480612    8.000000000   10.586504154<br> H    7.242519388    8.000000000   10.586504154<br>  &END COORD<br>&END SUBSYS<br>&END FORCE_EVAL<br>@SET GEO_MINI     BFGS<br>@SET GEO_MAXS     1000<br>&MOTION<br>  &GEO_OPT<br>    MINIMIZER ${GEO_MINI}<br>    MAX_ITER ${GEO_MAXS}<br>    MAX_DR 0.001<br>    MAX_FORCE 0.0000097<br>    RMS_DR 0.001<br>    RMS_FORCE 0.001<br>    &BFGS<br>      USE_MODEL_HESSIAN TRUE<br>      USE_RAT_FUN_OPT TRUE<br>      TRUST_RADIUS 0.1<br>    &END<br>  &END GEO_OPT<br><br>  &PRINT<br>    &RESTART<br>      LOG_PRINT_KEY T<br>      &EACH<br>        GEO_OPT ${OUT_STEPS}<br>      &END EACH<br>      ADD_LAST NUMERIC<br>    &END RESTART<br>    &TRAJECTORY<br>      LOG_PRINT_KEY T<br>      FORMAT ${OUT_FORM}<br>      UNIT ${OUT_UNIT}<br>      &EACH<br>        GEO_OPT ${OUT_STEPS}  <br>      &END EACH<br>      ADD_LAST NUMERIC<br>    &END TRAJECTORY<br>  &END PRINT<br>&END MOTION<br><br>@IF ( ${RESTART} == TRUE )<br>  &EXT_RESTART ON<br>    RESTART_DEFAULT F<br>    RESTART_FILE_NAME ${RESTARTFILE}<br>    RESTART_POS T<br>    RESTART_COUNTERS T<br>  &END EXT_RESTART<br>@ENDIF<br></div></blockquote></div></div>

<p></p>

-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/hm_x1O-CliA/unsubscribe" target="_blank">https://groups.google.com/d/<wbr>topic/cp2k/hm_x1O-CliA/<wbr>unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank">https://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/<wbr>optout</a>.<br>
</blockquote></div></div>