[gmx-users] extracting bond streching parameters

[Markus O Kaukonen]

Dear all,

I did 'pdb2gmx' using
> 1 Gromacs Forcefield with all hydrogens (proteins only)

after 'grompp' one gets binary topology file 'topol.tpr'
which one can look at by 'gmxdump -s topol.tpr'

So finally to the question:

How do the following lines in the output of 'gmxdump'
>functype[81]=BONDS,b0A= 1.00e-01,cbA= 3.74e+05,b0B=1.00e-01,cbB= 3.74e+05
transform to bond streching parameters
k^b_ij and b_ij for each atom type pair ij
(ie. to parameters in eq. 4.33-4.34 in pdf manual ver. 3.11)?

Have a nice midsommer!, Markus


-- 
-----------www= http://www.iki.fi/markus.kaukonen---------------
Markus.Kaukonen at helsinki.fi   \ HOME->   Verajakallionkatu 6 a 5
Helsinki Bioenergetics Group   \_______  02600 Espoo, Finland
Institute of Biotechnology             \ tel:+358-(0)9-5127122/h
Biocenter 3, P.O.BOX 65 (Viikinkaari 1) \____
FIN-00014 University of Helsinki              \______ H+ ______
tel: +358-(0)9-19159752/o, fax: +358-(0)9-19159920             \
--- Rikos ei kannata, eika maatalous Suomessa. (Paimio 1998) ---