[lammps-users] Question about bond style

Dear all,

I would like to know whether there is a bond style in the format:
f=C__*r__ij
this is the bond style used in normal dpd simulation actually.

I know there is a most similar bond style
f=C*(r-ro)^2
But they are not the same.

Is it the only way to develop it of my own or I can realise this by other means?
Thanks~

Dear all,

I would like to know whether there is a bond style in the format:
f=C*rij
this is the bond style used in normal dpd simulation actually.

I know there is a most similar bond style
f=C*(r-ro)^2
But they are not the same.

are you sure you are not mixing up potential and forces?

axel.

There is no bond style in LAMMPS that implements
that formula, but you could add one. See doc/Section_modify.html
and use bond_harmonic.cpp as a starting point.

Steve

Actually, I did mixed them up!