[gmx-users] settle for other solvents

[Anton Feenstra]

Anton Feenstra wrote:
> 
> David van der Spoel wrote:
> >
> [...]
> > You define the settle plane to be O-Me-Me and have the S as a dummy.
> > Then you reshuffle the mass of the S over the three atoms until the
> > moment of inertia is the same... Slightly easier, but still not
> > one-to-one mapping I'm afraid...
[...]
> 
> What about dummy atoms *with* mass? I think Berk explained to me once
> why and how that is impossible...

I just had a thought: You assign mass to a dummy atom. Coordinates
come from the dummy construction (based on other atoms), and velocity
from the difference with the previous coordinates. Then, mass and 
(change in) velocity will give you the intertial forces which you
could re-distribute just as the 'normal' forces are. Does this make
any sense at all, or am I out of my league here?


-- 
Groetjes,

Anton
 _____________ _______________________________________________________
|             |                                                       |
|  _   _  ___,| K. Anton Feenstra                                     |
| / \ / \'| | | Dept. of Pharmacochem. - Vrije Universiteit Amsterdam |
|(   |   )| | | De Boelelaan 1083 - 1081 HV Amsterdam - Netherlands   |
| \_/ \_/ | | | Room P275  Tel: +31 20 44 47608  Fax: +31 20 44 47610 |
|             | Feenstra at chem.vu.nl - www.chem.vu.nl/~feenstra/       |
|             | "Stay Cool, This is a Robbery" (Pulp Fiction)         |
|_____________|_______________________________________________________|