<div dir="ltr">Hi everyone,<div><br></div><div>I'm using the new cp2k 3.0 version with its k points sampling method to model the bcc-Fe.</div><div>I hope the magnetization per atom can be properly modeled, so that I may use it for large cells with point defect, where distribution of magnetization per atom will be explored.</div><div>I used the following input, which produces a magnetization of 3.6 muB/atom, which is larger than experiments and other simulations of around 2.2 muB/atom. (Similar in the case of fcc-Ni, which is 0.63 muB/atom > 0.52 muB/atom in exp)</div><div>Is it rather a problem of cutoff and k points or the problem of potential and basis sets?</div><div>I'm quite new to cp2k, wish this question is not too naive.</div><div><br></div><div>Many thanks for any suggestion.</div><div><br></div><div><br></div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT bcc-Fe</div><div>  RUN_TYPE ENERGY_FORCE</div><div>  PRINT_LEVEL LOW</div><div>  &FM</div><div>    TYPE_OF_MATRIX_MULTIPLICATION DBCSR_MM</div><div>    FORCE_BLOCK_SIZE T</div><div>  &END FM</div><div>  &DBCSR</div><div>    MM_DRIVER AUTO</div><div>  &END DBCSR</div><div>  FFTW_PLAN_TYPE MEASURE</div><div>&END GLOBAL</div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &SUBSYS</div><div>    &KIND Fe</div><div>      ELEMENT   Fe</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q16</div><div>    &END KIND</div><div>    &CELL</div><div>      A     -1.43 1.43 1.43</div><div>      B     1.43 -1.43 1.43</div><div>      C     1.43 1.43 -1.43</div><div>      PERIODIC XYZ</div><div>    &END CELL</div><div>    &COORD</div><div>      Fe    0.000000000    0.000000000    0.000000000</div><div>    &END COORD</div><div>  &END SUBSYS</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME  BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME  GTH_POTENTIALS</div><div>    SPIN_POLARIZED       T</div><div>    MULTIPLICITY 3</div><div>    RELAX_MULTIPLICITY 0.</div><div>    &PRINT</div><div>      &HIRSHFELD ON</div><div>         REFERENCE_CHARGE MULLIKEN</div><div>      &END HIRSHFELD</div><div>    &END PRINT</div><div>    &QS</div><div>      METHOD  GPW</div><div>      EPS_DEFAULT 1.0E-14</div><div>      EXTRAPOLATION USE_GUESS</div><div>    &END QS</div><div>    &MGRID</div><div>      NGRIDS 5</div><div>      CUTOFF 500</div><div>      REL_CUTOFF 100</div><div>    &END MGRID</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>    &END XC</div><div>    &SCF</div><div>      SCF_GUESS ATOMIC</div><div>      EPS_SCF 2.0E-7</div><div>      MAX_SCF 300</div><div>      &DIAGONALIZATION  ON</div><div>        ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>      &MIXING  T</div><div>        METHOD BROYDEN_MIXING</div><div>        NSKIP        4</div><div>        N_SIMPLE_MIX 4</div><div>        ALPHA 0.08</div><div>        BETA  0.25</div><div>        NBUFFER 12</div><div>      &END MIXING</div><div>      ADDED_MOS 6 6</div><div>      CHOLESKY  INVERSE_DBCSR</div><div>      MAX_DIIS  6</div><div>      &SMEAR T</div><div>        METHOD FERMI_DIRAC</div><div>        ELECTRONIC_TEMPERATURE [K] 300</div><div>      &END SMEAR</div><div>    &END SCF</div><div>    &KPOINTS</div><div>      SCHEME MONKHORST-PACK 12 12 12</div><div>      FULL_GRID T</div><div>    &END KPOINTS</div><div>  &END DFT</div><div>  &PRINT</div><div>    &FORCES ON</div><div>    &END FORCES</div><div>  &END PRINT</div><div>&END FORCE_EVAL</div></div><div><br></div><div><br></div><div><br></div><div>Yi Wang</div><div>-------------------------------</div><div>Yi Wang</div><div>Ph.D student, Nanjing University of Sci. & Tech.</div></div>