I am interested in fixing the cartesian coordinates of a few atoms in my simulation box. In addition, I would like to compute the force exerted on these fixed atoms. I noticed you can easily do this with the fix_freeze command, but it is limited to granular simulations. Is there a straightforward way I can do this with atomistic simulations (i.e. atom_style full)?