<br>Thank you very much for your response. I will definitely look into k-point grids and apply your recommendations. See if that changes things.<br><br>Most of the flags/values I used were from my supervisor's (other) student, who recommended these to me. I'm still very new to CP2K, so I am still learning what to use/not to use.<div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, January 18, 2021 at 3:33:21 PM UTC+2 fa...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hi,</div><div><br></div><div>you are running the optimization without a k-point grid, so you only consider the gamma point. To get a more reasonable result add<br></div><div>    &KPOINTS<br>      SCHEME MONKHORST-PACK 20 20 20<br>    &END KPOINTS<br></div><div>to the &DFT section. Adjust the number of kpoints to your needs. When using k-points you should probably change METHOD DIRECT_P_MIXING to METHOD BOYDEN_MIXING.</div><div><br></div><div>A few more hints:</div><div>with DIAGONALIZATION you dont need the outer_scf section and ADDED_MOS 100 seems excessive.</div><div>@SET SCF_CONV 1E-5 is rather inaccurat, I'd recomment decreasing that to 1e-6 or 1e-7 even.</div><div>Is there a specific reason to decrease @SET GRIDS 3 from the default? That should not be needed.</div><div>When performing simulations with a variable cell size it's recommended to include a &CELL_REF section. Set these lattice vectors longer than the maximum value that you expect during the simulation (for your real cell).</div><div><br></div><div>Cheers,</div><div>Fabian<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, 18 January 2021 at 14:16:16 UTC+1 ASSIDUO Network wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi there, I am trying to calculate the lattice constant of bulk Au using CP2K (cell optimization is fine). I have included my input file below.<br><br>With CP2K, I get a lattice constant around 4.6 Å, while with VASP I get a value of 4.17 Å. The literature lattice constant for bulk Au should be around 4.08 Å. Why do I get such a massive difference in value between CP2K and VASP/literature?<br><br>What would you advise I do with my calculation to improve my results? The largest basis set available for Au (with CP2K) is TZV2P. My supervisor wants me to use the GTH-PBE potential, so I can't change that.<br><br>Thank you for your help in advance :)</blockquote></div></blockquote></div>