<div dir="ltr">Dear Dr. Juerg,<div><br></div><div>Thank you very much.</div><div>The potential seems has been used by the BigDFT team for the delta parameter calculation on deltacodesdft website, with a small error of 2.7 meV/atom compared with Wien2k. According to the results of deltacodesdft, the potential may perform same well in different codes. So I guess maybe the basis is not suitable for bcc-Fe.</div><div>Is there any reference or tutorial you would like to suggest on optimize the basis set?</div><div><br></div><div>Best,</div><div><br></div><div>Yi Wang<br><br><blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi
<br>
<br>my guess would be 
<br>1) basis set
<br>2) pseudopotential
<br>
<br>regards
<br>
<br>Juerg
<br>------------------------------<wbr>------------------------------<wbr>--
<br>Juerg Hutter                        <wbr> Phone : ++41 44 635 4491
<br>Institut für Chemie C                FAX   : ++41 44 635 6838
<br>Universität Zürich                   E-<wbr>mail: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="1JXHdIBOBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">hut...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>------------------------------<wbr>------------------------------<wbr>---
<br>
<br>-----<a href="javascript:" target="_blank" gdf-obfuscated-mailto="1JXHdIBOBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a> wrote: -----To: cp2k <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="1JXHdIBOBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>>
<br>From: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="1JXHdIBOBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">wan...@gmail.com</a>
<br>Sent by: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="1JXHdIBOBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>
<br>Date: 05/16/2016 03:33PM
<br>Subject: [CP2K:7758] magnetization of bcc-Fe
<br>
<br>Hi everyone,
<br>I'm using the new cp2k 3.0 version with its k points sampling method to model the bcc-Fe.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.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)Is it rather a problem of cutoff and k points or the problem of potential and basis sets?I'm quite new to cp2k, wish this question is not too naive.
<br>Many thanks for any suggestion.
<br>
<br>
<br>&GLOBAL  PROJECT bcc-Fe  RUN_TYPE ENERGY_FORCE  PRINT_LEVEL LOW  &FM    TYPE_OF_MATRIX_MULTIPLICATION DBCSR_MM    FORCE_BLOCK_SIZE T  &END FM  &DBCSR    MM_DRIVER AUTO  &END DBCSR  FFTW_PLAN_TYPE MEASURE&END GLOBAL&FORCE_EVAL  METHOD Quickstep  &SUBSYS    &KIND Fe      ELEMENT   Fe      BASIS_SET DZVP-MOLOPT-SR-GTH      POTENTIAL GTH-PBE-q16    &END KIND    &CELL      A     -1.43 1.43 1.43      B     1.43 -1.43 1.43      C     1.43 1.43 -1.43      PERIODIC XYZ    &END CELL    &COORD      Fe    0.000000000    0.000000000    0.000000000    &END COORD  &END SUBSYS  &DFT    BASIS_SET_FILE_NAME  BASIS_MOLOPT    POTENTIAL_FILE_NAME  GTH_POTENTIALS    SPIN_POLARIZED       T    MULTIPLICITY 3    RELAX_MULTIPLICITY 0.    &PRINT      &HIRSHFELD ON         REFERENCE_CHARGE MULLIKEN      &END HIRSHFELD    &END PRINT    &QS      METHOD  GPW      EPS_DEFAULT 1.0E-14      EXTRAPOLATION USE_GUESS    &END QS    &MGRID      NGRIDS 5      CUTOFF 500      REL_CUTOFF 100    &END MGRID    &XC      &XC_FUNCTIONAL PBE      &END XC_FUNCTIONAL    &END XC    &SCF      SCF_GUESS ATOMIC      EPS_SCF 2.0E-7      MAX_SCF 300      &DIAGONALIZATION  ON        ALGORITHM STANDARD      &END DIAGONALIZATION      &MIXING  T        METHOD BROYDEN_MIXING        NSKIP        4        N_SIMPLE_MIX 4        ALPHA 0.08        BETA  0.25        NBUFFER 12      &END MIXING      ADDED_MOS 6 6      CHOLESKY  INVERSE_DBCSR      MAX_DIIS  6      &SMEAR T        METHOD FERMI_DIRAC        ELECTRONIC_TEMPERATURE [K] 300      &END SMEAR    &END SCF    &KPOINTS      SCHEME MONKHORST-PACK 12 12 12      FULL_GRID T    &END KPOINTS  &END DFT  &PRINT    &FORCES ON    &END FORCES  &END PRINT&END FORCE_EVAL
<br>
<br>
<br>Yi Wang--------------------------<wbr>-----Yi WangPh.D student, Nanjing University of Sci. & Tech.
<br>
<br>
<br>
<br>-- 
<br>
<br>You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>
<br>To unsubscribe from this group and stop receiving emails from it, send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="1JXHdIBOBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.
<br>
<br>To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="1JXHdIBOBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>
<br>Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.
<br>
<br>For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>
<br>
<br></blockquote></div></div>