Hello all,<div><br></div><div>I performed two cell optimizations on bulk Si, the first with a primitive cell and kpoints and the second using OT and a 250 atom supercell. Both structures started from the Materials Project (mp-149) lowest energy Fm-3m polymorph.</div><div><br></div><div>I found that while OT only took 3 cell opt steps, the kpoint calculation took 10 steps, which I found a little strange.</div><div><br></div><div>Upon investigation I found two things: </div><div>(a) the kpoint calculation, had a large initial pressure deviation of -92570 bar, compared to OT, which had an initial deviation of  -708 bar at DZVP level and -2934 at TZV2PX level. Again these are the same structure, just kpoints vs supercell. It is a little odd to me that DZVP is closer to converged pressure, but I don't think its too strange.</div><div>(b) the kpoint calculation relaxed from Fm-3m to Im-3m space group, which is not the ground state polymorph.</div><div><br></div><div>My question is whether or not this issue is due to CP2K's kpoint implementation, or if this is something more fundamental, and how (if at all) it can be resolved. Below is the abridged DFT section for the kpoint calculation, which I assume is where the problem lies, but I can provide more info if needed.</div><div><br></div><div>Best,</div><div>Nick</div><div><p>&DFT</p>
<p>                BASIS_SET_FILE_NAME BASIS_MOLOPT</p><p>                BASIS_SET_FILE_NAME BASIS_MOLOPT_UCL</p>
<p>                POTENTIAL_FILE_NAME GTH_POTENTIALS</p>
<p>                UKS True</p>
<p>                MULTIPLICITY 0</p>
<p>                CHARGE 0</p>
<p>                &QS</p>
<p>                        METHOD GPW</p>
<p>                        EPS_DEFAULT 1e-12</p>
<p>                        EXTRAPOLATION PS</p>
<p>                &END QS</p>
<p>                &SCF</p>
<p>                        MAX_SCF 100</p>
<p>                        EPS_SCF 1e-07</p>
<p>                        SCF_GUESS RESTART</p>
<p>                        MAX_ITER_LUMO 400</p>
<p>                        ADDED_MOS 100</p>
<p>                        &DIAGONALIZATION</p>
<p>                                &DAVIDSON</p>
<p>                                        PRECONDITIONER FULL_ALL</p>
<p>                                &END DAVIDSON</p>
<p>                        &END DIAGONALIZATION</p>
<p>                        &MIXING</p>
<p>                                METHOD BROYDEN_MIXING</p>
<p>                                ALPHA 0.2</p>
<p>                                NBUFFER 5</p>
<p>                        &END MIXING</p>
<p>                        &SMEAR</p>
<p>                                ELEC_TEMP 300</p>
<p>                                METHOD FERMI_DIRAC</p>
<p>                                FIXED_MAGNETIC_MOMENT -100.0</p>
<p>                        &END SMEAR</p>
<p>                &END SCF</p>
<p>                &KPOINTS</p>
<p>                        SCHEME MONKHORST-PACK 13 13 13</p>
<p>                        EPS_GEO 1e-06</p>
<p>                        FULL_GRID False</p>
<p>                        PARALLEL_GROUP_SIZE -1</p>
<p>                        SYMMETRY False</p>
<p>                        UNITS B_VECTOR</p>
<p>                        VERBOSE True</p>
<p>                        WAVEFUNCTIONS COMPLEX</p>
<p>                &END KPOINTS</p><p>&END DFT</p></div>