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

Add torsion potential to GFN-FF for rotation around triple bonded carbon #787

Merged
merged 8 commits into from
Mar 24, 2023

Conversation

Thomas3R
Copy link
Contributor

This potential was requested to get the correct torsion potentials for diphenylacetylene.

Comment on lines 1870 to 1872
else
! allocate with size() = 0
allocate(topo%sTorsl(6, nn), source=0)
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do you really need this?
you can just remove this and then check in gfnf_eg, if the topo%sTorsl allocated

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I changed it below, thanks.

Account for the calc time of the new torsional term
Lower the distance for carbon triple bond
Do not allocate topo%sTorsl if no torsion around triple bonded carbon was found.
Only apply special torsion potential if topo%sTorsl is allocated.
@Thomas3R Thomas3R requested a review from Albkat March 23, 2023 17:27
@Thomas3R Thomas3R merged commit 5980893 into grimme-lab:main Mar 24, 2023
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

Successfully merging this pull request may close these issues.

2 participants