Dear CP2K users, <div><br></div><div>I am trying the reproduce an old paper that calculates de binding energies of oxygen to a heme molecule. According to this reference, different multiplicities for oxygen and heme should be taken into account. I've tried to calculate the different energies for the different multiplicity states, but I consistently get the same results no matter the value I set to the flag MULTIPLICITY. </div><div><br></div><div>I am sure I am overlooking something but I cannot figure it out... I would be very grateful if you could help me fix this. </div><div><br></div><div>I have copied the input file at the end of this email. </div><div><br></div><div>I am using CP2K5.1 . </div><div><br></div><div>Thank you in advance:</div><div><br></div><div>Salome</div><div><br></div><div><p><b><font size="1">&GLOBAL</font></b></p>
<p><b><font size="1">  PROJECT heme_pbe3</font></b></p>
<p><b><font size="1">  RUN_TYPE ENERGY</font></b></p>
<p><b><font size="1">  PRINT_LEVEL LOW</font></b></p>
<p><b><font size="1">&END GLOBAL</font></b></p>
<p><b><font size="1"><br></font></b></p><p><b><font size="1">&FORCE_EVAL<br></font></b></p>
<p><b><font size="1">  METHOD QS</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">  &SUBSYS</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">    &CELL</font></b></p>
<p><b><font size="1">      ABC 20 20 20</font></b></p>
<p><b><font size="1">      PERIODIC NONE</font></b></p>
<p><b><font size="1">    &END CELL</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">    &TOPOLOGY                    ! Section used to center the atomic coordinates in the given box. Useful for big molecules</font></b></p>
<p><b><font size="1">      &CENTER_COORDINATES</font></b></p>
<p><b><font size="1">      &END</font></b></p>
<p><b><font size="1">      COORD_FILE_FORMAT xyz</font></b></p>
<p><b><font size="1">      COORD_FILE_NAME  ./heme.xyz</font></b></p>
<p><b><font size="1">    &END</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">    &KIND FE</font></b></p>
<p><b><font size="1">      BASIS_SET DZVP-MOLOPT-SR-GTH</font></b></p>
<p><b><font size="1">      POTENTIAL GTH-PBE-q16</font></b></p>
<p><b><font size="1">    &END KIND</font></b></p>
<p><b><font size="1">    &KIND N</font></b></p>
<p><b><font size="1">      BASIS_SET DZVP-MOLOPT-SR-GTH</font></b></p>
<p><b><font size="1">      POTENTIAL GTH-PBE-q5</font></b></p>
<p><b><font size="1">    &END KIND</font></b></p>
<p><b><font size="1">    &KIND C</font></b></p>
<p><b><font size="1">      BASIS_SET DZVP-MOLOPT-SR-GTH</font></b></p>
<p><b><font size="1">      POTENTIAL GTH-PBE-q4</font></b></p>
<p><b><font size="1">    &END KIND</font></b></p>
<p><b><font size="1">    &KIND H</font></b></p>
<p><b><font size="1">      BASIS_SET DZVP-MOLOPT-SR-GTH</font></b></p>
<p><b><font size="1">      POTENTIAL GTH-PBE-q1</font></b></p>
<p><b><font size="1">    &END KIND</font></b></p>
<p><b><font size="1">  &END SUBSYS</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">  &DFT</font></b></p>
<p><b><font size="1">    CHARGE 0</font></b></p>
<p><b><font size="1">    MULTIPLICITY 3</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">    BASIS_SET_FILE_NAME BASIS_MOLOPT</font></b></p>
<p><b><font size="1">    POTENTIAL_FILE_NAME GTH_POTENTIALS</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">    &POISSON</font></b></p>
<p><b><font size="1">      PERIODIC NONE</font></b></p>
<p><b><font size="1">      POISSON_SOLVER WAVELET</font></b></p>
<p><b><font size="1">    &END POISSON</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">    &QS</font></b></p>
<p><b><font size="1">      EPS_DEFAULT 1.0E-10</font></b></p>
<p><b><font size="1">      METHOD GPW</font></b></p>
<p><b><font size="1">    &END QS</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">    &MGRID</font></b></p>
<p><b><font size="1">      CUTOFF 500</font></b></p>
<p><b><font size="1">      REL_CUTOFF 80</font></b></p>
<p><b><font size="1">      NGRIDS 5</font></b></p>
<p><b><font size="1">    &END MGRID</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">    &SCF</font></b></p>
<p><b><font size="1">      SCF_GUESS ATOMIC</font></b></p>
<p><b><font size="1">      EPS_SCF 1.0E-05</font></b></p>
<p><b><font size="1">      MAX_SCF 100</font></b></p>
<p><b><font size="1"><br></font></b></p>
<p><b><font size="1">      &OT T</font></b></p>
<p><b><font size="1">        MINIMIZER  DIIS</font></b></p>
<p><b><font size="1">        PRECONDITIONER  FULL_ALL</font></b></p>
<p><b><font size="1">      &END OT</font></b></p><p><b><font size="1"><br></font></b></p><p><b><font size="1">      &OUTER_SCF T</font></b></p><p><b><font size="1">        EPS_SCF     1.000E-05</font></b></p><p><b><font size="1">        MAX_SCF  100</font></b></p><p><b><font size="1">         &END OUTER_SCF</font></b></p><p><b><font size="1">    &END SCF</font></b></p><p><b><font size="1"><br></font></b></p><p><b><font size="1">    &XC</font></b></p><p><b><font size="1">      &XC_FUNCTIONAL PBE</font></b></p><p><b><font size="1">      &END XC_FUNCTIONAL</font></b></p><p><b><font size="1">    &END XC</font></b></p><p><b><font size="1"><br></font></b></p><p><b><font size="1">  &END DFT</font></b></p><p><b><font size="1">&END FORCE_EVAL</font></b></p></div>