Skip to content

Apparent convergence issues #888

@amin-sagar

Description

@amin-sagar

Hello.

I tried doing a protein-peptide FEP (#682 )
I initially tried with the standard 11 replicas/lambda states.
mbar_overlap_matrix
replica_exchange_matrix
replica_state_timeseries

This seems to suggest not enough overlap between the replicas, so I increased the number to 15 and then 21 using

protocol_settings.simulation_settings.n_replicas = 21
protocol_settings.lambda_settings.lambda_windows = 21

Now, I get
mbar_overlap_matrix
replica_exchange_matrix
replica_state_timeseries

The timeseries plot seems to show exchanges but the transition matrix is mostly empty.
It doesn't seem like the transitions are 0 for the initial replicas.
Why could this be happening?
I would be really grateful for any suggestions.
Best,
Amin.

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