<div dir="ltr">I am running into similar issue with iron complexes with PBE functional. Someone please help.<br><br>On Sunday, October 23, 2016 at 10:51:17 AM UTC-4, J. Ye wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi, <div><br></div><div>I studied a Ni(OH)2 cluster, basically it is a Ni4 cluster. each Ni bond with 4 OH. Then there's two open site on top and bottom of Ni-4OH planar (Ni-4OH is square). I put two H2O molecule on the top and bottom of Ni (1.9 angstrom). The H2O move away from Ni (about 3 angstrom) after optimized by CP2K (PBE, D3), however, H2O stays at Ni (2.0 angstrom) after optimized by Gaussian( M06L). </div><div><br></div><div>Any CP2K expert know the reason about this different results? </div><div><br></div><div>Thanks very much!</div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT Ni4_h2o_cluster</div><div>  RUN_TYPE GEO_OPT</div><div>  PRINT_LEVEL LOW</div><div>&END GLOBAL</div><div>&FORCE_EVAL</div><div>  METHOD QS</div><div>  &SUBSYS</div><div>    &CELL</div><div>      A  40.0000   0.0000   0.0000</div><div>      B  0.00000  40.0000  0.0000</div><div>      C  0.00000  0.00000  40.0000</div><div>      PERIODIC NONE</div><div>    &END CELL</div><div>     &TOPOLOGY</div><div>      COORD_FILE_NAME ./Ni4_h2o_cluster.xyz</div><div>      COORD_FILE_FORMAT XYZ</div><div>      CONNECTIVITY OFF</div><div>    &END TOPOLOGY</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>    &KIND H</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH-q1</div><div>      POTENTIAL GTH-PBE-q1</div><div>    &END KIND</div><div>    &KIND O</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH-q6</div><div>      POTENTIAL GTH-PBE-q6</div><div>    &END KIND</div><div>    &KIND Zr</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH-q12</div><div>      POTENTIAL GTH-PBE-q12</div><div>    &END KIND</div><div>    &KIND Ni</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH-q18</div><div>      POTENTIAL GTH-PBE-q18</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME ~/input/cp2k/cp2k-2.4.0/tests/<wbr>QS/BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME ~/input/cp2k/cp2k-2.4.0/tests/<wbr>QS/GTH_POTENTIALS</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-14</div><div>    &END QS</div><div>    &MGRID</div><div>      CUTOFF 360</div><div>      NGRIDS 4</div><div>      REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &SCF</div><div>      SCF_GUESS RESTART</div><div>      EPS_SCF 1.0E-07</div><div>      MAX_SCF 1000</div><div>      &OT TRUE</div><div>       PRECONDITIONER FULL_ALL</div><div>       MINIMIZER CG</div><div>       ENERGY_GAP 0.001</div><div>     &END OT</div><div>      &PRINT</div><div>        &RESTART OFF</div><div>        &END RESTART</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>           DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div>             &PAIR_POTENTIAL</div><div>              TYPE DFTD3</div><div>              REFERENCE_FUNCTIONAL PBE</div><div>              PARAMETER_FILE_NAME ~/input/cp2k/cp2k-2.4.0/tests/<wbr>QS/dftd3.dat</div><div>              R_CUTOFF 10.0</div><div>             &END PAIR_POTENTIAL</div><div>          &END VDW_POTENTIAL</div><div>    &END XC</div><div>  &END DFT</div><div>&END FORCE_EVAL</div><div>&MOTION</div><div>  &GEO_OPT</div><div>    TYPE MINIMIZATION</div><div>    MAX_DR    1.0E-03</div><div>    MAX_FORCE 1.0E-04</div><div>    RMS_DR    1.0E-03</div><div>    RMS_FORCE 1.0E-04</div><div>    MAX_ITER 2000</div><div>    OPTIMIZER BFGS</div><div>   &END GEO_OPT</div><div>  &PRINT</div><div>     &CELL HIGH</div><div>     &END CELL</div><div>     &STRUCTURE_DATA ON</div><div>     &END STRUCTURE_DATA</div><div>   &END PRINT</div><div>&END MOTION</div></div><div><br></div></div></blockquote></div>