From hinsen ":at:" ibs.ibs.fr Wed Nov 6 15:22:28 1996 Received: from ibs.ibs.fr for hinsen ^at^ ibs.ibs.fr by www.ccl.net (8.8.2/950822.1) id OAA12339; Wed, 6 Nov 1996 14:12:42 -0500 (EST) Received: from lmshp2.ibs.fr (lmshp2.ibs.fr [192.134.36.228]) by ibs.ibs.fr (8.6.12/8.6.12) with SMTP id UAA10578; Wed, 6 Nov 1996 20:13:26 +0100 Message-Id: <199611061913.UAA10578 /at\ibs.ibs.fr> Received: by lmshp2.ibs.fr (1.37.109.4/16.2) id AA12939; Wed, 6 Nov 96 20:13:26 +0100 Date: Wed, 6 Nov 96 20:13:26 +0100 From: Konrad Hinsen To: walt%!at!%panix.com Cc: echamot:~at~:xnet.com, rochus:~at~:felix.anorg.chemie.tu-muenchen.de, chemistry-!at!-www.ccl.net In-Reply-To: (message from Walter Polkosnik on Wed, 6 Nov 1996 13:38:51 -0500 (EST)) Subject: Re: CCL:MD/MM combination > There seems to be some confusion here with understanding of the > equipartition theorm. The classical equipartition theorem assigns a total > energy (note not just kinetic energy) of kT/2 per each term in the system's > Hamiltonian (momentum or position) that is quadratic. Only for ideal gases > is the total energy kinetic, for other molecules this is not true. Also, > non quadratic terms will contibute to the energy of the system but not > necessarily in the simple kT/2 way. For empirical potential energy I don't see the problem. The equipartition theorem says that there will be kT/2 of (kinetic) energy for each degree of freedom whose momentum enters quadratically into the Hamiltonian, and kT/2 of (potential) energy for each degree of freedom whose coordinate enters quadratically into the Hamiltonian. The first part remains true even if the second part can't be applied. As far as I know, the equipartition theorem is used in MD only for the momenta, not for the coordinates. The form of the potential energy function is therefore irrelevant. For those who want it precisely, here's the equipartition theorem for the kinetic energy in its exact form: / d H \ kT < p ------- > = ---- \ i d p / 2 i (Of course it's a partial derivative, but ASCII didn't provide for this!) -- ------------------------------------------------------------------------------- Konrad Hinsen | E-Mail: hinsen \\at// ibs.ibs.fr Laboratoire de Dynamique Moleculaire | Tel.: +33-76.88.99.28 Institut de Biologie Structurale | Fax: +33-76.88.54.94 41, av. des Martyrs | Deutsch/Esperanto/English/ 38027 Grenoble Cedex 1, France | Nederlands/Francais -------------------------------------------------------------------------------