<div dir="ltr">Hi folks,<div><br><div>I'm new to CP2K and I want to measure adsorption energy of Ni atom (in further steps, cluster) on some oxide surfaces. Here are my questions:</div><div><br></div><div>1. When I try to calculate Energy of single Ni atom in gas phase by using Diagonalization method, it works with MULTIPLICITY 1, but once I set the MULTIPLICITY to 3, it doesn't converge. (Please check my input file below). I was wondering what's the reason for that? and how can I overcome this problem.</div><div><br></div><div>2. Is that reliable to calculate the energy of Ni atom/cluster with Diagonalization but the energy of complex (Ni on oxide surface) with OT and then compare them together to achieve the adsorption energy of metal? Do I need to use same method for all of the calculations? even if there is just one single Ni atom adsorbed on the surface? </div><div><br></div><div>3. the same question for POISSON_SOLVER: is that OK to use different Poisson solvers for metal atom/cluster and the complex? or again I need to use the same solvers to achieve the correct amount for adsorption energy?</div></div><div><br></div><div>I would highly appreciate your advises.</div><div><br></div><div>P.S:Here is my input file for the calculation of single Ni atom:</div><div><div>&GLOBAL</div><div>  PROJECT  m-Ni</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>    UKS</div><div>    MULTIPLICITY 3</div><div>    &MGRID</div><div><span style="white-space:pre">       </span>CUTOFF 800</div><div><span style="white-space:pre">    </span>REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &QS</div><div><span style="white-space:pre">  </span>METHOD GPW</div><div><span style="white-space:pre">    </span>EPS_DEFAULT 1.0E-12</div><div>    &END QS</div><div>    &POISSON</div><div><span style="white-space:pre">  </span>PERIODIC XY</div><div>        POISSON_SOLVER MT</div><div>    &END POISSON</div><div>     &SCF</div><div>     SCF_GUESS  ATOMIC</div><div>     EPS_SCF 1.0E-6</div><div>     MAX_SCF 500</div><div>     ADDED_MOS  20</div><div>     CHOLESKY INVERSE</div><div>     &SMEAR  ON</div><div>         METHOD FERMI_DIRAC</div><div>         ELECTRONIC_TEMPERATURE [K] 300</div><div>     &END SMEAR</div><div>     &DIAGONALIZATION</div><div>          ALGORITHM STANDARD</div><div>     &END DIAGONALIZATION</div><div>     &MIXING</div><div>          METHOD BROYDEN_MIXING</div><div>          ALPHA   0.1</div><div>          BETA    1.5</div><div>          NBROYDEN  8</div><div>     &END MIXING</div><div>     &PRINT</div><div><span style="white-space:pre">     </span>&RESTART</div><div><span style="white-space:pre">  </span> &EACH</div><div><span style="white-space:pre">            </span>QS_SCF 100</div><div><span style="white-space:pre">    </span> &END</div><div><span style="white-space:pre">     </span> ADD_LAST NUMERIC</div><div>        &END</div><div>     &END PRINT</div><div>    &END SCF</div><div>    &XC                       </div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>      &VDW_POTENTIAL</div><div><span style="white-space:pre">       </span>DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div><span style="white-space:pre">  </span>&PAIR_POTENTIAL</div><div><span style="white-space:pre">           </span>TYPE DFTD3(BJ)</div><div><span style="white-space:pre">                </span>REFERENCE_FUNCTIONAL PBE</div><div><span style="white-space:pre">              </span>CALCULATE_C9_TERM .FALSE.</div><div><span style="white-space:pre">             </span>PARAMETER_FILE_NAME ./dftd3.dat</div><div><span style="white-space:pre">               </span>R_CUTOFF 15.0</div><div><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>      A 11.35350 0.00000 0.00000</div><div>      B -5.67675 15.1390 0.00000</div><div>      C 0.000000 0.00000 30.0000</div><div>      PERIODIC XY</div><div>    &END CELL</div><div>    &COORD</div><div>        @INCLUDE 'Ni.coord'</div><div>    &END COORD</div><div>    &KIND C</div><div>      ELEMENT C</div><div>      BASIS_SET TZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>    &END KIND</div><div>    &KIND Ni</div><div>      ELEMENT Ni</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q18</div><div>    &END KIND</div><div>    &KIND O</div><div>      ELEMENT O</div><div>      BASIS_SET TZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q6</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>  &PRINT</div><div><span style="white-space:pre"> </span>&TOTAL_NUMBERS ON</div><div><span style="white-space:pre"> </span>&END TOTAL_NUMBERS</div><div>  &END PRINT</div><div>&END FORCE_EVAL</div><div>&MOTION</div><div><span style="white-space:pre">    </span>&GEO_OPT</div><div><span style="white-space:pre">  </span>  OPTIMIZER BFGS</div><div><span style="white-space:pre">     </span>  MAX_FORCE 4.500E-004</div><div><span style="white-space:pre">       </span>&END</div><div>&END MOTION</div></div><div><br></div></div>