Hi,

I'm a phD student in the university of Strasbourg, France.

First of all, thank you very much for LAMMPS. I'm really enjoying using it.

I have a several question concerning 'angle style - harmonic':

1. Why after calculating the sin (s) you write s=1/s.

2. I don't understand the calculation of the forces. Does the matrice 'a' represent the dynamical matrice?

how did you calculate it?

3. I want to use this harmonic potential but to make a difference between positive and negative angles.

In order to calculate the energy, I think I can calculate sin(theta) with the vectorial product, and then take

theta=acos(c) for positive values of sin(theta) and theta=-acos(theta) for negative values. But, I don't know

what to do with the functions that calculates the forces? do I need to change them also?

thank you ,

Nava

Vous voulez savoir ce que vous pouvez faire avec le nouveau Windows Live ? Lancez-vous !