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

Gulp optimisation and MD can produce odd H atom positions when starting from high energy structures. #162

Open
andrewtarzia opened this issue Jul 27, 2022 · 0 comments

Comments

@andrewtarzia
Copy link
Member

I have found many examples of structures after UFF4MOF optimisation or high T MD (using Gulp through stko), where H atoms near binding sites will be pointing the wrong way relative to their C atom. This issue normally resolves itself with subsequent GFN2-xTB optimisation or with higher T MD (not always the case for the MD, which suggests this is a local minima for the H atom).

Unsure the exact cause, but considering the initial "high energy" nature of stk molecules, with long bonds and bad inter-building block angles (even when using MCHammer of Collapser), this may not be so surprising and users should be careful.

As I do not know the cause, I assume it is due to high energy initial structures, I also do not know the fix from an stko perspective. Maybe it is the way the force field is described or the effect of typing (a Gulp issue).

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

1 participant