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

How to interpret wplots? #42

Open
DreamCykl opened this issue Jan 9, 2024 · 1 comment
Open

How to interpret wplots? #42

DreamCykl opened this issue Jan 9, 2024 · 1 comment

Comments

@DreamCykl
Copy link

Hi,

I've been using the pmx module for a while now to create single-topology merged molecules.
While I have a good understanding of most of the output, I havevn't found documentation explaining how the wplots should be interpreted.

I have two plots here, showing a transition from ligand A -> ligand B, in complex with protein.
image
image

Am I right in thinking that the overlap of the two histograms indicates good convergence? And more overlap reflects more overlap of the forward and reverse snapshots? With EQ FEP I believe you expect the forward and reverse energy values to converge towards the end of the time-series. Is this the same with non-EQ FEP? Or would good convergence / overlap of the forward/reverse work values be more present over the whole graph?

Thanks.

@vgapsys
Copy link
Member

vgapsys commented Jan 12, 2024

  1. Convergence looks fine. pmx analysis also gives a quantitative estimate for convergence
  2. Larger overlap of the work histograms corresponds to better convergence
  3. Two converged equilibrium FEP simulations at neighboring lambda windows should give similar work distributions to those that you have here
  4. I don't know what you mean with "present over the whole graph"

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

2 participants