[gmx-developers] Problem with dihedral restraint in gromacs 4.6

Anton Feenstra k.a.feenstra at vu.nl
Mon Feb 11 16:17:01 CET 2013


Could you please move this discussion to the users list?
This type of questions should be asked there.

On 11/02/13 16:07, Landry Charlier wrote:
> Last question (I hope).
>
> How can I choose the force constant?
>
> Indeed for the version 4.5, the force constant is defined by the option
> dihre_fc in the mdp input file. After this force is multiply by kfac
> contained in the dihedral_sections in the topology file
>
> Now this option is obsolete. The force constant is defined by kfac in
> gromacs4.6 ? If yes, kfac has now a unit. In kJ.mol-1.rad-2 ?
>
> Thank you
>
>
> Quoting "Shirts, Michael (mrs5pt)" <mrs5pt at eservices.virginia.edu>:
>
>> Hi, all-
>>
>>> Indeed, if i check the mdout.mdp file generated during the grompp
>>> step, there is no dihre option as in the version 4.5.
>>> So how can I turn on/off the dihedral restraints without the dihre
>>> option?
>>
>> Using the topology, like you are doing now.
>>
>>> For the moment, I modify the topology file by put or erase the
>>> dihedral_restraints section. The dihre option was more useful.
>>
>> The goal of this change was to harmonize all of the free energy and A
>> <-> B
>> state options; the mix of types resulted in bugs and complexity in the
>> underlying code.  Now, all types of transformation use the same
>> syntax.  If
>> there are issues with the interface, please let us know!
>>
>> Best,
>> ~~~~~~~~~~~~
>> Michael Shirts
>> Assistant Professor
>> Department of Chemical Engineering
>> University of Virginia
>> michael.shirts at virginia.edu
>> (434)-243-1821
>>
>>
>>
>> --
>> gmx-developers mailing list
>> gmx-developers at gromacs.org
>> http://lists.gromacs.org/mailman/listinfo/gmx-developers
>> Please don't post (un)subscribe requests to the list. Use the
>> www interface or send it to gmx-developers-request at gromacs.org.
>>
>
>
>
>


-- 
Groetjes,

Anton
  _____________ _______________________________________________________
|             |                                                       |
|  _   _  ___,| K. Anton Feenstra                                     |
| / \ / \'| | | IBIVU/Bioinformatics - Vrije Universiteit Amsterdam   |
|(   |   )| | | De Boelelaan 1081 - 1081 HV Amsterdam - Netherlands   |
| \_/ \_/ | | | Tel +31 20 59 87783 - Fax +31 20 59 87653 - Room P136 |
|             | Feenstra at few.vu.nl - www.few.vu.nl/~feenstra/         |
|             | "You Will Be Surprised At What Resides In Your        |
|             | Inside" (Arrested Development)                        |
|_____________|_______________________________________________________|



More information about the gromacs.org_gmx-developers mailing list