<div>Hello,</div><div><br></div><div>I just posted also on DFTB+ issues (https://github.com/dftbplus/dftbplus/issues/1001), including the input script for DFTB+</div><div><br></div><div>Xavier<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, April 5, 2022 at 3:22:02 PM UTC-5 Xavier Bidault 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>Hello,</div><div><br></div><div>I run a variable-cell optimization of a high-energy molecular crystal, beta-HMX, with CP2K 9.1 and DFTB+ 21.2, both using GFN1-xTB, which includes D3(BJ) correction. To compare apples to apples, I optimized the same supercell, using the Gamma point. The initial (experimental) volume for the supercell is 2069.80 A3. The resulting volume is:</div><div>1950.83 A3 with DFTB+</div><div>1759.92 A3 with CP2K</div><div><br></div><div>In order to check if I did something wrong in the configuration files, I used DFTB + D3(0) with the default DFTB parameter file from CP2K for both CP2K and DFTB+, and I find:</div><div>1467.79 A3 with DFTB+<br></div><div>1467.14 A3 with CP2K</div><div>So, DFTB + D3(0) is consistent between CP2K and DFTB+, but not GFN1-xTB.<br></div><div><br></div><div>IMO, the dispersion correction of xTB is wrong somewhere in CP2K. Looking at the difference between CP2K and DFTB+, could it be that D3(BJ) correction is counted twice in CP2K? Or with some wrong units we don't have access to?</div><div><br></div><div>I attached the CP2K input file for the GFN1-xTB calculation.</div><div><br></div><div>Please, let me know.</div><div><br></div><div>Thank you,</div><div>Xavier<br></div><div><br></div><div><br></div></blockquote></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/136300fe-59dd-43f6-b345-081e2ab63d5fn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/136300fe-59dd-43f6-b345-081e2ab63d5fn%40googlegroups.com</a>.<br />