How to write multi-linear harmonic potential in lammps input file

LAMMPS LAMMPS General Discussion
1

Hello all,

I would like to use this equation for bonded interaction of atoms.
F(r) =k1(r-ro) ; r<r1
=k1(r1-ro)+k2(r-r1) ; r2<r<r3
For the above equation, which bond_style should I mention in the Lammps input file?

Any suggestion would be highly appreciated.

Thank you.

2

There is no specific bond style for that kind of interaction, so you would have to use either bond style lepton or bond style table.