[gmx-users] extracting bond streching parameters
Markus O Kaukonen
mokaukon at cc.helsinki.fi
Thu Jun 20 10:35:32 CEST 2002
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) ---
More information about the gromacs.org_gmx-users
mailing list