question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

NAMD parser does not take temperature as argument

See original GitHub issue

The NAMD parser #7 (for FEP) does not take T as an argument and does not return u_nk in proper units of kT.

Issue Analytics

  • State:closed
  • Created 4 years ago
  • Comments:19 (9 by maintainers)

github_iconTop GitHub Comments

1reaction
mrshirtscommented, May 7, 2019

The BAR and MBAR operate on dimensionless quantities: kT must be in the same units as the energies. Otherwise beta*U becomes something different. . . If BOTH the energies and kT change units, then everything is valid, except of course the free energies will be in those different units.

Currently the software modules assume you are passing in the beta*U quantities (already nondimensionalized). So you multiply the results by whatever kT you divided the energies by before you passed them into BAR.

1reaction
orbeckstcommented, Apr 15, 2019

Thanks for confirming.

Currently, the temperature from the output file is not used, is it?

I don’t think that it would be correct to use the instantaneous temperature for kT – the Boltzmann factor comes from a Lagrange multiplier for system at constant T (i.e., constant average temperature) so this T should be the “target temperature” of the thermostat in the simulation. But maybe @mrshirts @davidlmobley want to correct me here.

Read more comments on GitHub >

github_iconTop Results From Across the Web

NAMD User's Guide
The NAMD User's Guide describes how to run and use the various features of the molecular dynamics program NAMD.
Read more >
Is there any way to take inputs for a option created using ...
OptionParser has been deprecated, you should use argparse.ArgumentParser instead. With ArgumentParser you'll be able to do add_argument('--temp' ...
Read more >
Syntax error in "args=parser.parse_args()"
I took on a project to create a temperature monitoring system using the influxdb and Grafana desktop. I used the example from circuits.dk....
Read more >
Three ways to parse arguments in your Python code
The first option is to use argparse which is a popular Python module specialized in command-line parsing. Another way is to read the...
Read more >
Expressions | Style Specification | Mapbox GL JS - Mapbox docs
This example uses the get operator to get the temperature value of each feature. ... However, some paint and layout properties do not...
Read more >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found