[gmx-users] criticism invited...

Paul Barrett barrett at biop.ox.ac.uk
Mon Mar 24 09:50:01 CET 2003


It happens every time, normally within about 300ps of restraint
removal.

Thanks to everyone for a great set of suggestions.

Paul

>Christoph Freudenberger wrote:
>> Paul Barrett wrote:
>> 
>>> My (alledged) solution: combine this hydrogen with its neighbouring
>>> oxygen by:
>[...]
>> I cannot judge how much this will effect the results in your special case,
>> but i don't think it's "good force field practice".
>
>I would tend to agree with that.
>
>> I would try to add the hydrogens mass into the oxygen
>> and define the hydrogen as a dummy keeping it's charge
>> and LJ's.
>
>Can't have a rotatable dummy, though. Try the addition of
>a dihedral angle potential David suggested, and/or shift
>mass from the oxygen to the hydrogen.
>
>Do you get this spinning OH group repeatedly? If it only
>happens eveery once in a while, you can use my 'upstream
>tweaking': write a .trr frame in slightly reduced precision
>(e.g. 6 in stead of 7 decimal places), and continue the
>run with that. I;ve posted a more complete description on
>the list only several days back.
>
>
>-- 
>Groetjes,
>
>Anton
>  _____________ _______________________________________________________
>|             |                                                       |
>|  _   _  ___,| K. Anton Feenstra                                     |
>| / \ / \'| | | Dept. of Pharmacochem. - Vrije Universiteit Amsterdam |
>|(   |   )| | | De Boelelaan 1083 - 1081 HV Amsterdam - Netherlands   |
>| \_/ \_/ | | | Tel: +31 20 44 47608 - Fax: +31 20 44 47610           |
>|             | Feenstra at chem.vu.nl - www.chem.vu.nl/~feenstra/       |
>|             | "If You See Me Getting High, Knock Me Down"           |
>|             | (Red Hot Chili Peppers)                               |
>|_____________|_______________________________________________________|
>
>
>_______________________________________________
>gmx-users mailing list
>gmx-users at gromacs.org
>http://www.gromacs.org/mailman/listinfo/gmx-users
>Please don't post (un)subscribe requests to the list. Use the 
>www interface or send it to gmx-users-request at gromacs.org.
>--------------------------------------------------------------------------------
>Return-Path: <gmx-users-admin at gromacs.org>
>Received: from oxmail.ox.ac.uk by nutmeg.biop.ox.ac.uk (MX V4.2 AXP) with SMTP;
>          Mon, 24 Mar 2003 08:09:54 +0000
>Received: from hawk.theophys.kth.se ([130.237.25.19]) by oxmail.ox.ac.uk with
>          smtp (Exim 4.10) id 18xN1m-00064f-01 for barrett at biop.ox.ac.uk; Mon,
>          24 Mar 2003 08:09:54 +0000
>Received: (qmail 18430 invoked from network); 24 Mar 2003 08:09:09 -0000
>Received: from localhost (HELO hawk.theophys.kth.se) (127.0.0.1) by 0 with
>          SMTP; 24 Mar 2003 08:09:09 -0000
>Received: (qmail 18316 invoked from network); 24 Mar 2003 08:08:10 -0000
>Received: from gallium.chem.vu.nl (130.37.144.37) by 0 with SMTP; 24 Mar 2003
>          08:08:10 -0000
>Received: from chem.vu.nl (far27.chem.vu.nl [130.37.31.139]) by
>          gallium.chem.vu.nl (8.9.1a/8.9.1) with ESMTP id KAA19598 for
>          <gmx-users at gromacs.org>; Mon, 24 Mar 2003 10:07:48 +0200
>Message-ID: <3E7B44B8.6040900 at chem.vu.nl>



More information about the gromacs.org_gmx-users mailing list