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

Question about equilibration protocol #29

Open
jasmin-guven opened this issue Oct 31, 2024 · 6 comments
Open

Question about equilibration protocol #29

jasmin-guven opened this issue Oct 31, 2024 · 6 comments

Comments

@jasmin-guven
Copy link

Hi!

I was just refreshing my memory of the BSS syntax for minimisation and equilibration, and noticed here that the tutorials first have a restrained NVT step, followed by a non-restrained NVT step, which then is followed by a restrained NPT step.

I wanted to ask if there is a reason why the position restraints are relaxed in between two position restraint steps? This seems different to similar tutorials from Amber, for example. I've also had some collaborators also point this out to me in the past, so I thought I'd ask here.

In my specific case of working with metals, the removal of restraints in between NVT and NPT allows the system to relax too much too early and can lead to problems with coordination.

Thanks a lot in advance!

@jmichel80
Copy link
Contributor

Hi @jasmin-guven I believe this protocol is based on @annamherz ‘s study . @annamherz can you comment ?

@jmichel80
Copy link
Contributor

Just to add that if your use case warrants it it’s simple to just update locally the equilibration protocol. The RBFE tutorial illustrates how a BSS RBFE pipeline may be deployed but there’s no guarantee that the protocol is optimal for any input dataset.

@jasmin-guven
Copy link
Author

Just to add that if your use case warrants it it’s simple to just update locally the equilibration protocol. The RBFE tutorial illustrates how a BSS RBFE pipeline may be deployed but there’s no guarantee that the protocol is optimal for any input dataset.

Yes I see what you mean, thanks Julien! I was just curious to see if there was a specific reason behind this protocol 😊

@annamherz
Copy link

Hello 😊 I think this is the protocol I used at the start (I adapted it slightly more after with the restraint weights) but maybe as I was using the benchmarking systems I never experienced any issues with it, so I did not spend a lot of time optimising it after my systems were equilibrated. I agree that now I would probably run it without the release of the restraints at the end of the NVT step as I think that's more robust.

@lohedges
Copy link
Contributor

Thanks for the clarification. We can leave this issue open as a reminder to rework the protocol for the next iteration of the tutorials.

@jasmin-guven
Copy link
Author

Thanks everyone!

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

4 participants