|
From: | Ulf Schiller |
Subject: | Re: [ESPResSo-users] referring to stored configs separately |
Date: | Tue, 26 Mar 2013 10:15:31 +0100 |
User-agent: | Mozilla/5.0 (X11; Linux i686 on x86_64; rv:17.0) Gecko/20130107 Thunderbird/17.0.2 |
Hi all, On 03/26/2013 09:49 AM, Olaf Lenz wrote:
-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi Koen! On 03/26/2013 08:18 AM, Koen Nickmans wrote:1. When a lennard-jones potential is shifted so that the potential is equal at a cutoff radius, how does espresso handle the calculation of forces around the cutoff point? I cannot find where this is discussed.The standard LJ potential in ESPResSo does have a discontinuity in the force at the cutoff radius. We do not handle the case otherwise. To my knowledge, this is only important when you choose the cutoff radius very small or want to simulate at very high accuracy. If you really need this, it should not be very hard to extend the LJ potential accordingly. We would be interested in that code if you write it.
Small addendum: If I'm not mistaken, it is a good idea to truncate the potential at the minimum. Otherwise energy conservation might be affected (NVE ensemble), and the minimum of the truncated potential may not be a point of zero force (which in turn might effect some energy minimization protocols). If the potential is cut off at a different point, one can 'taper' the potential to make its derivatives go to zero as well.
Best wishes, Ulf -- Dr. Ulf D. Schiller Building 04.16, Room 3006 Institute of Complex Systems (ICS-2) Phone: +49 2461 61-6144 Forschungszentrum Jülich, Germany Fax: +49 2461 61-3180
smime.p7s
Description: S/MIME Cryptographic Signature
[Prev in Thread] | Current Thread | [Next in Thread] |