Hi All, <div><br></div><div>I was wondering if anyone has tips for computing bulk metallic system in CP2K (specifically for bulk Cu and bulk Ni system). I need the electronic energy for bulk Cu and bulk Ni for some modeling on NP thermodynamic stability. </div><div><br></div><div>I've optimized the unit cell parameters for both Cu and Ni systems by varying the unit cell dimensions found on <a href="https://materialsproject.org/materials/mp-30/">Materials Project</a> with the same level of theory as my CP2K on metal-organic frameworks. When I try to compute different metallic properties though such as atomization and bulk Cu + O2 (g) --> bulk CuO my energies are way off.. </div><div><p dir="ltr" style=""><font color="#000000" style="">From <a href="https://www.webelements.com/copper/thermochemistry.html" style=""><span style="text-decoration-line: underline;">Web Elements for Ni</span></a>, the atomization energy was determined:</font><br></p>bulk Cu (s) --> Cu(0) (g) w/ deltaH=338 kJ/mol</div><div><br>Now using DFT, the follow atomisation energy was determined to be:<br><br></div><div>bulk Cu (s) --> Cu(0)  (g) w/ deltaE=1658.2 kJ/mol  </div><div><br></div><div>Either my bulk system is too stable, or my atomized Cu is too unstable. Here is my DFT section for the bulk Cu calculations. </div><div><br></div><div><div>   &DFT</div><div>     BASIS_SET_FILE_NAME BASIS_file</div><div>     POTENTIAL_FILE_NAME POTENTIALS_file</div><div>     UKS  T </div><div>     MULTIPLICITY  2 </div><div>     CHARGE  0</div><div>     &SCF</div><div>       MAX_SCF  1000</div><div>       EPS_SCF     9.9999999999999995E-07</div><div>       SCF_GUESS  ATOMIC</div><div>       &OT  T</div><div>         MINIMIZER  CG</div><div>         PRECONDITIONER  FULL_ALL</div><div>         ENERGY_GAP     1.0000000000000000E-03</div><div>       &END OT</div><div>       &OUTER_SCF  T</div><div>         EPS_SCF     9.9999999999999995E-07</div><div>         MAX_SCF  50</div><div>       &END OUTER_SCF</div><div>     &END SCF</div><div>     &QS</div><div>       EPS_DEFAULT     1.0000000000000000E-10</div><div>       METHOD  GPW</div><div>     &END QS</div><div>     &MGRID</div><div>       NGRIDS  5</div><div>       CUTOFF     3.6000000000000000E+02</div><div>       REL_CUTOFF     8.0000000000000000E+01</div><div>     &END MGRID</div><div>     &XC</div><div>       DENSITY_CUTOFF     1.0000000000000000E-10</div><div>       GRADIENT_CUTOFF     1.0000000000000000E-10</div><div>       TAU_CUTOFF     1.0000000000000000E-10</div><div>       &XC_FUNCTIONAL  NO_SHORTCUT</div><div>         &PBE  T</div><div>         &END PBE</div><div>       &END XC_FUNCTIONAL</div><div>       &VDW_POTENTIAL</div><div>         POTENTIAL_TYPE  PAIR_POTENTIAL</div><div>         &PAIR_POTENTIAL</div><div>           TYPE  DFTD3(BJ)</div><div>           PARAMETER_FILE_NAME dftd3.dat</div><div>           REFERENCE_FUNCTIONAL PBE</div><div>           CALCULATE_C9_TERM  F</div><div>         &END PAIR_POTENTIAL</div><div>       &END VDW_POTENTIAL</div><div>     &END XC</div><div>   &END DFT</div></div><div><br></div><div>Sample .inp and .out files are attached. </div><div><br></div><div>Any feedback would be much appreciated! Thanks for taking a look at my issues. </div><div><br></div><div>With much appreciation, <br><br>Stephen </div>