<div dir="ltr">Dear all, 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? &GLOBAL Â PROJECTÂ Â Â Â Â ABC Â RUN_TYPEÂ Â Â Â GEO_OPT Â PRINT_LEVELÂ MEDIUM Â WALLTIMEÂ Â Â Â 1000000 &END GLOBAL @SET RESTARTÂ Â Â Â Â FALSE @SET BASISFILEÂ Â Â ./BASIS_SET @SET PSEUDOFILEÂ Â ./GTH_POTENTIALS-05.07.13 !@SET WAVEFILE ---------------- @SET CUTOFFÂ Â Â Â Â Â 1300 @SET GRIDSÂ Â Â Â Â Â Â 6 @SET SCF_NCYCLESÂ 800 @SET SCF_OCYCLESÂ 200 @SET SCF_CONVÂ Â Â Â 1E-8 @IF ( ${RESTART} == TRUE ) Â @SET SCF_GUESS RESTART @ENDIF @IF ( ${RESTART} == FALSE ) Â @SET SCF_GUESS ATOMIC @ENDIF @SET SCF_MINIÂ Â Â Â CG @SET FUNCTIONALÂ Â PBE @SET OUT_FORMÂ Â Â Â XYZ @SET OUT_UNITÂ Â Â Â angstrom @SET OUT_STEPSÂ Â Â 1 &FORCE_EVAL Â &PRINT Â Â Â &FORCES Â Â Â &END Â &END Â METHOD Quickstep Â &DFT Â Â Â BASIS_SET_FILE_NAME ${BASISFILE} Â Â Â POTENTIAL_FILE_NAME ${PSEUDOFILE} @IF ( ${RESTART} == TRUE ) Â Â Â WFN_RESTART_FILE_NAME ${WAVEFILE} @ENDIF Â Â Â &MGRID Â Â Â Â Â CUTOFF ${CUTOFF}Â Â Â Â Â Â Â Â Â Â Â Â Â Â NGRIDS ${GRIDS}Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â REL_CUTOFF 60 Â Â Â &END MGRID Â Â Â &QS Â Â Â Â METHOD GPWÂ Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â EPS_DEFAULT 1.0E-10Â Â Â Â Â Â Â Â Â Â Â MAP_CONSISTENT TRUEÂ Â Â Â Â Â Â Â Â Â Â EXTRAPOLATION ASPCÂ Â Â Â Â Â Â Â Â Â Â Â EXTRAPOLATION_ORDER 3Â Â Â Â Â Â Â &END QS Â Â Â &SCF Â Â Â Â Â SCF_GUESS ${SCF_GUESS} Â Â Â Â Â EPS_SCF 1E-8 Â Â Â Â Â MAX_SCF 800 Â Â Â Â Â ADDED_MOS 900 Â Â Â Â Â CHOLESKY INVERSE Â Â Â Â Â &SMEAR ON Â Â Â Â Â Â Â METHOD FERMI_DIRAC Â Â Â Â Â Â Â ELECTRONIC_TEMPERATURE [K] 300 Â Â Â Â Â &END SMEAR Â Â Â Â Â &DIAGONALIZATION Â Â Â Â Â Â Â ALGORITHM STANDARD Â Â Â Â Â &END DIAGONALIZATION Â Â Â Â Â &MIXING Â Â Â Â Â Â Â METHOD BROYDEN_MIXING Â Â Â Â Â Â Â ALPHA 0.1 Â Â Â Â Â Â Â BETA 1.5 Â Â Â Â Â Â Â NBROYDEN 8 Â Â Â Â Â &END MIXING Â Â Â Â Â &OUTER_SCF Â Â Â Â Â Â Â EPS_SCF 5.0E-9 Â Â Â Â Â Â Â MAX_SCF 50 Â Â Â Â Â &END OUTER_SCF Â Â Â &END SCF Â Â Â &XC Â Â Â Â Â &XC_FUNCTIONAL ${FUNCTIONAL} Â Â Â Â Â &END XC_FUNCTIONAL Â Â Â Â Â &XC_GRID Â Â Â Â Â Â Â XC_DERIV SPLINE2 Â Â Â Â Â &END Â Â Â Â Â &VDW_POTENTIAL Â Â Â Â Â Â Â POTENTIAL_TYPE PAIR_POTENTIAL Â Â Â Â Â Â Â &PAIR_POTENTIAL Â Â Â Â Â Â Â Â Â TYPE DFTD3 Â Â Â Â Â Â Â Â Â PARAMETER_FILE_NAME ./dftd3.dat Â Â Â Â Â Â Â Â Â REFERENCE_FUNCTIONAL PBE Â Â Â Â Â Â Â Â Â D3_SCALING 1.0 1.217 0.722 Â Â Â Â Â Â Â &END PAIR_POTENTIAL Â Â Â Â Â &END VDW_POTENTIAL Â Â Â &END XC Â Â &END DFT &SUBSYS Â Â Â &KIND C Â Â Â Â Â Â Â ELEMENT C Â Â Â Â Â Â Â BASIS_SET TZV2PX-MOLOPT-GTH-q4 Â Â Â Â Â Â Â POTENTIAL GTH-PBE-q4 Â Â Â &END KIND Â Â Â &KIND O Â Â Â Â Â Â Â ELEMENT O Â Â Â Â Â Â Â BASIS_SET TZV2PX-MOLOPT-GTH-q6 Â Â Â Â Â Â Â POTENTIAL GTH-PBE-q6 Â Â Â &END KIND Â Â Â &KIND H Â Â Â Â Â Â Â ELEMENT H Â Â Â Â Â Â Â BASIS_SET TZV2PX-MOLOPT-GTH-q1 Â Â Â Â Â Â Â POTENTIAL GTH-PBE-q1 Â Â Â &END KIND Â &CELL AÂ Â Â 20.000000000Â Â Â 0.000000000Â Â Â 0.000000000 BÂ Â Â Â 0.000000000Â Â 20.000000000Â Â Â 0.000000000 CÂ Â Â Â 0.000000000Â Â Â 0.000000000Â Â 20.000000000 Â &END CELL Â &COORD Â CÂ Â Â 8.727500000Â Â 10.302462366Â Â 12.600119727 Â CÂ Â Â 8.000000000Â Â 11.479582093Â Â 10.695500000 Â CÂ Â Â 7.272500000Â Â 10.302462366Â Â 12.600119727 Â CÂ Â Â 9.177119727Â Â 11.029962366Â Â 11.423000000 Â CÂ Â Â 6.822880273Â Â 11.029962366Â Â 11.423000000 Â CÂ Â Â 8.727500000Â Â Â 5.697537634Â Â 12.600119727 Â CÂ Â Â 8.000000000Â Â Â 4.520417907Â Â 10.695500000 Â CÂ Â Â 7.272500000Â Â Â 5.697537634Â Â 12.600119727 Â CÂ Â Â 6.822880273Â Â Â 4.970037634Â Â 11.423000000 Â CÂ Â Â 9.177119727Â Â Â 4.970037634Â Â 11.423000000 Â CÂ Â Â 9.423000000Â Â Â 9.177119727Â Â 13.029962366 Â CÂ Â Â 9.423000000Â Â Â 6.822880273Â Â 13.029962366 Â CÂ Â Â 8.695500000Â Â Â 8.000000000Â Â 13.479582093 Â CÂ Â 10.600119727Â Â Â 8.727500000Â Â 12.302462366 Â CÂ Â 10.600119727Â Â Â 7.272500000Â Â 12.302462366 Â CÂ Â Â 9.423000000Â Â Â 6.822880273Â Â Â 6.970037634 Â CÂ Â Â 9.423000000Â Â Â 9.177119727Â Â Â 6.970037634 Â CÂ Â Â 8.695500000Â Â Â 8.000000000Â Â Â 6.520417907 Â CÂ Â 10.600119727Â Â Â 8.727500000Â Â Â 7.697537634 Â CÂ Â 10.600119727Â Â Â 7.272500000Â Â Â 7.697537634 Â CÂ Â Â 8.000000000Â Â Â 4.520417907Â Â Â 9.304500000 Â CÂ Â Â 8.727500000Â Â Â 5.697537634Â Â Â 7.399880273 Â CÂ Â Â 7.272500000Â Â Â 5.697537634Â Â Â 7.399880273 Â CÂ Â Â 6.822880273Â Â Â 4.970037634Â Â Â 8.577000000 Â CÂ Â Â 9.177119727Â Â Â 4.970037634Â Â Â 8.577000000 Â CÂ Â Â 8.000000000Â Â 11.479582093Â Â Â 9.304500000 Â CÂ Â Â 8.727500000Â Â 10.302462366Â Â Â 7.399880273 Â CÂ Â Â 7.272500000Â Â 10.302462366Â Â Â 7.399880273 Â CÂ Â Â 9.177119727Â Â 11.029962366Â Â Â 8.577000000 Â CÂ Â Â 6.822880273Â Â 11.029962366Â Â Â 8.577000000 Â CÂ Â Â 7.304500000Â Â Â 8.000000000Â Â Â 6.520417907 Â CÂ Â Â 6.577000000Â Â Â 6.822880273Â Â Â 6.970037634 Â CÂ Â Â 6.577000000Â Â Â 9.177119727Â Â Â 6.970037634 Â CÂ Â Â 5.399880273Â Â Â 8.727500000Â Â Â 7.697537634 Â CÂ Â Â 5.399880273Â Â Â 7.272500000Â Â Â 7.697537634 Â CÂ Â Â 6.577000000Â Â Â 9.177119727Â Â 13.029962366 Â CÂ Â Â 6.577000000Â Â Â 6.822880273Â Â 13.029962366 Â CÂ Â Â 7.304500000Â Â Â 8.000000000Â Â 13.479582093 Â CÂ Â Â 5.399880273Â Â Â 8.727500000Â Â 12.302462366 Â CÂ Â Â 5.399880273Â Â Â 7.272500000Â Â 12.302462366 Â CÂ Â 10.302462366Â Â 10.600119727Â Â 10.727500000 Â CÂ Â 11.029962366Â Â Â 9.423000000Â Â 11.177119727 Â CÂ Â 11.029962366Â Â Â 9.423000000Â Â Â 8.822880273 Â CÂ Â 10.302462366Â Â 10.600119727Â Â Â 9.272500000 Â CÂ Â 11.479582093Â Â Â 8.695500000Â Â 10.000000000 Â CÂ Â Â 5.697537634Â Â 10.600119727Â Â 10.727500000 Â CÂ Â Â 5.697537634Â Â 10.600119727Â Â Â 9.272500000 Â CÂ Â Â 4.970037634Â Â Â 9.423000000Â Â Â 8.822880273 Â CÂ Â Â 4.970037634Â Â Â 9.423000000Â Â 11.177119727 Â CÂ Â Â 4.520417907Â Â Â 8.695500000Â Â 10.000000000 Â CÂ Â Â 5.697537634Â Â Â 5.399880273Â Â 10.727500000 Â CÂ Â Â 5.697537634Â Â Â 5.399880273Â Â Â 9.272500000 Â CÂ Â Â 4.970037634Â Â Â 6.577000000Â Â Â 8.822880273 Â CÂ Â Â 4.970037634Â Â Â 6.577000000Â Â 11.177119727 Â CÂ Â Â 4.520417907Â Â Â 7.304500000Â Â 10.000000000 Â CÂ Â 10.302462366Â Â Â 5.399880273Â Â 10.727500000 Â CÂ Â 11.029962366Â Â Â 6.577000000Â Â 11.177119272 Â CÂ Â 11.029962366Â Â Â 6.577000000Â Â Â 8.822880273 Â CÂ Â 10.302462366Â Â Â 5.399880273Â Â Â 9.272500000 Â CÂ Â 11.479582093Â Â Â 7.304500000Â Â 10.000000000 Â OÂ Â Â 8.000000000Â Â Â 8.000000000Â Â 10.000000000 Â HÂ Â Â 8.757480612Â Â Â 8.000000000Â Â 10.586504154 Â HÂ Â Â 7.242519388Â Â Â 8.000000000Â Â 10.586504154 Â &END COORD &END SUBSYS &END FORCE_EVAL @SET GEO_MINIÂ Â Â Â BFGS @SET GEO_MAXSÂ Â Â Â 1000 &MOTION Â &GEO_OPT Â Â Â MINIMIZER ${GEO_MINI} Â Â Â MAX_ITER ${GEO_MAXS} Â Â Â MAX_DR 0.001 Â Â Â MAX_FORCE 0.0000097 Â Â Â RMS_DR 0.001 Â Â Â RMS_FORCE 0.001 Â Â Â &BFGS Â Â Â Â Â USE_MODEL_HESSIAN TRUE Â Â Â Â Â USE_RAT_FUN_OPT TRUE Â Â Â Â Â TRUST_RADIUS 0.1 Â Â Â &END Â &END GEO_OPT Â &PRINT Â Â Â &RESTART Â Â Â Â Â LOG_PRINT_KEY T Â Â Â Â Â &EACH Â Â Â Â Â Â Â GEO_OPT ${OUT_STEPS} Â Â Â Â Â &END EACH Â Â Â Â Â ADD_LAST NUMERIC Â Â Â &END RESTART Â Â Â &TRAJECTORY Â Â Â Â Â LOG_PRINT_KEY T Â Â Â Â Â FORMAT ${OUT_FORM} Â Â Â Â Â UNIT ${OUT_UNIT} Â Â Â Â Â &EACH Â Â Â Â Â Â Â GEO_OPT ${OUT_STEPS}Â Â Â Â Â Â &END EACH Â Â Â Â Â ADD_LAST NUMERIC Â Â Â &END TRAJECTORY Â &END PRINT &END MOTION @IF ( ${RESTART} == TRUE ) Â &EXT_RESTART ON Â Â Â RESTART_DEFAULT F Â Â Â RESTART_FILE_NAME ${RESTARTFILE} Â Â Â RESTART_POS T Â Â Â RESTART_COUNTERS T Â &END EXT_RESTART @ENDIF </div>