How does xtb constrain multiple bonds during the MD process? #1305
Replies: 1 comment
-
I'm sorry, I submitted this question repeatedly. |
Beta Was this translation helpful? Give feedback.
0 replies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
-
I need to constrain the distances between multiple atoms during the MD process. The 'constrain' section in the manual provides examples for constraining a single bond length, bond angle, or dihedral angle. I attempted to apply constraints to multiple atoms and increased the constraint force. Below is my input file:
$constraint
force constant=2.0
distance: 2, 9, 1.45
distance: 5, 6, 1.45
distance: 60, 61, 1.45
distance: 57, 64, 1.45
$end
However, it seems that the code isn't working as intended. I don't see a dedicated section listing bond lengths in the output file, similar to what's shown in the manual. Moreover, during the MD trajectory, I've observed significant deviations of these bonds from their preset values. In other discussions, I've come across methods where people indirectly achieve the constraint of multiple bonds by invoking xtb through programs like Orca (http://bbs.keinsci.com/thread-39106-1-1.html).
How should I structure the input file to constrain the distances between multiple atoms in xtb?
Thanks
Beta Was this translation helpful? Give feedback.
All reactions