Re: AMBER: About Lennard-Jones Parameters

From: David A. Case <>
Date: Fri, 15 Sep 2006 17:07:35 -0700

On Fri, Sep 15, 2006, Jianhui wrote:
> In parm99.dat file, listed are the van der waals radius sigma and the well
> depth epsilon of different elements, which is used for the same element
> Lennard-Jones interactions such as O-O. How do AMBER calculate the
> Lennard-Jones parameters for different element interaction such as Na-O?

This is actually not written down all that often, even in the force field
papers. Here is the relevant code (from

 * MathOpConvertNonBondToAC
 * Author: Christian Schafmeister (1991)
 * Convert a pair of nonbond parameters dE1, dR1 and dE2, dR2 to
 * a set of A and C for a non-bond interaction.
MathOpConvertNonBondToAC( double dE1, double dR1, double dE2, double dR2,
        double *dPA, double *dPC )
double dE, dR, dR6, dER6;

    dR = dR1 + dR2;
    dE = sqrt( dE1 * dE2 );
    dR6 = pow6(dR);
    dER6 = dE*dR6;
    *dPC = 2.0*dER6; /* C = 2*dE*dR^6 */
    *dPA = dER6*dR6; /* A = dE*dR^12 */

In words: the two R* values are added together; the epsilon is the geometric
mean of the two epsilon values. Then A anc C are computed from the formulas

...hope this helps....dac

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Sun Sep 17 2006 - 06:07:25 PDT
Custom Search