Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

crest with gfnff for intermolecular complexes has problem with charges #5

Closed
robpollice opened this issue Jun 2, 2020 · 1 comment

Comments

@robpollice
Copy link

robpollice commented Jun 2, 2020

Describe the bug
Crest crashes when GFN-FF is used for intermolecular complexes together with specification of charge and/or multiplicity. I am using the pre-compiled version 2.10.1 together with the pre-compiled version 6.3.0 of xtb. Here are the files: xtb.zip

To Reproduce
Steps to reproduce the behaviour:

  1. XYZ with more than one molecule (see mol.xyz in the example).
  2. Run crest mol.xyz -chrg 0 -uhf 0 -gfnff -nci -mquick -opt loose and the job crashes.
  3. The outputs crest.out and xtb.out are provided in the ZIP file.
  4. Run crest mol.xyz -gfnff -nci -mquick -opt loose and everything runs fine (after purging the folder of files from the previous run).
@pprcht
Copy link
Contributor

pprcht commented Jun 5, 2020

It seems to be that GFN-FF calculations fail in xtb if there is a .CHRG file with a specified molecular charge of 0 because it will produce a EEQ charge constrain error. This should be fixed with the upcoming xtb 6.3.1 version.

For now please just refrain from using -chrg 0 as a command line argument. The molecular charge of 0 is the default assumption and does not need specification.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants