Dear colleagues,<br /><br />I am comparing the MoS2 model with the PBE xc-functional in CP2K and CPMD with the same GTH pseudopotential. I have found some strange discrepancies.<br /><br />For comparison, I calculate energy differences while changing the lattice constant in the xy-plane (compression and tension). In these tests I use Gamma point only. Please see 3 plots with the Energy-vs-Strain attached.<br /><br />1) First, let us consider the bulk MoS2 in the 2H structure. CP2K and CPMD with the same GTH pseudopotential give different results. Other CPMD calculations with the ‘CPMD-native’ GTH model and the CPMD MT pseudopotential give similar results to CPMD-withGTHcp2k. Interestingly, VASP results are very close to CP2K.<br /><br />In order to decipher the puzzle, I have made calculations for pure Mo and for pure S.<br /><br />2) For the bcc Mo structure all the codes with different pseudopotentials give very close results.<br /><br />3) However, for the S6 structure of sulphur, the CPMD-withGTHcp2k and CPMD-withMT results agree well with VASP, but do not match CP2K. And the CPMD results for its ‘native’ GTH model show a significant offset. <br /><br />I also tried changing energy cutoffs, enabling various k-meshes, and enlarging the cell, however I couldn’t find the possible reasons for these discrepancies in CP2K and CPMD with the same GTH pseudopotential of MoS2 and for S6. <br /><br />There are example files for these comparisons attached. <br /><br />Please let me know your best judgement about the possible roots of these discrepancies.<br /><br />Best regards,<br />Nikita<br />

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/79e3b838-3ca5-4f31-a606-dbf12dec50cen%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/79e3b838-3ca5-4f31-a606-dbf12dec50cen%40googlegroups.com</a>.<br />