From 6d61408106c32a8468f657465efb1f3fae3bd7b7 Mon Sep 17 00:00:00 2001 From: Michele Invernizzi Date: Sun, 15 Nov 2020 21:35:54 +0100 Subject: [PATCH] ENERGY doc: warning about tail corrections (#644) * added warning about long tail bug * adjust lammps command --- src/colvar/Energy.cpp | 11 ++++++++--- 1 file changed, 8 insertions(+), 3 deletions(-) diff --git a/src/colvar/Energy.cpp b/src/colvar/Energy.cpp index 0e03d35038..cc570a26e4 100644 --- a/src/colvar/Energy.cpp +++ b/src/colvar/Energy.cpp @@ -34,17 +34,22 @@ namespace colvar { /* Calculate the total potential energy of the simulation box. -Total energy can be biased with umbrella sampling \cite bart-karp98jpcb or with well tempered metadynamics \cite Bonomi:2009p17935. +The potential energy can be biased e.g. with umbrella sampling \cite bart-karp98jpcb or with well-tempered metadynamics \cite Bonomi:2009p17935. Notice that this CV could be unavailable with some MD code. When it is available, and when also replica exchange is available, metadynamics applied to ENERGY can be used to decrease the number of required replicas. +\bug This \ref ENERGY does not include long tail corrections. +Thus when using e.g. LAMMPS `"pair_modify tail yes"` or GROMACS `"DispCorr Ener"` (or `"DispCorr EnerPres"`), +the potential energy from \ref ENERGY will be slightly different form the one of the MD code. +You should still be able to use \ref ENERGY and then reweight your simulation with the correct MD energy value. + \bug Acceptance for replica exchange when \ref ENERGY is biased -is computed correctly only of all the replicas has the same +is computed correctly only if all the replicas have the same potential energy function. This is for instance not true when -using GROMACS with lambda replica exchange of with plumed-hrex branch. +using GROMACS with lambda replica exchange or with plumed-hrex branch. \par Examples