<div>Dear Cp2k community,</div><div><br /></div><div>I want to calculate the RESP charges for a benzene ring with 2 Bromium atoms. Attached is my input script resp.inp, which worked so far for small molecules. However, I noticed that the magnitude of the partial charge is quite low for Br (q_i=-0.099), so I read up online and found that you need to include Halogen corrections to properly capture Br in particular, see the related paper Grimme 2017 from the cp2k website. In cp2k this is done via &XTB section and the USE_HALOGEN_CORRECTION keyword, correct?</div><div><br />I tried using the following section:</div><div><br /></div><div>&XTB<br /> COULOMB_INTERACTION T<br /> TB3_INTERACTION F<br /> CHECK_ATOMIC_CHARGES F</div><div> DO_EWALD T<br /> USE_HALOGEN_CORRECTION T</div><div> &PARAMETER<br /> DISPERSION_PARAMETER_FILE dftd3.dat<br /><br /> PARAM_FILE_NAME xTB0_parameters<br /> &END PARAMETER</div><div>&END XTB</div><div><br /></div><div>But there are multiple points where I am failing:</div><div><br /></div><div>The most important point is that I need a xTB_parameters file, but the only file I could find in cp2k repository was the xTB0_parameters file which begins with $info and every time I include it in my script, it says that $info as variable is not defined. However, as far as I understand I need some kind of file for the corrections, or no? If so where can I get the correct file?</div><div><br /></div><div>Second, am I even using the use_halogen_correction keyword correctly? In the end I want more realistic RESP charges for my molecule, especially on Bromium.</div><div><br /></div><div>I am very unfamiliar with xtb, so probably this is a very basic question. If you have any other input, I appreciate it as well! I attach my molecule xyz file as well as every other file for testing.</div><div><br /></div><div>Thank you very much in advance!</div><div><br /></div><div>Sincerely, </div><div>Andreas Doell</div><div><br /></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 visit <a href="https://groups.google.com/d/msgid/cp2k/51bcff8d-699b-497a-bffb-28460c9fd2d5n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/51bcff8d-699b-497a-bffb-28460c9fd2d5n%40googlegroups.com</a>.<br />