Skip to content

Question about equilibration protocol  #29

Description

@jasmin-guven

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!

Activity

  1. jmichel80 commented on Oct 31, 2024

    @jmichel80
    Contributor

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

  2. jmichel80 commented on Oct 31, 2024

    @jmichel80
    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.

  3. jasmin-guven commented on Oct 31, 2024

    @jasmin-guven
    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 😊

  4. annamherz commented on Oct 31, 2024

    @annamherz

    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.

  5. lohedges commented on Oct 31, 2024

    @lohedges
    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.

  6. jasmin-guven commented on Oct 31, 2024

    @jasmin-guven
    Author

    Thanks everyone!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions