I am doing an MD simulation in LAMMPS of molecules, where I need to compute the polarizability of the atoms at every time step and on the fly. It would be really helpful if anyone could provide a lead on this front. Please let me know if you need any further details.
Thank you.