<div dir="ltr">Dear CP2K experts/users,<div>I'm working on catalysts and have recently started to work with CP2K. For the first step of my journey through CP2K I'm trying to replicate the results for CO2 adsorption on Anatase (001) surface that's been done previously by using VASP in our group.</div><div>Unfortunately, the results that I get for adsorption energy from CP2K are far from the results published in literature (done by other packages such as VASP). </div><div>I'm struggling with finding the source of problem but after trying different modifications (Slab size, Cutoff, Basis sets, Dispersion corrections, etc) still couldn't achieve the correct amount of adsorption energy. </div><div>- The best result that I've achieved is about -38.08 kcal/mol (with DZVP basis sets for C & O) (I tested DZVP, DZVP-SR, TZVP, TZV2P basis sets, but the least variation from the results of literature was for DZVP) </div><div>- The result for adsorption energy achieved by VASP for the same structure is -34.41 kcal/mol (without zero point energy correction)</div><div>- I'm using CP2K version 2.6. the newer versions are not available for me at the moment.</div><div>I would really appreciate if someone could take a look to my input, xyz and output (excel table) files and let me know about any possible problem with them..</div><div><br></div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT c-dzvp</div><div>  RUN_TYPE GEO_OPT</div><div>  PRINT_LEVEL LOW</div><div>  EXTENDED_FFT_LENGTHS</div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL</div><div>  METHOD Quickstep            </div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div>    LSD</div><div>    MULTIPLICITY 1</div><div>    CHARGE 0</div><div>    &MGRID</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>CUTOFF 1200</div><div><span class="Apple-tab-span" style="white-space:pre">  </span>REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &QS</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>METHOD GPW</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>EPS_DEFAULT 1.0E-10</div><div>    &END QS</div><div>    &POISSON</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>PERIODIC XYZ</div><div>    &END POISSON</div><div>    &SCF                </div><div>      SCF_GUESS</div><div>      EPS_SCF 1.0E-6</div><div>      MAX_SCF 140</div><div>      &OT</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>PRECONDITIONER FULL_ALL</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>MINIMIZER CG</div><div>      &END OT</div><div>      &OUTER_SCF</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>MAX_SCF 20</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>EPS_SCF 1.0E-6</div><div>      &END OUTER_SCF</div><div>    &END SCF</div><div>    &XC                </div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>      &VDW_POTENTIAL</div><div><span class="Apple-tab-span" style="white-space:pre">  </span>DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>&PAIR_POTENTIAL</div><div><span class="Apple-tab-span" style="white-space:pre">          </span>TYPE DFTD3(BJ)</div><div><span class="Apple-tab-span" style="white-space:pre">               </span>REFERENCE_FUNCTIONAL PBE</div><div><span class="Apple-tab-span" style="white-space:pre">             </span>CALCULATE_C9_TERM .FALSE.</div><div><span class="Apple-tab-span" style="white-space:pre">            </span>PARAMETER_FILE_NAME /dftd3.dat</div><div><span class="Apple-tab-span" style="white-space:pre">               </span>R_CUTOFF 15.0</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>&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 15.138 15.138 30.000</div><div>      PERIODIC XYZ </div><div>    &END CELL</div><div>    &TOPOLOGY</div><div>      COORD_FILE_FORMAT XYZ</div><div>      COORD_FILE_NAME  ./c.xyz</div><div>    &END</div><div>    &KIND C</div><div>      ELEMENT C</div><div>      BASIS_SET DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>    &END KIND </div><div>    &KIND Ti</div><div>      ELEMENT Ti</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q12</div><div>    &END KIND</div><div>    &KIND O</div><div>      ELEMENT O</div><div>      BASIS_SET DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q6</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>  &PRINT</div><div><span class="Apple-tab-span" style="white-space:pre">    </span>&TOTAL_NUMBERS ON</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>&END TOTAL_NUMBERS</div><div>  &END PRINT</div><div>&END FORCE_EVAL</div><div>&MOTION</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>&GEO_OPT</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>  OPTIMIZER BFGS</div><div><span class="Apple-tab-span" style="white-space:pre">    </span>  MAX_FORCE 4.500000E-004</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>&END</div><div><span class="Apple-tab-span" style="white-space:pre">     </span>&CONSTRAINT</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>  &FIXED_ATOMS</div><div><span class="Apple-tab-span" style="white-space:pre">          </span>LIST 0..95</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>  &END FIXED_ATOMS</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>&END CONSTRAINT</div><div>&END MOTION</div></div><div><br></div><div><br></div><div><br></div><div> </div></div>