Analyze compute the energy difference between the wrong states
See original GitHub issueThe function analyze.analyze()
extracts from the ncfile the number of states without taking into account the presence of the fully interacting state, and ends up computing the free energy difference between the fully interacting state and the next to last state for each phase here
# Accumulate free energy differences
entry = dict()
entry['DeltaF'] = Deltaf_ij[0,nstates-1]
entry['dDeltaF'] = dDeltaf_ij[0,nstates-1]
entry['DeltaH'] = DeltaH_i[0,nstates-1]
entry['dDeltaH'] = dDeltaH_i[0,nstates-1]
entry['DeltaF_restraints'] = DeltaF_restraints
data[phase] = entry
In general, since we treat the fully interacting state as special, we need to be careful during the analysis.
I’ll soon open a PR with the fix.
Issue Analytics
- State:
- Created 7 years ago
- Comments:5 (5 by maintainers)
Top Results From Across the Web
Guidelines for the analysis of free energy calculations
tool covers analysis of alchemical calculations generally, in- cluding free energy estimates ... To compute accurate free energy differences between states,.
Read more >Ionization Energy
Ionization energy is the quantity of energy that an isolated, gaseous atom in the ground electronic state must absorb to discharge an ...
Read more >Analysis of Situations in Which Mechanical Energy is ...
A common Physics lab involves the analysis of a pendulum in its back and forth motion. The transformation and conservation of mechanical energy...
Read more >Error Analysis
Every lab report must have an error analysis. For many experiments, significant figure rules are sufficient. For a brush up on significant figure...
Read more >Conservation of Energy Principle - an overview
The first law of thermodynamic defines internal energy; it is associated with all thermodynamic systems, and since heat is an energy form, it...
Read more >Top Related Medium Post
No results found
Top Related StackOverflow Question
No results found
Troubleshoot Live Code
Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start FreeTop Related Reddit Thread
No results found
Top Related Hackernoon Post
No results found
Top Related Tweet
No results found
Top Related Dev.to Post
No results found
Top Related Hashnode Post
No results found
Top GitHub Comments
Ah, yes, good catch. It may be better to just replace
[0,nstates-1]
with [0,-1] to make it independent of how much shifting we do at the head of the row.Closing this to open a better titled issue.