<div dir="ltr">Hi all, <div><br></div><div>I simulated the silicon with 64 atoms and calculated the energy-volume plot and fitted with murnaghan.</div><div>I found the lattice parameter 5.48 A ( lattice parameter in experiment is 5.44 A), That is OK but the bulk modulus is 63 GPa that is very low. (B_exp = 100 GPa)</div><div>I really don' t know why cp2k gets B so low???</div><div><br></div><div>I tried with DZVP basis set and PBE functional with 280 Ryd cut off for density.</div><div>I also checked with TZVP, it is not changing so much. It is 66 GPa.</div><div>Can anybody help me with that?<br></div><div>You can find my inputfile in below.</div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT si64-10.88A</div><div>  PRINT_LEVEL LOW</div><div>  RUN_TYPE GEO_OPT</div><div>&END GLOBAL</div><div>&MOTION</div><div>  &CELL_OPT</div><div>    TYPE GEO_OPT</div><div>    OPTIMIZER CG</div><div>    MAX_ITER 20</div><div>    EXTERNAL_PRESSURE [bar] 1.0</div><div>    MAX_DR 0.02</div><div>    RMS_DR 0.01</div><div>    MAX_FORCE 0.0002</div><div>    RMS_FORCE 0.0001</div><div>    KEEP_ANGLES</div><div>    &CG</div><div>      &LINE_SEARCH</div><div>        TYPE 2PNT</div><div>        &2PNT</div><div>        &END</div><div>      &END</div><div>    &END</div><div>  &END</div><div>  &GEO_OPT</div><div>    MAX_ITER 200</div><div>    MINIMIZER BFGS</div><div>  &END</div><div>  &PRINT</div><div>    &STRESS</div><div>    &END STRESS</div><div>    &RESTART</div><div>      &EACH</div><div>        CELL_OPT 1</div><div>      &END</div><div>    &END</div><div>  &END</div><div>&END</div><div>&FORCE_EVAL</div><div>  METHOD QS</div><div>  STRESS_TENSOR ANALYTICAL</div><div>  &DFT<br></div><div>    BASIS_SET_FILE_NAME ./BASIS_MOLOPT</div></div><div>    POTENTIAL_FILE_NAME ./GTH_POTENTIALS</div><div>    &MGRID</div><div>      CUTOFF 280</div><div>    &END MGRID</div><div>    &QS</div><div>      METHOD GPW</div><div>    &END QS</div><div>    &SCF</div><div>      SCF_GUESS RESTART</div><div>      MAX_SCF 20</div><div>      EPS_SCF 1.0e-5</div><div>      &OUTER_SCF</div><div>      EPS_SCF 1.0E-5</div><div>      &END</div><div>      &OT ON</div><div>#        MINIMIZER DIIS</div><div>         ROTATION</div><div>      &END OT</div><div>      &PRINT</div><div>        &RESTART</div><div>          BACKUP_COPIES 0</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>            CALCULATE_C9_TERM .TRUE.</div><div>            REFERENCE_C9_TERM .TRUE.</div><div>            LONG_RANGE_CORRECTION .TRUE.</div><div>            PARAMETER_FILE_NAME ./dftd3.dat</div><div>            REFERENCE_FUNCTIONAL PBE</div><div>            R_CUTOFF 8.</div><div>            EPS_CN 0.01</div><div>         &END PAIR_POTENTIAL</div><div>      &END vdW_POTENTIAL</div><div>    &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &TOPOLOGY</div><div>      COORD_FILE_NAME ./GEOMETRY.xyz</div><div>      COORDINATE XYZ</div><div>#      CONNECTIVITY OFF</div><div>    &END TOPOLOGY</div><div>    &KIND Si</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>#      BASIS_SET SZV-MOLOPT-SR-GTH</div><div>#       BASIS_SET TZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>    &END</div><div>    &CELL</div><div>      ABC    10.88  10.88 10.88</div><div>    &END</div><div>  &END SUBSYS</div><div> &END</div></div>