Force constants

Dear Lammps Users,

Does anyone know how to calculate second-order and third-order (cubic) force constants?

Best regards,

Nima Pirouzmand

You might want to give a little more context.