[gmx-users] Potential energy reported by gmx energy command

atb files atbfiles at zoho.com
Wed Jul 4 11:20:19 CEST 2018


Thanks Prof. David, now my energies are making some sense and they are following the trend that, with increase in particle number RMSD decreases by (N)^1/2 for LJ fluid. If -nmol option is available, which may be used to calculate per particle averages, why default unit of "KJ/mol" has been used, which now looks very confusing. Please help Yogi ---- On Fri, 29 Jun 2018 01:31:47 +0530 David van der Spoel <spoel at xray.bmc.uu.se> wrote ---- Den 2018-06-28 kl. 09:03, skrev atb files: > Hello experts, I did a LJ particles simulation, where I simulated two systems, 100 particle and 400 particle. For 100 particles energy is: Energy                      Average   Err.Est.       RMSD  Tot-Drift ------------------------------------------------------------------------------- Potential                  0.126288      0.029   0.615731 -0.0324041  (kJ/mol) Kinetic En.                 337.039       0.24    27.6568   0.443176  (kJ/mol) Total Energy                337.165       0.22    27.6615   0.410771  (kJ/mol) For 400 particle energy is: Energy                      Average   Err.Est.       RMSD  Tot-Drift ------------------------------------------------------------------------------- Potential                  0.502517      0.042    1.34584   0.220501  (kJ/mol) Kinetic En.                    1359       0.82    56.1468   -1.24836  (kJ/mol) Total Energy                 1359.5       0.78    56.1548   -1.02787  (kJ/mol) Since both systems have exact same operating conditions, only particle number has been changed, the per particle energy should be same, right? Here I see the energies has become 4 times in the second case. Then what is meaning of reporting energies in "KJ/mol" ? There is nothing like my system is not fully equilibrated, I have thoroughly checked, the both systems are well equilibrated. Please help in understanding these energies. Thank you > gmx energy does not know what your system is composed of. However the -nmol 400 flag may help. -- David van der Spoel, Ph.D., Professor of Biology Head of Department, Cell & Molecular Biology, Uppsala University. Box 596, SE-75124 Uppsala, Sweden. Phone: +46184714205. http://www.icm.uu.se -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.


More information about the gromacs.org_gmx-users mailing list