[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