Hello,<div><br></div><div>I have been trying to recreate a VASP calculation of CuO in CP2K 7.1 (I'm working on compiling CP2K 2022.1) with poor results.  I've also seen this calculation using Quantum Espresso with good success too.</div><div><br></div><div>The calculation uses PBE+U and I'm using broken symmetry to set up the antiferromagnetic spins. </div><div><br></div><div>In the paper (https://pubs.acs.org/doi/10.1021/acs.jpcc.5b10431), they are able to get pretty good agreement with the experimental geometry (Table 1) even without the +U but regardless of my +U value, the optimized cell is terrible and the band gap is too low.</div><div><br></div><div>I have tried different U_MINUS_J values, 0 to 5, and 3 seems to be the "best".  The high values lead to instability. I've also tried without broken symmetry but results are inconclusive (I need to play with this more).</div><div><br></div><div>Does anyone have any suggestions on why the calculation might be giving odd/wrong results?</div><div><br></div><div>My input, xyz coordinates, and output are attached.</div><div><br></div><div>If more information on my end is needed, please don't hesitate to ask.</div><div><br></div><div>This has been troubling me for a month now so any help would be greatly appreciated! </div><div><br></div><div>Best regards,</div><div>Frank</div>

<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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/5271528c-4306-4345-a439-100d546e84ddn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/5271528c-4306-4345-a439-100d546e84ddn%40googlegroups.com</a>.<br />