[gmx-users] angle force constant in mdp file?

[David]

On Mon, 2005-06-13 at 12:08 -0700, David Mobley wrote:
> Dear Gromacs users and team,
> 
> Concerning restraints, does anyone know why there isn't an option in
> the mdp files currently for changing angle restraint force constants?
> For example, dihre-fc and disre_fc are valid options and change the
> multiplier of the force constant specified in the topology file, but
> there doesn't seem to be a corresponding parameter for the angle force
> constant.
> 
> Does such a parameter exist? If not, does anyone know why not and if
> this would be expected in the future?

It does not. Don't know why. Science isn't always rational.
Anyway you can probably emulate it's functionality by using in your
topology

#define MULTFACTOR 1000.0
#define MYFC(x) (MULTFACCTOR*x)

And then using rather than your normal force constants

[ my_function ]
1 2 3 4 1 MYFC(3)
2 3 4 5 1 MYFC(4)

etc.



> 
> I'm just interested in doing a bunch of runs with slightly different
> restraint forces, and it is (obviously) more annoying to have to go in
> and manually edit the topology file for each different combination of
> restraint forces than to be able to specify this in the input file.
> 
> Thanks,
> David
> _______________________________________________
> 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.
-- 
David.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,          75124 Uppsala, Sweden
phone:  46 18 471 4205          fax: 46 18 511 755
spoel at xray.bmc.uu.se    spoel at gromacs.org   http://xray.bmc.uu.se/~spoel
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++