[gmx-developers] calculating LJ (SR)
spoel at xray.bmc.uu.se
Mon Feb 16 06:51:10 CET 2004
On Mon, 2004-02-16 at 07:17, Ranjit Sudhakar Ranbhor wrote:
> I am calculating LJ potential (SR ie 1-5 and onwards) in one small
> program. While doing so I am using C6 and C12 parameters given under [
> nonbond_params ] in file ffG43a1nb.itp. I think under nonbond_params all
> combinations are not present. I am getting less energy (85% of Gromacs)
> than Gromacs calculated energy. Is there any rule for combination or I can
> simply make all combinations for C6 and C12?
First make a matrix from what's under [ atomtypes ] then overwrite the
matrix elements with [ nonbond_params ]. The procedure is coded
somewhere in grompp:
Alternatively you can just read a tpr file using gromacs routines.
> Ranjit Ranbhor
> Research Scholar
> Indian Institute of Technology-Bombay
> Biotechnology Group
> School of Biosciences and Bioengineering
> Powai, Mumbai 400076
> Web: http://www.btc.iitb.ac.in/~ranjit
> gmx-developers mailing list
> gmx-developers at gromacs.org
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-developers-request at gromacs.org.
David van der Spoel, PhD, Assist. 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
More information about the gromacs.org_gmx-developers