[gmx-users] Multiple torsional parameters
feenstra at chem.vu.nl
Wed Aug 20 08:59:11 CEST 2003
Pedro Alexandre Lapido Loureiro wrote:
> I am implementing the lipid force field of Smondyrev & Berkowitz (J Comp Chem
> 20, 531 - 1999).
> The problem is they use two expressions for one dihedral, such as:
> LC2 LC2O 1 0.0 5.8576 3
> LC2 LC2O 1 180.0 4.184 1
> I have built my .top file in this way.
> When I run grompp, there is a warning saying the former parameters will be
> overridden and the latter will be used.
> How can I solve this?
> (Sorry if a matter like this has been already discussed in the list...)
You are probably defining two dihedral *types* for the same pair of atoms.
That won't work. What you will need to do, is add another dihedral to the
.rtp, .itp or .top of your lipid, and specify explicit parameters (or an
explicit type, like in the G43a* ff's) for both of them.
| | |
| _ _ ___,| 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) |
More information about the gromacs.org_gmx-users